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rom this fountain (the free will of god) it is those 
laws, which we call the laws of nature, have flowed, 
in which there appear many traces of the most wise 
contrivance, but not the least shadow of necessity. 
These therefore we must not seek from uncertain con- 
jectures, BUT LEARN THEM FROM OBSERVATIONS AND EXPER- 
IMENTAL. He who is presumptuous enough to think that 
he can find the true principles of physics and the laws 
of natural things by the force alone of his own mind, and the internal 
light of his reason, must either suppose the world exists by necessity, 
and by the same necessity follows the law proposed; or if the order of 
Nature was established by the will of God, the [man] himself, a miserable 
reptile, can tell what was fittest to be done. 



Isaac Newton 
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Summary 



Based on the assumption of the Einstein equivalence principle and the principle of gen- 
eral covariance general relativity describes the gravitational field successfully as a purely 
geometrical property of four dimensional spacetime on Riemannian manifolds. However, 
despite the so far remarkable accuracy in its experimental verification, general relativity 
remains a classical theory. So the necessity of finding a quantum mechanical description 
of the gravitational field implies the need to embed or to modify the above principles in 
a more general framework, which is one of the major challenges in modern theoretical 
physics. 

In this thesis we investigate specific predictions of such a class of more general the- 
ories, the so called nonmetric theories of gravity. Within this framework theories based 
on e.g., a metric-affine geometry of spacetime predict that in a gravitational field a pair 
of orthogonal linear polarisation states of light propagates with different phase veloc- 
ities. This gravity-induced birefringence could in principle be measured in local test 
experiments and, hence, violates the Einstein equivalence principle. Therefore we have 
used polarization measurements in solar spectral lines, as well as in continua and lines of 
various isolated magnetic white dwarfs and of cataclysmic variables (interacting binary 
systems) to constrain the essential coupling constants for this effect predicted by metric- 
affine gravity and other prototypes of nonmetric theories. These measurements provide 
an empirical formula which predicts the upper limit on the metric-affine coupling con- 
stant, measured for a particular celestial body as a function of its Schwarzschild radius 
and its physical, stellar radius. By modelling the lightcurves of a certain cataclysmic 
variable system, the results could, in principle, be interpreted as a direct detection of 
gravitational birefringence, although alternatives exist. 

This thesis provides the first systematic search for signals of gravitational birefringence 
in astronomical polarimetric data. As an outlook I propose further promising tests which 
also have the potential for setting strong upper limits on gravity-induced birefringence. 
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Chapter 1 
Introduction 



The Einstein equivalence principle plays the role of a key element in the development of 
new improved theories of gravity. Although being an important building block in Ein- 
stein's general relativity, theoretically predicted violations of its validity are an important 
feature in alternative, nonmetric gravitation theories if they are to incorporate quantum 
mechanical principles. Hence, the intention of this chapter is to motivate the conviction, 
grown within the last few years, that violations of the equivalence principle must be an 
essential part of every theory of gravity which pays attention to the quantum mechanical 
character of matter. 

After a brief historical outline of the weak and the Einstein equivalence principle and 
its implications, this chapter presents a theoretical framework which admits the analysis 
as well as the development of experimental tests for a broad class of gravitation theories. 
This purpose requires a critical examination of the underlying, mostly classical, concepts 
and notions. So to say as a side effect one is led to the possibility of looking at EEP 
violations as violations of spacetime symmetries in the spirit of modern quantum field 
theory. Indeed the principle of gauge symmetries, taken from the Standard model of 
elementary particle physics is used as a cornerstone in the mathematical formalism of 
the metric-afline gauge theory of gravity (MAG), which is the second theory considered 
here besides Moffat's nonsymmetric gravitation theory (NGT). Metric-affine theories can 
be regarded as extensions of Einstein-Cartan type theories. Becoming nonmetric when 
the additional gravitational potentials couple directly to matter, MAG as well as NGT 
predict that a gravitational field singles out an orthogonal pair of polarization states of 
light that propagate with different phase velocities. This gravity-induced birefringence 
implies that propagation through a gravitational field can alter the polarization of light 
and, so, violate the Einstein equivalence principle. Quantitative predictions for this phase 
shift are given which are used in the following chapters for setting strong limits on this 
effect by utilizing astrophysical spectropolarimetry of compact stellar objects. 

The polarization of an electromagnetic wave is completely and consistently described 
by a system of four real valued quantities, called Stokes parameters. Since the possible 
influence of gravity- induced birefringence on polarized light shows up in an alteration of 
these parameters, a brief introduction to this topic is given at the end of this chapter. 
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CHAPTER 1. INTRODUCTION 



1.1 Equivalence Principles 

The significance of the principles of equivalence for the development of modern physics 
can hardly be overestimated. For example, Galilei's famous free fall experiments per- 
formed from the leaning tower of Pisa marked the beginning of the development from the 
medivial, aristotelic way of science, up to modern physics pQ. Also, Newton realized very 
soon that his new ideas about the principles of motion and a universal gravitational force 
are basically founded on the equivalence between gravitational and inertial mass, so that 
he performed numerous pendulum experiments to have an experimental justification for 
his new laws. The importance this equivalence had for Newton can easily be estimated 
from the fact that he later devoted the opening paragraph of his Principia [2 to it. What 
Newton and also Galilei had introduced into modern physics is today known as the 



Weak Equivalence Principle (WEP): 

In a gravitational field all bodies fall with the same acceleration regardless of their mass 

or internal structure. 



The weak equivalence principle currently belongs to the physical predictions with the 
most accurate empirical underpinning. Beginning with the torsion balance experiments, 
performed by Baron von Eotvos and collaborators in the 19th. century an accuracy of 
approximately 10~ 9 is currently reached (in comparison to 10~ 3 of Newton's pendulum 
tests), while an accuracy of 10~ 12 — 10~ 15 is theoretically expected for free fall experiments 
in orbit, e.g. for STEP Thus the refinement of WEP tests still continues. For a 
detailed summary of the current experimental status and historical overviews see [2*4*] . 

Historically, the next important theoretical development after Newton with respect 
to the equivalence principle was given by Einstein in the context of his theory of general 
relativity in 1915 [B]. While the WEP formulated in the language of special relativity 
demands that in a sufficiently small, free falling laboratory all mechanical laws of physics 
are the same as if gravity is absent, Einstein generalized this statement from mechanical 
to all laws of physics which is formulated as the 



Einstein Equivalence Principle (EEP): 

All physical laws of special relativity are valid in the presence of a gravitational field in 
an infinitesimally small, free falling laboratory. 



Together with the principle of general covariance, the EEP provides the foundation of 
general relativity and hence of the idea that gravity is a phenomenon of curved spacetime. 
For this reason the necessity of having a solid experimental verification as well as a critical 
analysis of the underlying concepts and possible connections between the WEP and the 
EEP becomes very clear. These topics will be developed within the next sections. 
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1.1.1 Schiff's Conjecture 

The clear distinction that was made in the early days between the basic concepts of the 
WEP and the EEP has become increasingly blurry today. Test masses are composed of 
atoms where the constituents, the protons, neutrons and electrons, interact via the mass- 
energy of the electromagnetic, strong and weak interaction. Validity of the WEP in this 
context implies that all nongravitational fields couple in the same universal way to gravity 
so that, finally, measurements of the freefall accelerations turn out to be profound tests of 
the EEP as well as gravitational redshift measurements. This plausibility argument also 
further supports Schiff's Conjecture, originally invented by Leonard I. Schiff in 1960 [H]: 
Any complete, self- consistent theory of gravity that embodies WEP necessarily embodies 
EEP. 

Here a theory is defined as being "complete" when it is capable of making definite 
predictions about the result of any experiment, within the scope of a theory of gravity, 
that the current technology is able to perform. In this sense Milne's kinematic relativity 
[TU] must be considered as an uncomplete theory since it makes no gravitational red- 
shift prediction. A theory is called "self-consistent" if the prediction for the outcome 
of any experiment within its scope is unique and does not depend on the way it was 
derived. According to this definition Kustaanheimo's various vector theories ^l] must 
be seen as being inconsistent since the results for light propagation are different for light 
viewed as waves and for light viewed as particles. Following the argumentation above 
general relativity provides an example of where Schiff's conjecture is validated since it 
describes gravity by a second-rank symmetric tensor to which all matter fields couple 
universally. 

It is obvious that Schiff's conjecture, if valid, would have a strong impact on gravita- 
tional research since, e.g., Eotvos experiments could then be seen as direct tests of the 
EEP and the idea of gravity as a phenomenon of curved spacetime. It is generally recog- 
nized, that a rigorous proof in a mathematical sense of such a conjecture is impossible, 
since such a proof would require an at least moderately deep understanding of all gravi- 
tation theories that satisfy WEP, including those not yet invented, and never destined to 
be invented |Sj. Nonetheless, a number of plausibility arguments have been formulated 
within the past decades. The most general and elegant of these consists of a simple cyclic 
gedanken experiment under the assumption of energy conservation and was first formu- 
lated by Dicke in 1964 [7] and subsequently developed by Nordvedt (1975) [Sj and Haugan 
(1979) [9J. A more qualitative argumentation given by Thorne, Lee and Lightman 
is founded on Lagrangian-based theories of gravity and is similar to the one mentioned 
above. 

1.2 Theoretical context for analyses of EEP tests 

The relativity revolution and the quantum revolution are certainly among the greatest 
successes of 20th century physics. Both have changed our view of space, time and matter 
in a radical way but the underlying concepts of the theories they produced are unfortu- 
nately fundamentally incompatible. For example general relativity is a purely classical 
theory where each particle has simultaneously a definite position and momentum in a 



6 



CHAPTER 1. INTRODUCTION 



given spacetime point, whereas quantum mechanics tells us that this is only approxi- 
mately true for macroscopic objects within a region where the spacetime curvature is 
much greater than the position uncertainty. This conceptual tension becomes even more 
obvious by including the weak equivalence principle. The WEP, statet in an alternative 
formulation, says that if an uncharged test body is placed at an initial event in spacetime 
and given an initial velocity there, then its subsequent trajectory will be independent of 
its internal structure and composition [7]. Here, an uncharged test body is meant to 
describe an electrically neutral test mass with negligible self-gravitational energy that is 
small enough that inhomogenities and therefore tidal effects of the external gravitational 
field can be ignored. So, given two test masses which should be used for testing the WEP, 
the locality of this principle requires that the volume of space between both trajectories 
must go to zero before the statement becomes exact [T^j. That is exactly the point where 
the WEP comes into conflict with the uncertainty principle since this limiting process 
causes an infinite uncertainty in their momenta and, hence, makes any prediction about 
the trajectory impossible. A simple gedanken experiment which reveals a violation of the 
WEP in this sense was given by R. Chiao in [Tl|: Two perfectly elastic balls with differ- 
ent chemical compositions were dropped from the same height above a perfectly elastic 
table. WEP predicts that the vertical trajectories as well as the subsequent oscillations 
are identical and indistinguishable since the total amount of kinetic and potential energy 
remains constant. However the time-energy uncertainty in quantum mechanics indeed 
allows the balls to propagate into the classicaly forbidden regions above their turning 
points. Since tunneling depends on the mass and therefore on the chemical composition 
of the object this effect would represent a quantum violation of the weak equivalence 
principle. This example of a possible quantum violation of the WEP underlines again 
the discomfort many physicists feel with having two fundamental theories, both so far 
experimentally verified with an outstanding precision but without a satisfying common 
interface. It would certainly go beyond the scope of this thesis to summarize, even in 
a brief way, all the approaches physicists and philosophers have taken within the last 
70 years to resolve this conflict. The interested reader may therefore have a look at the 
exellent reviews by, e.g., Rovelli [TH] and Carlip [THj which also provide a huge list of 
references for further reading. 

However, from the viewpoint of a relativist a key role in understanding the unification 
problem is certainly played by the EEP and its incorporated WEP (see PHEEj)- Since 
theories which predict violations of the EEP are numerous and the experimental guidance 
is so far negligible, it is important to establish a systematic theoretical framework for 
analysing various experiments and theoretical concepts of different gravitational theories. 
For this reason the first step must consist of providing careful definitions of general 
concepts and notions that every valid theory of gravity has to obey. This procedure can 
easily be regarded as pedantic and even as superfluous since most of the notions from 
everyday gravitational physics and experience seems to be more than obvious without 
further need of explanations. Indeed the next section will show that even the distinction 
between what is a gravitational and what is a nongravitational phenomenon is highly 
nontrivial. Taking into consideration that every theory in physics for historical reasons 
is build up on notions of everyday experience, one has to be very careful by using these 
concepts in more sophisticated theories. Nevertheless, from this starting point general 
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schemes for analysing gravitation theories will be developed. This schemes encompass 
that class of theories which predict violations of the EEP, which are relevant to this work. 

1.2.1 Basic concepts and notions 

Spacetime and Gravitational Theories: Following the notions and definitions 
given by Thorne, Lee and Lightman all gravitation theories can be regarded as a 
subclass of the more general spacetime theories. A spacetime theory basically possesses 
a mathematical formalism which is constructed from a 4-dimensional manifold and from 
geometric objects defined on that manifold [2Z|. Two different mathematical formalisms 
will be called different representations of the same theory if the predictions they produce 
are identical for every experiment. The geometric objects of a particular representation 
are called its variables. The equations which these variables have to satisfy will be called 
the physical laws of the representations, e.g. in the case of general relativity the physical 
laws are the Einstein field equations. 

Gravitational phenomenon: Certainly, the above general scheme is able to encom- 
pass a rich variety of different theories for various physical phenomena. To restrict ourself 
to gravitation theories one simply has to demand that the physical laws the spacetime 
theory provides, must correctly match with generally recognized laws based entirely on 
gravitational phenomena like Keppler's law. This definition immediately requires a clear 
distinction between what is a gravitational and what is a nongravitational phenomenon. 
Already Thorne, Lee and Lightman ^2] mentioned that there seems to be a variety of 
ways in which such a distinction could be made. They suggested to define gravitational 
phenomena as "those which are either absolute or 'go away' as the amount of mass-energy 
in the laboratory (isolated from external influences) decreases". In other words, they sug- 
gested that gravitational phenomena are either prior geometric effects or generated by 
mass-energy. Concerning the first issue one has to reply that the interpretation of gravity 
as a geometric phenomenon is entirely based on the validity of the EEP and, so, is inap- 
propriate for a general theory of gravitational theories. Concerning the second point it 
is important to note that this definition also includes electromagnetic phenomenon since 
the electromagnetic charge is so far restricted to massive elementary particles. If the 
amount of mass-energy could be totally removed from a shielded laboratory, then also all 
electromagnetic phenomena would vanish. It is therefore more appropriate to define the 
classical gravitational phenomena as those which are generated purely by mass-energy, 
regardless of charges. Taking also quantum mechanical properties of matter like charges 
and spins into consideration could therefore certainly lead to a modification of the above 
definition, important for quantum gravity approaches. 

LOCAL NONGRAVITATIONAL TEST EXPERIMENT: An experiment performed in an arbi- 
trary spacetime point is called local and nongravitational if the following conditions are 
satisfied 

• Performed in the center of a freely-falling laboratory. 

• Inhomogenities of the external field can be ignored. 

• Self-gravitational effects are negligible. 
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In addition the laboratory must be impermeably shielded against external electromag- 
netic and other (real or virtual) particle fields. To make sure that external inhomogenities 
of the gravitational fields are unimportant, one has to perform a sequence of experiments 
with decreasing size, until the experimental results reaches asymtotically a constant value. 
An example of a local nongravitational test experiment is a measurement of the fine struc- 
ture constant a, while a Cavendish experiment is not. 

1.2.2 The EEP revised 

Using the last definition it is possible to give an alternative formulation of the EEP 
which is capable of providing a large variety of new sophisticated tests and also reveals 
the important and far reaching symmetries which are inherent to the EEP. 

Einstein Equivalence Principle: 

1. WEP is valid. 

2. The outcome of any local nongravitational experiment is independent of the velocity 
of the freely-falling reference frame in which it is performed. 

3. The outcome of any local nongravitational experiment is independent of where and 
when in the universe it is performed. 



The second aspect in this definition demands that in two frames moving relative 
to each other, all the nongravitational laws of physics must make the same predictions 
for identical experiments. This is therefore called Local Lorentz Invariance (LLI) 
[2U I2H]- The third point in the above formulation of the EEP requires a homogenity 
of spacetime since the outcome of any local nongravitational test experiment must be 
independent of the spacetime location of the laboratory and is therefore called Local 
Position Invariance (LPI). It is important and interesting to make clear that the 
local position invariance not only refers to the position in space but also to the position 
in time. Validity of LPI forces the fundamental constants of nongravitational physics like 
the fine structure constant a or the weak and strong interaction constants not to change 
throughout the lifetime of the universe. For a detailed review and references on this topic 
see !23JEB|- 

The time variation of fundamental nongravitational contants or the gravity-induced 
birefringence of light in gravitational fields, which is the main topic in this thesis, are only 
two examples of tests of certain aspects of the EEP. Indeed many of the experiments in 
gravitational and nongravitational physics are direct tests of the symmetries defined by 
the principles of equivalence. Nongravitational test experiments respond to their external 
gravitational environment during their free-fall and, so, the presence or absence of local 
Lorentz or local position invariance is entirely determined by the form of the coupling of 
the gravitational field to matter [2E]- Therefore, after defining the group of gravitational 
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theories which incorporate the EEP, a formalism capable of representing the coupling 
between gravitational and matter fields for a whole class of gravitational theories will 
be presented. For this purpose, Lagrangian field theory provides a natural setting for 
general considerations. 

1.2.3 Metric theories of gravity 

The validity of the EEP is a crucial distinctive feature regarding the classification of 
various gravitational theories. If EEP is valid then, according to the second and the 
third point of the definition given in Sect. 1.2. 2, the laws of physics which govern a cer- 
tain experiment must be independent of the velocity of the free falling laboratory (local 
Lorentz invariance) and also of its position in spacetime (local position invariance) which 
demands time-independent physical constants. The only laws which are known to satisfy 
these requirements are those of special relativity. According to the first point, validity 
of WEP, it follows that test bodies within the laboratory are moving unaccelerated on 
locally straight lines which can be regarded as geodesies of a metric g, i.e. geodesies in 
a curved spacetime |2"H 123] . These inferences from the EEP are commonly summarized 
in the 

Metric Postulates 

• Spacetime is endowed with a symmetric metric g. 

• The trajectories of freely falling bodies are geodesies of that metric. 

• In local freely falling reference frames, the nongravitational laws of physics are those 
of special relativity. 



Every theory which embodies the EEP necessarily includes the metric postulates and 
is therefore called a metric theory of gravity. As a consequence, in every metric theory 
all nongravitational fields couple in the same way to a single second rank symmetric 
tensor field which is called universal coupling. This means that the metric itself can be 
viewed as a property of spacetime itself rather than as a field over spacetime. For this 
reason one can say, that EEP serves not only as the foundation of general relativity but 
of the more general idea of gravity as a curved spacetime phenomenon. However, it is 
important to note that it tells us nothing about how spacetime is curved, i.e. how the 
metric is generated. Actually for this reason it is possible that besides the metric other 
gravitational fields such as scalar, vector or tensor fields could exist which only modifies 
the way in which matter and nongravitational fields generate the metric. Nevertheless, 
in order to preserve universal coupling only the metric acts back on matter and non- 
gravitational fields in the way, prescribed by EEP. For example in general relativity the 
metric is generated directly by the stress-energy tensor of matter and nongravitational 
fields, whereas in the Brans-Dicke theory |31j . besides general relativity the most famous 
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representative from the class of metric theories, matter and nongravitational fields first 
generate a scalar field (p. Then acts together with matter and other fields to generate 
the metric but it only couples indirectly to matter so that the theory remains metric. 

From these two examples one can see that the main feature which distinguishes dif- 
ferent metric theories is the number and the kind of the additional gravitational fields 
they contain and, in turn, the equations which govern the evolution and the structure 
of these fields. Whether or not a theory of gravity exhibits the symmetries defined by 
EEP depends therefore entirely on the manner in which the theory couples the metric 
to matter and nongravitational fields. This aspect later becomes very important with 
respect to nonmetric couplings which lead to gravity-induced birefringence of polarized 
light. 

Nonmetric theories of gravity like Moffat's nonsymmetric gravitation theory (NGT) 
[SB] or the metric-affme gauge theory of gravity (MAG) (SEj violate, by definition, one or 
more of the metric postulates and hence it is more necessary than surprising that they 
predict novel couplings between gravitational and nongravitational fields. An appropiate 
framework for general considerations of theses aspects is provided by Lagrangian field 
theory which will be discussed in the next section. 

1.2.4 Lagrangian based theories of gravity 

According to Thorne, Lee and Lightman |12j . a generally covariant representation of a 
spacetime theory is called Lagrangian-based if 

1. There exists an action principle that is extremized only with respect to variations 
of all dynamical variables (For a detailed definition of dynamical variables see p. 8 
of Thorne et al.). 

2. The dynamical laws of the representation follow from the action principle. 

A theory is called Lagrangian-based if it possesses a generally covariant, Lagrangian- 
based representation. Examples are general relativity as well as the Brans-Dicke theory 
and also NGT and MAG. An action principle of the form 



encompasses metric as well as nonmetric theories of gravity. Here, ip m and ip g denotes 
the corresponding (quantummechanical) matter and gravitational fields respectively of a 
given theory. Then the final objective is, in the end, to break down the explicit structure 
of a particular action. Like in conventional Langrangian field theory this strategy is 
supported by the existence of certain symmetries which enforce definite restrictions on 
the form of the action. In the case of gravitational theories the symmetries of LLI and 
LPI are consequences of EEP so that many experiments in gravitational physics are direct 
tests of the structure of a given lagrangian density. The investigation of this aspect can 
be deepend by splitting a given lagrangian density into a purely gravitational part C g 
and a nongravitational part C ng . 





C = C g + C 



ng 
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EEP 
Violated 



Figure 1.1: Starting from the two known approaches which can lead to a relativistic 
theory of gravity, EEP and action principle, different ways could result in a theory which 
violates EEP. Details are given in the text. 



The gravitational density depends entirely on the gravitational potentials and their 
derivatives, so that its structure determines the dynamics of the free gravitational fields 
in the theory. The nongravitational part C ng also depends on the gravitational poten- 
tials and their derivatives but additionally on the matter fields and their derivatives so 
that the form of C ng specifies the coupling between matter and gravity, i.e. how matter 
responds to gravity and how matter acts as a source of gravity [32J. 

Matter equations of motion which predict the outcome of a local, nongravitational 
test experiment are determined only by the form of C ng and, therefore are derived from 
the action principle 

5j£ ng (ij m ,^)d*x = . (1.3) 

In this equation, only the matter fields ip m are variable in a true sense, whereas the 
external gravitational fields in which the experiment is performed are viewed as static. 
In the case of metric theories, the symmetries defined by EEP forces C ng to a metric 
form, i.e. C ng have to couple a single symmetric tensor gravitational field universally to 
all nongravitational fields. For nonmetric theories C ng can admit several EEP violating 
couplings so that detailed investigations of the structure of this lagrangian density with 
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respect to experimental tests requires a broader framework which is able to encompass 
classes of conceivable nonmetric couplings. This will be discussed within the next sections. 

A broad class of experiments in gravitational physics are those where the self-gravita- 
tional effects are not negligible. The underlying symmetries are defined by the Strong 
equivalence principle [7] which, as a generalization of EEP, also refers to local grav- 
itational test experiments and demands that the symmetries are analogous to those of 
EEP. Then, the equations of motion follow from the action principle 

5 J £(Vw,^°+Vi e Vz = (1.4) 

where the matter field ip m , and the local gravitational fields ip® are varied but the 
external field ip^ is assumed to be static. General relativity provides an example of a 
theory which exhibits the symmetries defined by the strong equivalence principle. 

Having laid down the basic classes of relativistic theories of gravity, the remaining part 
of this section is devoted to the question in which class one could expect violations of WEP 
and EEP respectively. Basically only two ways are known in which a set of gravitational 
laws could be combined with the special relativistic, nongravitatational laws of physics. 
Classical starting point of the first approach is the EEP where gravity is described by 
one or more fields of tensor ial character including the metric. By requiring that in 
local Lorentz frames the nongravitational laws of physics take their special relativistic 
form one arrives with the principle of general covariance at a metric theory of gravity 
which, by construction, includes EEP. Starting point of the second approach is a general 
action principle with a relativistic lagrangian where one arrives at a relativistic theory 
of gravity in the way described aboved. Assuming a universal coupling, there are again 
two possibilities: Taking the usual second rank symmetric tensor one arrives at a metric 
theory which obeys EEP. Whereas taking a nonsymmetric metric, Will has shown 
that even if coupled to matter fields in a universal way, those theories violate WEP and 
thus EEP. Rejecting the idea of a universal coupling leads to a nonmetric theory which, 
by definition, violates at least one of the metric postulates. Now, since Schiff 's conjecture 
states that any complete, self-consistent theory of gravity that obeys WEP also includes 
EEP, it therefore suggests that nonmetric, relativistic, Lagrangian based theories should 
always violate WEP. 

1.2.5 General theoretical frameworks 

Although every experiment performed so far is in nearly perfect agreement with gen- 
eral relativity, its conceptual problems (i.e. incompability with quantum mechanics, see 
Sect. 1.2) have led to the development of numerous alternative theories of gravity within 
the last decades. Since the outcome of a certain experiment is, in most cases, not only 
relevant for a special but for a whole class of theories with similar characteristics it has 
become essential to have theoretical frameworks for the classification and comparison of 
various approaches. Within such a framework one can seek for conceptual differences 
and similarities and also compare predictions for a variety of experiments. Among these 
schemes the Dicke framework could be seen as a starting point for enhanced models like 
the Tife/z-formalism and the x^-formalism which have become crucially important for 
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designing and interpreting experiments that have the ability to reveal possible violations 
of the EEP. 

Dicke Framework 

The Dicke framework, given 1964 in Appendix 4 of Dicke's Les Houches lectures [7j 
imposes several very fundamental constraints that every acceptable theory of gravity has 
to obey. Assuming that nature likes things as simple as possible he suggests that the only 
geometrical concepts introduced a priori in a spacetime theory are those of a differentiable 
four-dimensional manifold, with each point in the manifold corresponding to a physical 
event. The manifold need not have either a metric or an affine connection, whereas one 
would hope that experiments will lead to the conclusion that it has both. Furthermore, 
the dynamical equations should be constructed in a generally covariant form to avoid 
arbitrary subjective elements in the equations of motion which merely reflects properties 
of a particular coordinate system. 

After formulating these mathematical constraints, Dicke requires two aspects to be 
fulfilled by every viable theory of gravity: 

1. Gravity must be associated with one or more fields of tensorial character. Since 
nature abhors complicated situations, interactions involving one field will occur 
before involving two or more, which favours the viewpoint that gravity is associated 
with only one field. 

2. Having in mind the close connection between variational principles and conserva- 
tion laws, all dynamical equations that govern gravity must be derivable from an 
invariant variational principle. 

The Dicke framework can be seen as the prototype of all successive classification 
schemes. The assumptions and constraints placed on all acceptable gravitational theories 
are often used as a basic building block in the development of enhanced schemes like the 
THefj, and the xg-formalism as will be shown below. The main achievement of the Dicke 
framework is that it supports the design and interpretation of experiments which address 
fundamental questions on the nature of gravity, like what types of fields (scalars, vectors, 
tensors of various rank) are associated with gravity. For specific questions especially about 
the motion of electromagnetically charged particles in gravitational fields the THe/i- 
formalism was developed. 

Tife/i-formalism 

The Tife/i-formalism, developed in 1973 by Lightman and Lee |Sj encompasses all met- 
ric and many nonmetric theories of gravity. It describes motions and electromagnetic 
interactions of charged particles in an external static spherically symmetric (SSS) grav- 
itational field which is described by a potential U. The equations of motion of charged 
particles in the external potential are characterized by two arbitrary functions T(U) and 
H(U) while the response of the electromagnetic field to the external potential, the gravi- 
tationally modified Maxwell equations, are characterized by the functions e(U) and fi{U). 
Clearly, the explicit forms of the phenomenological gravitational potentials T, H, e and 
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/i varies from theory to theory. One of the most important results of this formalism is 
that certain combinations of these functions reflect different aspects of EEP j2Hj which 
could be shown by means of the corresponding action. 

Within the Tife/i-formalism the nongravitational laws of physics can be derived from 
an action I^q for a structureless test particle a with restmass m 0a , charge e a and elec- 
tromagnetic fields coupled to gravity, given by the sum 

J NG = J + knt + hm (1.5) 

where the motion of a free (neutral particle) with coordinate velocity v% = dx^/dt on the 
particle world line x%(t) is described by 

/o = -E m o« f{T-Hvl) l ' 2 dt . (1.6) 

The interaction of the particle a with the electromagnetic field follows from 

/mt = E e - / A Xdt (1-7) 

where represents the electromagnetic vector potential F^ v = A v ^ — A^ u . Finally, the 
coupling of the electromagnetic to the gravitational field is given by 

iem= {^Y 1 J{eE 2 - l i- x B 2 )d i x . (1.8) 

with E = VA - dA/dt and B = V x A. 

Basically, I^q violates Local Lorentz Invariance [21] • A metric theory is obtained if 
and only if e and \x satisfies 

e = /i = (H/T) 1/2 or T' 1 He~ 1 /j l '~ 1 = 1 (1.9) 

for all U. The quantity (T^f/e -1 /! -1 ) 1 / 2 plays the role of the ratio of the speed of light 
c light t° the limiting speed of neutral test particles Co Analogously, one gets that I^q is 
Local Position Invariant if and only if 

eT i/2 H -\/2 = constant ( L10 ) 

fiT 1/2 H- 1/2 = constant (1.11) 

independent of the position P(x M ) 24J. Nonconstant combinations like above therefore 
directly lead to preferred-positions effects like gravitational birefringence. 

Xg-formalism 

Invented by W.-T. Ni in 1977 |33j . the xg-formalism provides similarly to the THefi- 
formalism a framework for the analysis of electrodynamics in a background gravitational 
field. However, one of the main differences is that the xg-formalism is not restricted to 
static, spherically symmetric gravitational fields. The \ of its name refers to a tensor 
density which provides a phenomenological representation of the gravitational fields. The 
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structure of the xg-formalism is in agreement with the basic assumptions and constraints 
of the Dicke framework. Furthermore it is assumed that in the absence of gravity, the 
nongravitational Lagrangian density C^q reduces to the special relativistic form. Re- 
specting these assumptions together with demanding electromagnetic gauge invariance 
and linearity of the electromagnetic field equations the most general Lagrangian density 
becomes 

£ N G = -^X aM5 F a pF lS . (1.12) 

Within this general formalism, the independent components of the tensor density x al3lS 
comprise 21 phenomeno logical gravitational potentials which allows one to represent grav- 
itational fields in a very broad class of nonmetric theories. In the case of metric theories, 
the tensor density is constructed alone from the metric tensor 

= X ^Z-g (pong* _ g/CtgPf} . (L13) 



2 

The Lagrangian density has then the usual metric form like the one from general rel- 
ativity and, therefore, incorporates the coupling between the electromagnetic field and 
a metric gravitational field. As already mentioned, nonmetric theories involve various 
other gravitational potentials in addition or instead of the metric one. In this case the 
Lagrangian density of the xg-formalism describes the coupling between the electromag- 
netic field and the various, mostly nonmetric, gravitational potentials which could lead 
to EEP violating effects. 



1.3 Nonsymmetric Gravitation Theory 

The idea of describing the gravitational field by means of a nonsymmetric second rank 
tensor field traces back to the work of Einstein and Strauss in 1946 who tried to formulate 
a unified theory of gravitation and electromagnetism [SHUH!!- By using the decomposition 

9^ = 90^+9^] (1-14) 
where ^ ^ 

g^v) = -{g^v + g„»), S[h = -(g^ - g vv ) (1.15) 

their intention was to interpret the symmetric guv) a s an expression of the gravitational 
field and the antisymmetric part gu^ as an expression of the electromagnetic field. Unfor- 
tunately, it was soon realized that g^ v ] could not describe physically the electromagnetic 
field without serious contradictions, so that the nonsymmetric ansatz vanished for more 
than 30 years. 

In 1979 Moffat picked up this issue again and published the first version of his non- 
symmetric gravitation theory (NGT) (SH| which is based on a non-Riemannian geometry 
according to ()1.14j) and the analog expression for the affine connection 

= ^(ixu) + rf Mi ,j • (i.i6) 

The motivation of NGT is to construct the most general classical description of space- 
time that contains general relativity as a special, low-energy case. The main conceptual 
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difference to Einstein- Strauss theory is that the nonsymmetric field structure now rather 
describes a generalization of Einstein gravity than a unified field of gravitation and elec- 
tromagnetism. We consider here the physical implications which follows from the NGT 
version, described by Moffat in 1990 [13], and set sharp limits on the essential coupling 
constant £ 2 , responsible for violations of EEP. Although it became clear in the meantime 
that this theory as well as a later published modification [33] suffers from serious problems 
like ghost poles, tachyons and higher order poles (see [3H] for a list of references) NGT 
nevertheless serves as a prototype for a whole class of nonmetric gravitational theories 
which predict spatial anisotropy and birefringence. Setting sharp and reliable limits on 
C 2 is therefore not only a further test of NGT but rather adresses the question on the 
physical relevance of gravity-induced birefringence in principle. 
Defining the contravariant tensor g^ v in terms of the equation 



the Lagrangian density with matter sources is given by 

^NGT — £r + £m 

with 

C R = gTRAW) - 2Av^ 

and 

c 4 3 

Here, A is the cosmo logical contant and g Miy = yf—gg^, while the other constants 
their usual meaning. R^iW) denotes the NGT contracted curvature tensor 



2 

defined in terms of the unconstrained nonsymmetric connection 

2 

W x = T A 6 x Wu 



where 



2 

The full NGT field equations with matter source therefore become 



G„ v (W) = — r T llv + Kg lxv , 
K [H =4ttS 



where 

G^(W) = R^-~g^R 



1.17) 

1.18) 
1.19) 

1.20) 
have 

1.21) 

1.22) 
1.23) 

1.24) 
1.25) 

1.26) 



So, in addition to a conserved nonsymmetric energy-momentum tensor T^ u (and in con- 
trast to general relativity), NGT also contains a conserved- vector-current density S 1 ^, 
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where the current conservation arises by Noether's theorem from the invariance of the 
Lagrangian density ()1.18jl under the transformations of an Abelian U(l) group, i.e. 

g^U* = 47rSM ,M = ■ (1.27) 
It suggests itself to interpret as the conserved particle number of the fluid 

s" = £/iWf ■ ( L2 8) 

i 

Here, ff is a coupling constant for each species % of fermions, rij is the constant fermion 
particel number and w M = dx^/dr denotes the proper-time velocity of the particle. The 
so-called NGT charge £ 2 is defined as 

i 2 = J S°d 3 x , (1.29) 

having the dimension of [length] 2 . Since i 2 arises from a conserved current, it has been 
postulated that it is proportional to conserved particle number |43| . In order to under- 
stand the exceptional position of fermions in NGT one has to recall that in the NGT 
scheme, the nonsymmetric tensor leads to a nontrivial extension of the homogeneous 
Lorentz group £0(3, 1) of general relativity to the local gauge group GL(A,R). This 
extension describes the most general transformations of the linear frames that contains 
the homogeneous Lorentz group as a subgroup [10]. Furthermore, the group GL(4,R) 
contains only infinite-dimensional spinor representations which leads to the idea that 
particles are described as extended objects in contrast to general relativity where point 
particle fermions are conventionally described by finite, nonunitary representations of the 
homogeneous Lorentz group S0(3, 1) [101 E]. Therefore, fermions play an important role 
in NGT and are given special consideration as a source of the (antisymmetric part of the) 
gravitational field |42j . 

Despite various arguments, like beauty and simplicity of the theory, one could have 
in favour of NGT it is merely a, perhaps elegant, hypothesis because of the complete lack 
of any experimental verification. One therefore has to look for the physical implications 
of this idea which implies constraining the essential coupling constant I 2 since history is 
full of elegant hypothesis, later contradicted by nature. 

1.4 Metric- Affine Gravity 

Questioning critically the foundations of a certain physical theory is often the first, 
promising step for getting a deeper insight into the basic mechanisms ruling the cor- 
responding phenomena. In this sense it is certainly important to realize that the weak 
equivalence principle as it is formulated by Newton only gives a prediction for the be- 
haviour of macroscopic objects and, in the language of general relativity, their interaction 
and influence on the structure of spacetime. Microscopic properties of matter like the 
spin angular momentum of particles are totally neglegted as they average out on the 
macroscopic scale which justifies the accusation that general relativity is blind to the 
microscopic structure of matter. So, the hypothesis is near at hand that spin angular 
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momentum might be the source of a gravitational field too, since it certainly characterizes 
matter dynamically in the microphysical realm. 

Although this issue seems to be evident, the question immediately arises what kind 
of "interface" or ansatz one has to choose in order to extend the concepts of general 
relativity to emcompass quantum mechanically relevant observables. An answer to this 
is very likely given by looking at gravity in a gauge theoretical way which is favoured by 
the successfull description of the other three fundamental interactions by means of gauge 
theories of underlying local symmetry groups. In this sense Utiyama [IB] has shown in 
1956 that general relativity could be recovered by gauging the Lorentz SO(l,3) group. 
Nevertheless this procedure is unsatisfactory since the conserved current associated to 
the Lorentz group via the Noether theorem is merely the angular momentum current 
which alone cannot represent the source of gravity. This problem was solved by Sciama 
and Kibble in 1961 JTTt EHl EH] who prooved that it is really the Poincare group as the 
semi-direct product of the translation and the Lorentz groups, which underlies gravity. 
This scheme now allows spin angular momentum to be included. In analogy to the 
coupling of energy momentum to the metric, spin is coupled to a geometrical quantity 
which is related to rotational degrees of freedom in spacetime. This concept leads to a 
generalization of the Riemannian spacetime of general relativity to the Riemann-Cartan 
spacetime U4. In a C/4 space the affine connection F^a is not symmetric in the lower 
indices a and (3 which leads to a nonzero torsion tensor 

ZV= 5(1^-1%) , (1-30) 

initially introduced by E. Cartan [20], as the antisymmetric part of the affine connection. 
Analogously to Riemannian spacetime it is required that the local Minkowski structure is 
preserved, i.e. that the line element is invariant under parallel transfer. The deformation 
of length and angle standards during parallel transport is measured by the so-called 
nonmetricity one-form, defined by 

Q a p := Dg aP . (1.31) 

In Riemann-Cartan spacetime U4 as well as in Riemannian spacetime V4 of general rela- 
tivity and in the Minkowski spacetime R4 of special relativity the nonmetricity vanishes 

Q a p = . (1.32) 

The geometrical structure of U4 is that of an n-dimensional differentiable manifold M n 
where at each point of M n , there is an n-dimensional tangent vector space Tp(M n ). The 
local vector basis e a can be expanded in terms of the local coordinate basis <9, := d/dx l 

e a = e\di (1.33) 

where a, (3 = 0, 1, 2, . . . , (n—1) are anholonomic or frame indices and i, j,k, . . . , (n—1) are 
holonomic or coordinate indices. In the cotangent space Tp(M n ) there exists a one-form 
or a coframe 

^ = e/dx i . (1.34) 
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Defining further the one- form r\ a ^ = *($ Q A $p A with the Hodge star * and the 
three-form r; a = * q 9 a the field equations of Einstein-Cartan theory read 

T a (1-35) 
i 2 r a( 3 (1.36) 

where A denotes the cosmo logical constant, £ 2 Einstein's gravitational constant and R^ 1 
the curvature- 2-form. J2a denotes the canonical energy-momentum current of matter. 
While the first equation relates curvature to energy momentum, the second equation 
provides a link between the spin angular momentum tensor r a p to Cartan's Torsion. 
It is now obvious how Einstein-Cartan gravity, general relativity and special relativity 
are connected. Starting from a Riemann-Cartan spacetime U4 the usual Riemannian 
spacetime of general relativity is recovered by neglegting torsion. If additionally curvature 
vanishes one gets the Minkowski spacetime R 4 of special relativity. 

U 4 ^ V 4 ^ R 4 (1.37) 

The Einstein-Cartan theory of gravity is a viable theory since it is in agreement 
with all experiments performed so far [HH E21 EH1 • Under usual conditions the spin r Q( g 
averages out and can be neglegted which in turn, according to the second field equation, 
implies vanishing torsion and, so, general relativity is recovered. The additional spin-spin 
contact interaction shows up only at extremely high matter densities (~ 10 54 g/cm 3 ) and 
therefore has not been seen so far since even typical neutron stars have only densities of 
the order of (10 15 g/cm 3 ). 

The reason why the Einstein-Cartan theory has been considered here is, that it pro- 
vides the simplest model of the so-called metric-afline gauge theory of gravity (MAG) [54, 
153} 1ST)] , representing the most general canonical gauge theory. Metric-afline theories are 
build upon the more general affine group A(n, R) which is the semidirect product of the 
translation group and the group of linear transformations, i.e. A(n, R) = T n <zGL(n, R). 
This transformation group acts on an affine n- vector according to 

£_>£' = Ae + r (1.38) 

where A = {A a p} G GL(n,R) and r = {r Q } G R a . The transformations fHffijl of the 
affine group represent a generalization of the Poincare group where the pseudo-orthogonal 
group SO(l,n — 1) is replaced by the general linear group GL(n,R). The reason for 
introducing A(n, R) is the assumption that physical systems are indeed invariant under 
the action of the entire affine group and not only invariant under its Poincare subgroup. 
The physical symmetries which are added by the general affine invariance are the dilation 
and shear invariance. Both are of physical importance since dilation invariance is a crucial 
component of particle physics in the high energy regime and shear invariance was shown 
to yield representations of hadronic matter. Further, the corresponding shear current 
can be related to hadronic quadrupole excitations [5BJ EH] Therefore, the additional 
appearance of symmetries in the high energy regime could be taken as an indication that 
metric-afline gravity has played an important role at the early stages of the universe and 
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reduces to general relativity and translational invariance in the low-energy limit after 
some symmetry breaking mechanism |56{ 157] . 

Metric affine gravity uses the metric g a p, the coframe and the linear connection 
Tj 3 to represent independent gravitational potentials. This is summarized in Table 1. 



Potential 


Field strength 


metric g a /3 
coframe d a 
connection F a P 


nonmetricity Q a p = —Dg a p 

torsion T a = Dd a 

curvature Rj 3 = dTj 3 - IV A Yj 3 



Tab. 1: Gauge fields in metric-affine gravity. 



In contrast to NGT, metric-affine theories do not allow for an antisymmetric part in 
the metric tensor, since it does not lend itself to a direct geometrical interpretation [5^] . 
Although the symmetric tensor is referred to as the metric, metric-affine gravity becomes 
nonmetric when the new gravitational potentials or their derivatives (the nonmetricity, 
torsion and curvature gravitational fields) couple directly to matter, as they generally do. 
Nonmetric couplings to the electromagnetic field are what can lead to gravity-induced 
birefringence which will be discussed in the next section. 

1.5 Electrodynamics in a background gravitational 
field 

1.5.1 Birefringence in nonsymmetric theories 

Moffat s nonsymmetric gravitation theory (NGT) can be regarded as the prototype of 
a diverse class of Lagrangian-based nonmetric theories where a nonsymmetric tensor 
field does not couple universally to matter stress energy. Especially the coupling of the 
antisymmetric part of the nonsymmetric gravitational field to the electromagnetic field 
leads to a polarization dependent delay and, so, to an alteration of a polarization signal 
which is a consequence of the breakdown of EEP in these theories |21[l22j. The following 
analysis was first published by Gabriel et al. 51. from whom most of the notations are 
adopted. 

The electromagnetic field equations which govern the propagation of light through a 
nonsymmetric gravitational field can be derived from an action principle. A general form 
for this action, quadratic in both the electromagnetic field strength F^ u = and the 
inverse metric, was given by Mann et al. |2H] 
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Iem = -y^ J d^x^gTg^g^ (ZF,„F aP + (1 - Z)F au F^ + YF m F vfi ) 



1.39) 



The matrix g^ v denotes the inverse of the nonsymmetric gravitational field g^ defined 



by g^ a g v 



Y and Z are constants while T is a scalar function which 



cannot depend on the electromagnetic field and which must be unity in the Einstein- 
Maxwell limit g^v] — > 0. This implies that T = J r ( v /— g/y/— 7), where g = detg^ and 
7 = det^). 

Within this scheme a static, spherically symmetric gravitational field like that of the 
Sun is described by an isotropic coordinate system, centered on the sun. The symmetric 
part of the field takes the form g 00 = — T(r), g/ 0i \ = 0, and guj\ = H{r)5ij where T and 
H are functions of r = \x\. The theories which are encompassed by this formalism like 
NGT provide a representation of the antisymmetric part of the gravitational field that 
can be expressed in an isotropic coordinate system as g[ q = L{r)n,i and g^ = 0, where 
rii = Xi/r. At this point it should be mentioned that the polarization dependent delay 
that this analysis reveals reflects this special form of the nonsymmetric gravitational 
field. Gabriel et al. point out in their paper that a similar analysis, based on a more 
general representation also reveals polarization dependence. 

Employing this special representation together with definitions for the electric and 
magnetic fields via Ej = Ej and Fjk = ejkiBi this yields 



57T 



with 



d x 



eE 2 + Xea(n- E) 2 - -B 1 



// 



(n-Bf 



1.40) 



and 



a = 



n 



H- 


1/2 




L 2 ' 


7f_ 




1 - 


TH 



H- 


1/2 


[i- L2 l 






TH 



T 

2L 2 
TH 
L 2 
TH 



X= 1-Y-Z 



1/2 



TH 



1.41) 

1.42) 

1.43) 
1.44) 

1.45) 



In order to derive the electromagnetic field equations from the general action (jl.40j) 
the special form of the coupling between the nonsymmetric gravitational field and the 
electromagnetic field in terms of X, Y, Z and T must be given. Gabriel et al. used the 
formF = 1 — ZandjF= (\—L 2 /TH) 1 ^ 2 = y/—g / y/— 7 by demanding e = [i in accordance 
with NGT. However, it is important to note that similar analyses of theories having other 
values of Y and Z also reveals a polarization dependence in delay measurements. 
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Using this special coupling, the action (jl.40|) reduces to 

/em = ^/ d 4 x [tE 2 - -^{B 2 - fi(n • B) 2 )^ (1.46) 

This action differs from the usual Einstein-Maxwell action mainly because of the presence 
of the fl(n • B) 2 term. Will [26] has pointed out, that such a term may produce pertuba- 
tions in the energy levels of an atomic system, depending on the relative orientation of 
the system's wave function to the direction n. Such pertubations could be constrained 
by ultraprecise energy-isotropy experiments of the Hughes-Drever type jTTl [TH], using 
trapped atoms and magnetic resonance techniques JH] which is presently being inves- 
tigated. Given this action J, the field equations that govern the propagation of light 
through the nonsymmetric gravitational field are 

VxE + — = (1.47) 
at 

V-B = (1.48) 

V • (eE) = (1.49) 



V x 



9(eE) „ 



Qn(n ■ B) 



(1.50) 



The first two pairs follow from the usual definitions of E and B in terms of electromagnetic 
vector potentials while only the last two pairs follow from the action ([1.46)1 . 

To investigate the propagation of electromagnetic waves with respect to polarization 
dependend delay, the appropiate representations of a locally plane wave are 

E = A E e**, B = A B e iq> (1.51) 

Denoting as the gradient of the phase function 

d M $ = (d$/dt, V$) = (-w, k) (1.52) 

the eikonal equation which governs the propagation of a locally plane wave in the limits 
of geometric optics can be derived by inserting the representations (|1.51|) into the above 
set of Maxwell equations, ignoring all derivatives other than those of the phase function. 
Therefore, one gets 

k 2 A B - efiu 2 A B - On • A B [A; 2 n - k(n ■ k)] = (1.53) 

under the assumptions that the wavelength A is much smaller than the typical scale 
on which any of the fields e, /i, Q, k, Ae and Ab vary significantly. This equation is 
similar to the usual dispersion relation, exept the term which is proportional to Q. Since 
the speed of an electromagnetic wave ()1.51|) is given by u/k, it is apparent from the 
structure of (jl.5r>)) that this speed depends on the orientations of k and A# relative to n. 
Therefore, the velocity of an electromagnetic wave propagating through a nonsymmetric 
gravitational field depends on its orientation which, in turn, implies a violation of EEP! 
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Figure 1.2: Visualization of gravitational birefringence in the case of Moffat's nonsym- 
metric gravitation theory (NGT). 



The speed of a linearly polarized wave with its magnetic field lying perpendicular to 
n follows directly from the eikonal equation with n ■ Ab so that one gets 

d = uj/k = (e/i)- 1/2 . (1.54) 

Light propagating in any other direction travels at c± only if Ab is perpendicular to the 
plane spanned by n and k. The same wave, this time with the magnetic field vector 
lying in the plane spanned by n and k, propagates according to k 2 {l — VLsm 1 6) = e/icu, 
where 6 denotes the angle between n and k. The coordinate speed of a wave, having this 
polarization is 

[1 -ft sin 2 #1 1/2 , . 

c e = . 1.55 

e/i 

Since both waves, parallel and perpendicular to the plane, are independent solutions of 
the modified Maxwell equations, they can propagate independently through spacetime 
without changing their polarization. However, because they travel at different veloci- 
ties the "medium" has now two refractive indices, one for each such eigenmode so that 
spacetime has become (linearly) birefringent. 

Now, since light having any other polarization can be viewed as a coherent superposi- 
tion of the two eigenmodes where the difference in the propagation velocities between the 
components changes their initial relative phase and, thus, the light's polarization. The 
accumulated phase shift is calculated by integrating the eikonal equation along a straight 
line, extended far past the Sun (or other astronomical body). The result for light with 
frequency u, given by Gabriel et al. [21], is 

A$ = -co / nsin 2 6(t)dt . (1.56) 
2 Jt 
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The integration of (|1.56|) requires a ray parametrization x(t) = b + k t where the unit 
vector ko denote the ray direction and b is the impact vector that connects the center 
of the sun with the closest point on the ray. When b is smaller than the radius R of the 
Sun (or the star), the portion of the ray inside the object is, of course, of no interest. The 
integration of ()1.56|) is performed from the Sun's surface with t = (R 2 ~ b 2 ) 1 ^ 2 along a 
straight line up to an observer in an infinite distance t\ — oo. 

Restricting the investigations of birefringence to Moffat's NGT with Q = i^/r 4 the 
integration yields 

^ = X& { 16(1 - „ 2 )3/2 - I - ^13^) - 8(1 _ ^2)3/2 — > 

where // denotes the cosine of the heliocentric angle # between the ray's source and the 
Sun's center. Although in a different context, the basic calculations which lead to (jl.57j) 
are given in Appendix B. 

The accumulated phase shift is inversely proportional to the observed wavelength A 
and the radius R. Therefore with respect to experimental tests of gravitational bire- 
fringence, only those objects can be utilized which emit polarized radiation at short 
wavelengths from a sufficiently small radius for a given mass. For a source located at the 
center of the solar disc (/x = 1), the phase shift vanishes, whereas for light emitted from 
the solar limb (// = 0) we have A$ = S^ijjWXR 3 . 

1.5.2 Birefringence in metric-affine gravity 

Soon after Gabriel et al. j2l] has shown that nonsymmetric gravitation theories like 
Moffat's NGT predict gravitational birefringence and, so, a violation of the Einstein 
equivalence principle, Haugan and Kauffmann (2Hj remarked that this phenomenon could 
be extended to the far more diverse class of nonmetric theories encompassed by the \g- 
formalism. Although already discovered in 1984 by Ni this result was overlooked up 
to that time because no gravitation theories predicting such birefringence were known 
and the available techniques for testing these predictions were not sufficient. 

Starting with the nongravitational Lagrangian density (jl.12)) of the xgAoxmaXvspca. 
in the limit of weak gravitational fields, Haugan and Kauffmann [20J gave a general 
prediction for the accumulated phase shift A</> by using the methodology of geometric 
optics. The detailed calculations are given in Appendix A. Basically, the relative phase 
shift is simply a function of the difference in velocity between the two eigenmodes denoted 
by 5c/ c plus a small second order correction from the tensor density x al3lS ■ 

A$ = oj J dt + 0(5 X 2 ) (1.58) 

The explicit form of ()1.58|) depends of course on the phenomenological representation of 
the gravitational potential and their couplings which is provided by x al3 ^ s . In the case of 
metric-affine gravity one first has to answer the question which of the gravitational fields 
one has to couple to electromagnetism. Nonmetricity is a rather exotic possibility since 
it is assumed that it only plays a relevant role in the very high energy regime like the 
early stages of the Big Bang. It is therefore more suggestive to think of torsion which 
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couples to electromagnetic fields. However, the form that such coupling could have is 
not exactly clear and only very little has been done is this direction so far. Indeed there 
are numerous nontrivial ways in which torsion could couple to the electromagnetic field. 
However, we decided to use 



since this form is equivalent to a particular fourth-rank 5x with tensorial character and 
such a term could, as we have shown in [84J, induce birefringence. In analogy to the 
nonsymmetric charge £ 2 , the strength of the coupling is described by a possibly material 
dependent constant k 2 (see Appendix A) having the dimension of length. Our intention 
is to set strong limits on k 2 and, so, to decide about the physical relevance of gravity- 
induced birefringence. Since different astrophysical objects (Sun, white dwarfs, active 
galatic nuclei) may have different chemical compositions, it is important to set and to 
compare limits on k 2 for a variety of different objects which is one of the objective targets 
of this thesis. 

Since we are going to look for birefringence in the spherically symmetric gravitational 
fields of stars, we are interested in static and spherically symmetric solutions of the metric- 
affine field equations with respect to torsion. Such a solution was given by Tresguerres in 
1995 jSZHHH] (see ifCTSjl in Appendix A) which can be split into a nonmetricity dependent 
and a nonmetricity independent part which is assumed to couple to the electromagnetic 
field via (jl.59j) . Using the method of Haugan and Kauffmann [20] the phase shift 1)1.58)1 
takes after tantalizing computations, given in detail in Appendix A, the form 



where u denotes the light's circular frequency and M* the mass of the star in geometrized 
units. Using the same integration technique as outlined in the case of the nonsymmetric 
theories this leads to the explicit form 



The concept of torsion has become an important tool in many present time gravitation 
theories. Quite recently it has been suggested to identify torsion with the field strength of 
a second rank symmetric Kalb-Ramond tensor field which also appears in the low-energy, 
effective field theory limit of string theory jSnHSOj- The rank and symmetry of this field are 
similar to that of the torsion field in metric-affine theories so that it is not unreasonable to 
expect analogous couplings to the electromagnetic field and, consequently, birefringence. 
Explicit calculations on this subject has not been done so far. 



5£ EM = k 2 *(T a AFY(T a AF) 



(1.59) 




(1.60) 





1.6 Description of Polarized Radiation 



If the concept of gravity-induced birefringence is in fact realized by nature, the required 
ingredients for a chance to observe it are certainly given by a source of strong gravitational 
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fields, emitting a reasonable amount of polarized electromagnetic radiation. Since every 
polarized wave can be decomposed into two orthogonal modes, a speed difference between 
them due to birefringence would lead to a phase shift and, hence, to an alteration of the 
initial polarization state. A measurement of this effect, or at least establishing upper 
limits on it, therefore requires some basic knowledge of astrophysical spectropolarimetry. 
A brief introduction to this subject is given here so that the results given in the following 
chapters are understandable. In this section I will not go into details concerning the 
diverse generation processes of polarized radiation since it is from a pedagogical point of 
view by far more useful to do this later in the individual chapters when the corresponding 
knowledge is needed. Instead I will briefly discuss the three main types of polarization 
which are, in practice, best described by means of Stokes parameters. 

For this purpose, we consider the time harmonic representation of a plane, electro- 
magnetic wave, i.e. when each Cartesian component of E is of the form 

a cos(r + 6) = 9£(a e - l(rU) ) (1.62) 

with t — tot — kx. Since the electric and magnetic components are related via 

Bi = v^— ^ , (1.63) 

it is sufficient to consider in the following only the components of E. Hence, a transversal 
wave propagating in the z-direction can be written as 

E x = cii cos(r + 5i) 

E y = a 2 cos(r + 5 2 ) (1-64) 
E z = . 

On the basis of this equations one can distinguish between three types of polarization 
states by considering the nature of the curve which is described by the end points of 

(EMD. ' 

Elliptic Polarization: The most general form of elliptic polarization is recovered by 
squaring and adding the components of (jl.64|) which yields 

(E,\ 2 + (Ey\ 2 _ 2 E.Ey cos6 = s . n 2 S (L65) 

with 6 = 62 — 61. Since the associated determinant is not negative, this equation describes 
an ellipse where, in general, the axes of the ellipse do not coincide with the x and y 
direction. Instead, the ellipse is characterized by (a) the tilt angle ip between the x-axis 
and the major axis and (b) by the ratio between major and minor axis described by the 
quantity tan/5 = £i/£ 2 - 

Circular Polarization: The case of circular polarization is recovered if the polarization 
ellipse degenerates into a circle. A necessary condition for this is that a x = a 2 = a. In 
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Figure 1.3: Polarization ellipse. The geometry and, hence, the polarization state of a light 
beam is completely characterized by means of £i and £2 and the tilt angle ip. Positive 
helicity is defined for an electric vector rotating counterclockwise. 

addition always one component of E has to be zero while the other has its maximum 
which is fulfilled by 



so that (jl.65j) reduces to the equation of a circle E% + E 2 = a 2 . The wave is said to have 
positive helicity if sin S > and negative helicity for sin 5 < 0. 

Linear Polarization: A linearly polarized wave is obtained if the ellipse (jl.65J) reduces 
to a straight line. This is the case for 



such that E x /E y = (—l) m a 2 /a 1 . 

1.7 Stokes Parameters 

The geometry of the polarization ellipse can be completely characterized by means of 
three independent parameters, e.g. the amplitudes a 1; a 2 and the phase difference 5, or 
the major and minor axes £1 and £2 and the orientation angle (p. For the practical use 
in astronomy it is convenient to use a set of four real valued parameters, the so-called 
Stokes parameters I, Q, U and V, first invented by Sir G.G. Stokes in 1852. In terms of 
the amplitudes and phases in (jl.64j) they can be defined as 



S = 5 2 - 61 = mvr/2 (m = ±1, ±3, ±5 . . .) 



(1.66) 



5 = 5 2 -5 l = mit (m = 0, ±1, ±2, ±3. . .) 



(1.67) 



a\ + a\ — Iq + I 90 — J45 + / 135 — Iai rc (+) + Idrc( — ) 



(1.68) 
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Figure 1.4: Pictorial representation of the Stokes parameters. The observer is taken to 
face the radiation source (adapted from Landi Degl'Innocenti [HI])- 

Q = a 2_ G 2 =/o _J 90 ( L69 ) 

U = 2a 1 a 2 cos(<5 2 - 6{) = J 45 - A35 (1-70) 
V = 2aia 2 sin(5 2 - 5i) I circ (-) . (1.71) 

For completely polarized radiation only three of them are independent so that we have 

I 2 = Q 2 + U 2 + V 2 . (1.72) 

The Stokes parameters have the great advantage that they are quadratic in the ampli- 
tudes and, hence, easily obtained from a telescope which is equipped with a polarizer. 
A very evident, operational definition can be given by defining a reference direction in a 
plane perpendicular to the light beam of interest. Setting the transmission axis of an ideal 
polarizer along this reference direction, a measurement at the exit of this polarizer yields 
the value Iq. This procedure is repeated three times after rotating the polarizer clockwise 
by the angles 45°, 90° and 135°, respectively, obtaining the values I45, / 90 and J 135 . The 
linear polarizer is then replaced by an ideal filter for positive circular polarization which 
gives I C irc(+) at exit and, afterwards, by an ideal filter for negative circular polarization, 
measuring Id rc (—)- Then, the operational definition of the Stokes parameters is given by 
(jl.68j) - (jl.71)l which is pictorally summarized in Fig ll.41 Following this definition the 
fractional degree (or percentage) of linear polarization is given by (Q 2 + U 2 ) 1 / 2 /I while 
the fractional degree of circular polarization is simply |V|/7. 



Chapter 2 

Solar Observations 



Since Eddingtons famous measurement of the bending of light at the solar limb during a 
total eclipse in 1919 [H2], the sun is counted among the most important objects concerning 
experimental tests of theories of gravity. Due to it's relative spatial proximity it is possible 
to determine the relevant initial parameters for a particular experiment like the mass and 
the distance with high accuracy. This, together with the high gravitational potential of 
the sun, opens the possibility to look for often tiny effects which makes the difference in 
predictions between competing theories. 

In this sense we utilize solar polarimetric data to constrain birefringence induced 
by the gravitational field of the Sun and set limits on the coupling constants £ 2 and 
k 2 required by NGT and metric-affine gravity. The initial parameter mentioned above 
which is relevant for the major part of this present work represents a prediction about 
the fractional percentage of Stokes profiles with anomalous symmetry properties in solar 
polarimetric data. Based on intensive numerical simulations of the creation of Stokes 
profiles in the solar atmosphere as well as on observations, it is found that the fraction of 
such profiles is always less than 10% of all observed profiles (if this number is sufficiently 
large), independent of the spatial resolution of the instrument. Since a sufficiently strong 
gravity-induced birefringence could produce any desired amount of anomalous profiles 
up to 100%, this value serves as an upper limit for our subsequent analysis. To test this 
approach, we have developed a second technique which is independent of the previous 
assumptions and gives comparable constrains. 

The following chapter starts with a description of the two independent statistical 
approaches, the so-called 'Stokes asymmetry technique' and the 'profile difference tech- 
nique'. The employed data sets, recorded 1995 at the Izana Observatory in Tenerife and 
in 2000 in Locarno, are described in section 2.2 together with the relevant technical details 
of the instruments. We have measured the line profiles of the full Stokes vector in four 
spectral lines which gave us a total of 1120 spectra. The Stokes asymmetry technique 
yields 1% < (57.1km) 2 in line Cr I 5247.56A in the case of NGT and k% < (1.1km) 2 in 
the same line for metric-affine gravity, respectively. These results are consistent with the 
profile difference technique. In the case of NGT, the result is seven orders of magnitude 
smaller than the l\ value favored by Moffat (Moffat 1979 |3E|). 
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2.1 Technique 

We follow two strategies to test for gravitational birefringence. One of these was outlined 
by Solanki & Haugan (1996) [63 . It is summarized and its implementation is described 
in Sec. 2.1.1. The other technique is new and is described in Sec. 2.1.2. 



2.1.1 Stokes asymmetry technique 

The strategy proposed by Solanki & Haugan (1996) [S3] makes use of the symmetry 
properties of the Stokes profiles produced by the Zeeman splitting of atomic spectral 
lines. In the absence of radiative transfer effects in a dynamic atmosphere net circular 
polarization, Stokes V, is antisymmetric in wavelength and net linear polarization aligned 
at 45° to the solar limb, Stokes U, is symmetric. Gravitational birefringence changes the 
phase between orthogonal linear polarizations and thus partly converts Stokes V into 
Stokes U and vice versa. However, U produced from V by gravitational birefringence 
still has the symmetry of V and can thus be distinguished from the Zeeman signal. 

Let A$ be the phase shift which accumulates at the central wavelength of a spectral 
line between Stokes V and U as light propagates from a point on the solar surface to 
the observer. Formulae for A$ as predicted by metric affine theories were already given 
in chapter 1. For Moffat's NGT (Moffat 1979 38 ) a corresponding expression has been 
published by Gabriel et al. (1991) |2*T]. Further, let subscripts 'a' and V signify the 
antisymmetric and symmetric parts of the Stokes profiles, respectively, and the subscripts 
'src' and 'obs' the Stokes profiles as created at the source and as observed, respectively. 
Then, 

£4,src _ ^a,obs sin A$ + t/ a)0bs cosA$ 
^s.src 14,obs sinA$ + [/ Sj0bs cosA$ 

^s,src ^s,obs cosA$ + ?7 Si0bs sinA$ 
Vk.src Va )0bs cosA$ + U afihs sinA$ 

Thus for observed symmetric and antisymmetric fractions of U and V Eqs. ()2.1j) and ()2.2j) 
predict the ratios U a /U s and V s /V a at the solar source. 

If the solar atmosphere were static these ratios would vanish (C/ a ,src = ^a ;S rc = 0), so 
that any observed U a or V s would be due to either A$ or noise: £7 aj0 bs = Va,srcsinA$, 
K,obs = ^s,srcsinA$. The solar atmosphere is not static, however, and consequently 
the Stokes profiles do not fulfill the symmetry properties expected from the Zeeman 
effect even for rays coming from solar disc centre, which are unaffected by gravitational 
birefringence. This asymmetry has been extensively studied, in particularly for Stokes V, 
which most prominently exhibits it (e.g. Solanki & Stenflo 1984 [64J, Grossmann-Doerth 
et al. 1989 jHS], Steiner et al. 1999 |H], Martinez Pillet et al. 1997 HE])- Although most 
profiles have V s /V a < 0.2, a few percent of V profiles exhibit V s /V a values close to unity, 
even at solar disc centre. Such profiles occur in different types of solar regions, e.g. the 
quiet Sun (Steiner et al. 1999 [66J), active region neutral lines (Solanki et al. 1993 |69j) 
and sunspots (Sanchez Almeida & Lites 1992 [ZD])- The magnitude of V s /V a decreases 
rapidly with increasing V a and profiles with V s /V a >l are all very weak. They are often 



(2.1) 
(2.2) 



2.1. TECHNIQUE 



31 



associated with the presence of opposite magnetic polarities within the spatial resolution 
element and a magnetic vector that is almost perpendicular to the line of sight, situations 
which naturally give rise to small V (e.g. Sanchez Almeida & Lites 1992 (20], Ploner et 
al. 2001 [H]). 

The observed Stokes U asymmetry is on average smaller than the V asymmetry. 
This is true in particular for extreme asymmertic values, i.e. (^a,obs/^s,obs)max "C 
(^s,obs/^a,obs)max- Since this relation also holds at cos 9 = ji = 1 where 9 denotes 
the heliozentric angle between the source on the solar surface and the line-of- sight, it is 
valid for source profiles as well. Thus, Sanchez Almeida & Lites (1992) point out 
that Stokes U retains ^a,obs/^sobs ^ 1 throughout a sunspot, although V$/V a > 1 is 
invariably achieved at the neutral line. The reason for the smaller maximum asymmetry 
lies in the fact that Stokes U senses the transverse magnetic field. Since velocities in the 
solar atmosphere are directed mainly along the field lines they generally have a small line- 
of-sight component when U has a significant amplitude. Sizable line-of-sight velocities 
are needed, however, to produce a significant asymmetry (Grossmann-Doerth et al. 1989 
|65j). Another reason for the smaller maximum U asymmetry is that, unlike Stokes V, 
it does not distinguish between oppositely directed magnetic fields. 

Thus it is not surprising that in the following analysis Stokes U provides tighter limits 
than Stokes V. Another reason is that due to the on average stronger observed V profiles 
asymmetries introduced in U (through gravitationally introduced cross-talk from V) are 
larger than the other way round. However, we also analyse Stokes V as a consistency 
check. 

In order to seperate the asymmetry produced by solar effects from that introduced by 
gravitational birefringence, one strategy to follow is to consider large amplitude Stokes 
profiles only. Another is to analyse data spanning a large range of fi values, since gravita- 
tional birefringence follows a definite centre-to-limb variation, as predicted by particular 
gravitation theories. Finally, the larger the number of analysed line profiles, the more 
precise the limit that can be set on A$. Better statistics not only reduce the influence of 
noise, they are also needed because for a single profile gravitational birefringence can both 
increase or decrease V s /V a and U a /U s . The latter may become important if the source 
profiles are strongly asymmetric. Thus a small observed V s /V a or U a /U s is in itself no 
guarantee for a small gravitational birefringence. However, since almost all source pro- 
files are expected to have V s /V a <C 1 and U a /U s <C 1, on average we expect gravitational 
birefringence to increase these ratios. 



2.1.2 Profile difference technique 

The profile difference technique relies on the fact that A$ is expected to be a strong 
function of fi which is confirmed in the two concrete cases of NGT and metric afline 
theories (see chapter 1). This means that for any sufficiently large NGT charge £ & or 
equivalent metric afline parameter k a mixture of Stokes V and U profiles will be observed 
across the solar disc irrespective of the relative numbers and strengths of V and U profiles 
leaving the solar photosphere. Thus, irrespective of the value of < |V^ rc | > — < \U STC \ > 
for sufficiently large £ & or k < |V^b s | > — < |£/ b s | > will tend to zero. The averaging is 
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Figure 2.1: Top: | V^^s I — l^obsl vs - ^© an d Z 1 f° r l^srcl = (left) and \U SIC \ = (right). 
Bottom: I V^bs I — | t/ bs I averaged over all \i for the cases above. 

over all fi values and the total number of profiles is assumed to be very large. 

This effect is illustrated Fig. 12.11 In in the top of Fig l2.ll | V^^s I — l^obsl i s plotted vs. 
£q and fi for the extreme cases |V^ rc | = (Fig l2.ll top left) and \U SVC \ = (Fig l2.ll top 
right). Other combinations of |V^ rc (/^)| and \U ST c{^)\ give qualitatively similar results. 

The iV^bsl — I C/ bs I surface oscillates even more rapidly with increasing £ Q and with 
decreasing /1. Lines of equal iVJ^sl ~~ l^obsl are strongly curved in the £ Q — fi plane. 
These two points combine to lead to decreasing < |V^b s | > — < l^obsl > with increasing 
£q. This is shown in the bottom of Figs l2.ll left and right for the cases illustrated in 
the figures above. As expected the < |T^bsl > _ < l^obsl > vs - ^0 curves exhibit a 
rapidly damped oscillation around zero. This effect can be used to set upper limits on 
gravitational birefringence if the observations exhibit a < | V^bsl > — < l^obsl > that 
differs significantly from zero. As is later shown, this is indeed the case. 
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2.2 Observations and data 

Two sets of data have been analysed in the present work. They are described below. 
2.2.1 Data obtained in 1995 

Observations were carried out from 7 - 13th Nov. 1995 with the Gregory Coude Telescope 
(GCT) at the Izana Observatory on the Island of Teneriffe. For the polarimetry we 
employed the original version of the Zurich Imaging Polarimeter (ZIMPOL I), which 
employs 3 CCD cameras, one each to record Stokes I±Q, I±U and 7± V simultaneously 
(e.g. Keller et al. 1992 [721). 

The recorded wavelength range contains four prominent spectral lines, Fe I 5247.06A, 
Cr I 5247.56A, Fe I 5250.22A & Fe I 5250.65A. Three of these spectral lines are among 
those with the largest Stokes amplitudes in the whole solar spectrum and are also un- 
blended by other spectral lines (Solanki et al. 1986 [73] )• Blending poses a serious 
problem since it can affect the blue-red asymmetry of the Stokes profiles. By analysing 
more than one such line it is possible to reduce the influence of hidden blends and noise. 
Nowhere else in the visible spectrum are similar lines located sufficiently close in wave- 
length that they can be recorded simultaneously on a single detector. Also, compared 
to other lines with large Stokes amplitudes the chosen set lies at a short wavelength. 
This is important since the influence of gravitational birefringence on line polarization 
is proportional to 1/A. The sum of the above properties make the chosen range almost 
uniquely suited for our purpose. 

In order to image all 4 spectral lines of interest onto a single CCD we introduced 
reduction optics between the image plane of the spectrograph and the detectors. They 
produce an image-scale reduction by a factor of 3.2. The final spectral resolving power 
A/AA corresponded to 210'000. The spatial scale corresponding to a pixel was 1.13" (or 
860 km on the sun). However, the spatial resolution of the data is limited by turbulence 
in the Earth's atmosphere, the so-called seeing. This varied somewhat in the course of the 
observing run, so that the estimated angular resolution of the observations lies between 
1.1" and 3". 

The modulator package, composed of 2 piezoelectric modulators oscillating at frequen- 
cies of around 50'000 & lOO'OOO and a linear polarizer (glass prism), was placed ahead of 
the entrance slit to the spectrograph, but was nevertheless (unavoidably) located after 2 
oblique reflections in the telescope. Oblique reflections produce cross-talk between Stokes 
parameters, i.e. they partially convert one form of polarization into another. Since we 
are trying to observe, or at least set limits on "cross-talk" between Stokes U and V due 
to gravitational birefringence we took some trouble to reduce the instrumental cross-talk 
to the extent possible. A first step was the choice of the telescope. With only two 
oblique reflections, whose relative angles change only slowly in the course of a year, the 
GCT is relatively benign compared to most other large solar telescopes. Secondly, a 
half-wave plate was introduced between the two oblique reflections. Sanchez Almeida et 
al. [JjJ EE] have pointed out that a half-wave plate at that location should, under ideal 
circumstances, completely eliminate all instrumental cross-talk. To test the efficiency 
of the half-wave plate in suppressing instrumental cross-talk between Stokes Q, U and 



34 



CHAPTER 2. SOLAR OBSERVATIONS 



V we first carried out a series of observations of a sunspot near the centre of the solar 
disc both with and without a half-wave plate introduced in the light path. Such tests 
were necessary since the half-wave plate available at the GCT is not optimized for the 
observed wavelength. Note that at solar disc centre (// = 1) gravitational birefringence 
disappears, so that we test for instrumental cross-talk only. The half-wave plate was 
indeed found to significantly reduce instrumental cross-talk. Remaining cross-talk was 
removed during data reduction using a numerical model of the telescope that includes an 
imperfect half- wave plate (adapted from a model kindly provided by V. Martinez Pillet). 
The parameters of the model were adjusted slightly using the observations of the sunspot 
umbra close to disc centre. We estimated that the residual cross-talk after this procedure 
is at the level of a few percent. Since Stokes V, Q, U profiles generally have amplitudes 
of 0.1I C or less, the influence of the cross-talk is of the same order as the noise, which is 
roughly 1 — 2 x 10^ 3 / c , where I c is the continuum intensity. Photon noise is by far the 
largest contributor to this noise level. At this level instrumental cross-talk caeses to be 
of concern for our analysis. 

The Zimpol polarimeters are unique in that they combine CCD detectors with a very 
high modulation frequency (> 50'000 Hz) and hence preclude distortion of the Stokes 
profiles and cross-talk between them due to seeing fluctuations. This again improves the 
accuracy of the profile shapes and hence the accuracy of our results. 

A total of 106 recordings were made at different locations on the solar disc in an 
attempt to cover a large range of /i homogenously. Particular emphasis was placed on 
observations close to the limb since gravitational birefringence is expected to be largest 
for such rays. 

Since only a single sunspot was present on the solar disc during the observing run most 
recordings refer to faculae and network features, i.e. magnetic features with lower Stokes 
Q, U, V signals. In Fig l2.2l a sample Stokes /, Q, U, V spectrum of a facular region near 
the solar limb is plotted. The 4 analyzed spectral lines are identified. These data were 
fully reduced following the tedious, but well-tested procedures described by Bernasconi 
(1997) [ZH]. 

2.2.2 Data set of March 2000 

In order to improve the statistics and the [i coverage a second observing run was carried 
out in March 2000 with the Coude - Gregory Telescope in Locarno, Switzerland. This 
telescope is almost identical to the GCT on Tenerife and the parameters such as spectral 
resolution, noise level etc. are very similar to those of the 1995 observations. 

The second generation, ZIMPOL II polarimeter (Gandorfer & Povel 1997, Povel 1998 
jZHEHj) was employed for the polarization analysis and data recording. It simultaneously 
records three of the four Stokes parameters, either Stokes J, Q, V or /, U, V on a single 
CCD detector chip. Observations in these two modes were interlaced, such that alternate 
exposures record Stokes J, Q, V and J, U, V, respectively. Exposures of the same Stokes 
parameters were then added together to reduce noise. Thus the final data set consists of 
all four Stokes parameters. The only differences with respect to the recordings made in 
1995 are that the number of spatial pixels is reduced and that the noise level of Stokes V 
in the newer data is a factor of a/2 lower than of Stokes Q and U, whereas Stokes Q, U 
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Figure 2.2: Sample spectrum of Stokes /, Q, U and V. 
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and V had the same noise level in the earlier recordings. Due to the superior modulation 
scheme implemented in ZIMPOL II Stokes Q and U achieve a noise level of 10 _3 / c after 
roughly the same exposure time during the observations made in 2000 as during the 
earlier campaign. 

These observations were carried out on the day of the equinox, at which time the two 
mirrors producing oblique reflections of the beam ahead of the modulator package are 
oriented such that their polarization cross-talk cancels out. Hence for these observations 
the instrumental cross-talk is essentially zero and no further treatment of the data for 
this effect is required. 

The Sun was very active at the time of these observations with many active regions 
harboring sunspots and faculae present on the solar disc. Since active regions generally 
give larger amplitude Stokes signals we concentrated on observing them. A total of 7 
exposures were made. The typical seeing during these observations was estimated to be 
2 — 3", while the spatial pixel size was 1.13". 




Figure 2.3: SOHO MDI intensity image from 22nd March 2000. 
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2.3 Data analysis 



Each exposure gives us the line profiles of the 4 spectral lines in Stokes I,Q,U and V at 
a set of 94 (128 in the 1995 data) positions on the solar disc. Once the reduction and 
calibration procedure is completed we select from a given frame those spectra for further 
analysis for which the S/N ratio for either Stokes U or V in at least one of the four 
spectral lines is above 12 in the 1995 data and above 15 in the 2000 data. 

This criterion gave us a total of 4480 profiles from the four lines and were analysed fur- 
ther. The non-orthogonal wavelet-packets smoothing scheme of Fligge & Solanki (1998) 
|79*j was employed to enhance the S/N ratio by a factor of 1.5 - 2 without significantly 
affecting the profile shapes. 

Then a set of parameters was determined of all Stokes profiles of all 4 spectral lines. 
Of these parameters only the signed amplitudes of the blue and red wings, % and a r 
(i.e. of the blue and red Zeeman cr-component) are of relevance for the present study. In 
Fig. 2 we plot a Stokes U and V profile of the Fe I line at 5250. 22A. % r (F) and % T (U) 
are indicated in the figure. Using these we can form the symmetric and antisymmetric 
parts of the Stokes V and U profile amplitudes, Vs = (% + a r )/2, V& = (% — Or)/2, 
U s = (% + a r )/2 and U a = (% — a r )/2, respectively, which enter Eqs. (|2.1j) and (j2.2jl . 
In the following all V and U values and parameters are normalized to the continuum 
intensity, although this is often not explicitely mentioned. 
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Figure 2.4: Simultaneously measured profiles of Stokes V and Stokes U of the Fe I line 
at 5250.22A. 
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Figure 2.5: Measured amplitude asymmetries for Stokes U and Stokes V. 
2.3.1 Stokes asymmetry technique 

A measure of the asymmetry of a Stokes profile is given by the ratio SV = V s /V a , 
respectively 5U = U a /U s (e.g. Solanki & Stenfio 1984 |Mj, Martinez Pillet et al. 1997 
In Fig. 12.51 we plot these quantities vs. \V a \ and \U S \, respectively. Each point in 
these plots refers to a Stokes profile of the Fe I 5250. 65A line. Although our observations 
cover a range of /i values, Fig. 12.51 is very similar to corresponding figures based on data 
obtained near solar disc centre (Grossmann-Doert et al. 1996 [HZj, Martinez Pillet et al. 
1997 [HE]), where the influence of gravitational birefringence is expected to be negligible. 
For large amplitudes (V a , U s ) the relative asymmetry (V s /V a ,U a /U s ) is small, while for 
weaker profiles it shows an increasingly large spread. For the weaker profiles this spread 
is mainly due to noise as can be judged from the solid and dashed curves in Fig. 12.51 
which outline the la and 3a spread expected due to photon noise, respectively. The 
curves reveal that Stokes profiles with amplitudes (Va, Us) smaller than one percent of 
the continuums intensity are so strongly affected by noise that they are of little use for 
the present purpose. This leaves us with 1966 individual profiles for further analysis. 
In Fig |2.6l we plot a histogram of the number of these profiles as a function of \x. The 
distribution is uneven, being determined by the position on the solar disc of magnetic 
features at the time of the observations. 

The further analysis is made more complicated by the fact that a V s and a U a signal 
can be produced not just by gravitational birefringence, but also by radiative transfer 
processes acting in the dynamic solar atmosphere, as described in Sect. 3. To circumvent 
this problem we consider all profiles satisfying the criterion that |Va| or \U S \ > 0.01. For 
all these profiles | <5V^b s | < 0.7 and |<5£7 bsl < 0-6- A similar picture is also obtained at 
the center of the solar disc. Hence one way to limit £@ is to require |oV^ rc | < 1 and 
|<5£/ src | < 1 for all positions on the solar disc. This condition is strengthened by the fact 
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Figure 2.6: Histogram of the number of profiles as a function of /x. 

that I^V^bsl and |^t/ bsl decrease with decreasing /x, (Stenflo et al. 1987 [HI], Martinez 
Pillet et al. 1997 jEH|), whereas gravitational birefringence increases towards the limb, so 
that one would expect exactly the opposite behaviour if gravitational birefringence had 
a significant effect on the Stokes V or U profiles. 

In Fig l2.7l we plot the maximum (C/a,src/£4,src) va hie predicted for each of the 4 
spectral lines, based on all analysed data, vs. £ 2 Q . The horizontal line represents the value 
log(£/a,src/£4,src = 1), a limit above which this ratio is not observed at fj, ~ 1. Clearly, 
as t% increases (t/ a ,src/t4,src)max initially remains almost equal to (f/ a , bs/ f/ s,obs)max, 
but begins to increase for £Q>(50km) 2 , becoming > 10 at £ 2 & < (100 km) 2 and finally 
oscillating randomly around (f/a,src/^s,src)max 

of 100 - 1000. All 4 spectral lines exhibit 
a similar behaviour, implying that the influence of noise is very small. The largest effect 
of gravitational birefringence is exhibited by the two most strongly Zeeman split lines, 
Cr I 5247.56A, Fe I 5250.22 A, which also produce the largest Stokes V and U signals, 
while the line with the smallest splitting, Fe I 5250.65 A provides the weakest limit. 

It is in principle sufficient to limit £ Q by requiring that none of the observed spectral 
lines has (?7a,src/£4,src) > 1 within the range of allowed £ Q values. This gives £ 2 Q < 
(57.1km) 2 . A limit obtained similarly from Stokes V is both larger l 2 & < (80.3 km) 2 and 
less reliable, since we cannot completely rule out V^ iSrc /Va iSrc > 1 to be present, although 
we expect such profiles to be very rare, among large amplitude profiles. 

Fig |2.8l shows the analogous plot for metric affine theories. As in the above case for 
Moffat's NGT it can be seen that the strongest limit for gravity-induced birefringence is 
exhibited in the Cr I 5247.56A and Fe I 5250.22A line. We get k 2 < (1.1 km) 2 from Cr I 
5247.56A and k 2 < (0.9 km) 2 from Fe I 5250.22A. 
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Figure 2.7: Maximum of (C/ a ,src/^s,src), on a 
lines vs. £|. 
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Figure 2.8: Maximum of (C/ a)Sr c/C/ S)Src ), on a logarithmic scale, values for all 4 spectral 
lines vs. the metric affine parameter k. 
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Figure 2.9: Fraction of profiles with (£/ a , sr c/£4,src) > 1 for NGT (top) and MAG (bot- 
tom). 



An alternative test is to determine the fraction of profiles with (£/a,src/^s,src) > 1 (Fig l2.9|) 
Fig |2.9l reveals that initially no U profile satisfies the criterion, above £ 2 & = (57.5 km) 2 1.5 
% of the profiles does. This number keeps increasing with £ & , before finally oscillating 
around 70% at large £®. Thus 10% of all data points have (£/ a ,src/^s,src) > 1 for £ 2 & = 
(74.0 km) 2 , 20% for l 2 Q = (79.4 km) 2 . We are not aware of any solar observations of Stokes 
U with SU > 1 for which instrumental cross-talk is negligible (see Sanchez Almeida & 
Lites 1992 [7U], Skumanich et al. 1990 |E2|). £% < (80km) 2 is thus a very conservative 
upper limit. The lower picture shows the analogous plot for metric affine theories. Also, 
initially no profile satisfies the criterion, above k 2 = (1.33 km) 2 1.7 % of the profiles 
does. This number keeps increasing with £ Q , before finally oscillating around 70% at 
large k. Thus 10% of all data points have (?7a,src/^s,src) > 1 for k 2 = (1.95 km) 2 , 20% 
for k 2 = (2.43 km) 2 . k 2 < (2.5km) 2 is thus a very conservative upper limit. 
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Figure 2.10: Plot of the observable (Stokes V, Stokes U) mixture for initial phase dif- 
ferences of Ott, 0.l7r, 0.257T and 0.47T. Note that A$ values bigger than 0.57T give cyclic 
results. 

2.3.2 Profile difference analysis 

We now apply the technique outlined in Sect. 2.1.2 to our data. Due to the limited 
number of profiles and their irregular distribution over // (see Sect. 2.3.1) any limit 
on gravitational birefringence will be less tight than what is achievable with ideal data 
presented in Sect. 2.1.2. In order to be able to compare the present results with literature 
values we first set limits on the ^©-parameter in Moffats NGT, since earlier work has 
concentrated on constraining this theory. 
In FigEH3we plot 



< l^obsl > - < l^obsl > 

< lK,b S l > + < l^obsl > 



vs. £ Q , (2.3) 



for different initial phase differences A<3> between the orthogonal modes of line Fe I 
5250.65A. The averaging has been done over the fi values at which observations are 
available. This line is chosen, since it gives the tightest limits on £ 2 & . The thick horizon- 
tal line represents the value obtained from observations. Obviously above l 2 & = (178 km) 2 
the curve obtained from theory always lies below the observed value, hence ruling out 
such £q values. Note that it is sufficient to test for l 2 & < (305 km) 2 , since this upper limit 
has been set using independent data and another technique by Solanki & Haugan (1996) 



Similarly we can also set a limit on the coupling constant k of the metric affine 
theories. Using this technique we obtain k 2 < (0.69 km) 2 measured in line Fe I 5250. 65A. 
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2.3.3 Brief history of constraints on £q 

Up to this point we have presented our results with respect to sharp constraints on the 
coupling constants £ 2 & and k 2 . In order to get an estimate of the quality of these results 
compared to previous upper limits we will present here a brief overview on the history of 
constraints on £ 2 & achieved with a variety of different techniques. 

The first upper bound on the value of the NGT charge £ 2 Q was given by Moffat himself 
in 1982 [39J. Based on a new determination of the quadrupole moment of the sun he 
concluded that this result would lead to a deviation of 1.6% from Einstein's prediction 
for the precession of the perihelion of Mercury. Moffat claimed that NGT could fit the 
measured precession with the new quadrupole moment by using a value of £ Q = 3. 1-10 3 km 
close to the upper bound of £ Q < 2.9 • 10 3 km obtained in 1982 by using an average value 
on the available Mercury data. 

The next significant improvement was made by Gabriel et al. in 1991 |22j . By using 
for the first time the NGT prediction of gravity-induced birefringence they presented a 
quantitative prediction for the phase difference between the orthogonal modes. Conse- 
quently, birefringence leads to a depolarization of the Zeeman components of spectral lines 
emitted from extended, magnetically active regions on the sun. Since the solar-physics 
data which were available at that time limited the extend of such depolarization, they 
could only conclude that the Sun's antisymmetric charge must be less than (535 km) 2 . 

This result was based on extremely conservative assumptions regarding the polar- 
ization of light emitted by the observed solar feature and an inapplicable requirement 
regarding the determination of the magnetic filling factor, which is defined as the frac- 
tion of the aperture covered by magnetic field. In 1995 Solanki and Haugan jHS] improved 
the limit on £ 2 & by using more realistic values of the parameters that characterize the so- 
lar polarization source and a refined analysis procedure which yielded an upper bound of 
£% < (305 km) 2 . 

The current work represents the state of the art with respect to upper limits on grav- 
itational birefringence induced by the nonsymmetric field of the sun. Since the physical 
consequences of the nonmetric NGT description relevant for our purpose are proportional 
to £q this new constraint on the Sun's NGT charge is three orders of magnitude smaller 
than the previous limit by Solanki and Haugan and 7 orders of magnitude smaller than 
the value favoured initially by Moffat. This gives a further significant restriction to the 
viability of NGT. 
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Figure 2.11: Left figure: Sketch of the Solar Probe trajectory near the sun as seen from 
Earth. Right figures: Expected phase shifts for a signal emitted at 3R & from the sun. 
Top: Optical wavelength. Bottom: Soft X-ray. 





2.4 Possible tests using the Solar Probe spacecraft 

The possibility of setting strong limits on gravitational birefringence with the tests which 
we have outlined so far suffered from the major drawback that they all depend strongly 
on our knowledge of solar magnetic features. For example we have assumed that none 
of the observed spectral lines has (t/ a ,src/^s,src) > 1 which is also based on current 
observations with a certain spatial resolution. At the moment nobody can say if this 
statement will hold if the instrumental resolution could be highly increased in the future, 
although simulations of solar magnetoconvection suggest that £/ a ,src < U ssrc remains 
valid even if the resolution is improved by an order of magnitude. 

Taking this into account, it is obvious that drastic improvements of our limits for C. 2 & 
could be achieved by using an artificial source for the polarized radiation with well defined 
properties, placed in close vicinity to the sun. Such a source could possibly be provided 
by the Solar Probe spacecraft [83J, which is part of NASA's mid-term plans. Designed 
mainly as a mission to explore the Sun's Corona, its trajectory (sketched in Fig. a) lies 
in a plane perpendicular to the ecliptic with a target perihelion distance at 3-4 solar 
radii. This close perihelion could provide a unique opportunity to look for gravitational 
birefringence if it is possible to receive a well defined polarized signal from an instrument 
onboard the spacecraft when it passes near the solar limb. 

To see, if we could expect a reasonable gravitational effect in the signal, we simply 
have to take the case \i = in the phase shift formula (jl.57j) and replace the solar radius 
R Q with (R p + 1) • Rq where R p is the perihelion distance of the spacecraft in units of 
the solar radius. This leads to 
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Figure 2.12: Expected phase shifts for a signal emitted at 3R Q from the sun when the 
spacecraft passes in front of the solar disc on an alternative trajectory. The wavelengths 
are the same as in Fig l2.11l 

Of course this result depends strongly on the wavelength of the signal. The shorter the 
wavelength the bigger is the possible phase shift that we could expect. This is illustrated 
in Fig. (|2. lib ), which shows the phase difference in units of 7r for a range of £ Q values 
with R p = 3. Optical wavelength is used in the upper diagram and X-ray in the lower. 
One can see that for optical wavelength a gravitational effect on the polarization of the 
signal would be extremely difficult to observe under realistic conditions with a certain 
noise level. Whereas for wavelengths of the order 10 _10 m, we could expect a reasonable 
transformation of linear to circular polarization and vice versa for £ 2 Q values which are 
clearly below our current upper limit. However, at the moment it is out of reach from the 
technological point of view, to place a source of polarized X-rays in the near vicinity of 
the sun. Nevertheless it is worth to keep this promising possibility for very precise tests 
of the Einstein equivalence principle in mind. 



2.4.1 Alternative Spacecraft Trajectories 

Interestingly the planned spacecraft trajectory is also a very suitable one with respect 
to possible tests of gravity-induced birefringence. This is due to the fact that when the 
spacecraft reaches its perihelion position the phase shift is maximized by the mimimal 
distance to the sun, i.e. highest gravitational potential, and also by the maximum distance 
that the signal could travel through the gravitational field of the sun. 

If the craft reaches on a different trajectory its closest position to the sun which is 
also closer towards earth this would mean a shorter distance for the signal to travel and, 
hence, a smaller total phase shift. This is shown for a spacecraft passing in front of the 
visible solar disc at distance R to the sun. The phase shift of a signal with wavelength A 
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emitted at this position is given by 

( 3n - ( 1 3 \ 3 R \ 

1 ^ " X [l6Rj ~ 2 {2JRJ+W + 4R 2 (R 2 o + R 2 ) J ~ 8^I arctan i^J • 

(2.5) 

The detailed calculation is given in Appendix B. Fig l2. 121 shows the result for the same 
wavelength as in Fig J2. Ill and /x = 0. The cumulative phase shift for the alternative 
trajectory is smaller by a factor of 15 although the distance to the sun is the same in both 
cases. The crucial point is that in the first case the signal covers a bigger distance within 
that region of the gravitational field where it is strong enough to induce a reasonable 
phase shift. This relevant distance has an extend of approximately 5-6 solar radii. A 
different trajectory where the spacecraft passes the visible limb behind the sun as seen 
from earth would give no improvement since the phase shift which accumulates from 
behind the sun up to a position above the poles has an opposite sign compared to the 
phase shift of a signal which is emitted above the pole, propagating towards earth. This 
can be seen from equation (B.8) in Appendix B which gives the phase shift for a signal 
running from the surface up to a distant point R. If the signal runs in the opposite 
direction we have to switch the integration limits which gives a negative sign for the 
phase shift. 



2.5 Discussion and Conclusions 

Using two techniques (the Stokes asymmetry technique and the new profile difference 
technique) we have improved previous limits on £q given by Solanki & Haugan (1996) 
jHS] by nearly one order of magnitude. It will be difficult to set much tighter limits on 
gravitational birefringence than those found here using solar data in the visible spectral 
range. To obtain a significant improvement on the basis of solar data one would need to 
observe at shorter wavelengths. The line at the shortest wavelength that is strong enough 
to provide a hope of detecting Stokes U and V at sufficient S/N is Lya at 1216A. The 
maximum gain that one could expect relative to the current analysis is a factor of 

^!* = §0 4.32 . (2.6) 
A LyQ 1216 

Although the idea of a possible utilization of future space missions for enhanced new tests 
is very appealing because of its high potential for further improvements of the current 
limits on £ 2 & and k 2 , the technological realization of such an experiment is currently out of 
reach. For setting stronger limits on gravitational birefringence or, if birefringence really 
has a physical relevance, for having the chance of a direct detection of this effect it is 
therefore more promising to proceed the investigations with more compact astrophysical 
objects. This will be the content of the following chapters. 



Chapter 3 

Magnetic White Dwarfs 



The conjecture that gravity-induced birefringence is most distinctive for high gravita- 
tional potentials suggests itself from the structure of the phase shift formulas in NGT 
and MAG. For this reason the next logical step in this work is to continue the analysis 
with suitable, compact objects like magnetic white dwarfs. These stars provide a versatile 
tool for testing predictions of nonmetric theories of gravity because of their Megagauss 
magnetic fields and high surface gravity. Currently 65 or ~ 5% of all known white dwarfs 
are classified as isolated magnetic stars with field strengths in the range 3 x 10 4 — 10 9 G 
with the field strength distribution peaking at 1.6 x 10 7 G [87 . 

However, within this small class only a fraction of stars is of interest for our purpose. 
In addition to wavelength resolved polarimetric data we also need information about the 
mass and the radius of an object as well as measurements of the magnetic field structure 
to set strong limits on gravitational birefringence. With these restrictions we are left 
with just four magnetic white dwarfs which display high levels of circular polarization 
between 8 and 22 percent: Grw +70°8247, REJ0317-853, PG2329+267 and 40 Eridani B. 
The degree of circularly polarized light of wavelength A, {V/F) ^ s , reaching the observer 
from any of these stars is to first order simply a function of the strength of the local 
longitudinal magnetic field [85 a and the coupling constant k 2 . Since gravitation theories 
based on a metric-affine geometry of space-time predict a depolarization of light emitted 
from extended astrophysical sources we have looked for the largest k 2 that predicts values 
of {V/F)^ larger than or equal to that observed. This sets an upper limit on k. In the 
case of Grw +70°8247 such a limit was already set by Solanki et al. 86J based on NGT 
predictions. We redo the analysis using both an improved model of the magnetic field 
distribution of the star, as well as the predictions of metric affine theory, which had not 
been predicted earlier. 

However, we encounter a self-consistency problem since the relevant source properties 
like masses, inclination angles, etc. we use have been determined with models which, 
of course, neglect the possible influence of birefringence and, so, need not be valid. We 
circumvent this problem by assuming worst-case properties of the source that minimize 
depolarization caused by gravity-induced birefringence. 
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3.1 Polarization Modelling Technique 

The search for a possible influence of gravity-induced birefringence on astronomical sig- 
nals, i.e. the modelling technique which is used to set strong upper limits, depends of 
course heavily on the amount of appropiate information that we can gain on the source of 
polarized radiation. This becomes decisively clear in the transition from the techniques 
which we have used to evaluate the solar observations of the last chapter to the magnetic 
white dwarfs. Here, a crucial point is basically played by the spatial resolution of the 
observations, since the individual calculation of the phase shift A$ as a function of /x for 
a given light source requires, at least theoretically, an infinite resolution of the source. 

So, due to the finite spatial resolution of observations, the effect the phase shift has 
on a light ray depends on whether one observes a pointlike source, e.g. small sunspots, or 
an extended source like white dwarf magnetospheres. In the case of a pointlike source, all 
light received from it suffers the same phase shift A$(/x p ). Introducing Stokes parameters 
to describe polarized light, with Stokes Q defined to represent the difference between 
linear polarization parallel and perpendicular to the part of the stellar limb closest to the 
source, one therefore finds a crosstalk between Stokes U and Stokes V. This crosstalk is 
such that although the observed values U \y S and V \^ s differ from the values emitted by 
the source, U STC and V S tc, the composite degrees of polarization remains equal: (U^ s + 
^obs) 1 ^ 2 = (^src + Ksrc) 1,/2 - ^ an extended source covering a range of /x values is observed 
then light emitted from different points suffers different phase shifts and, so, adds up 
to an incoherent superposition. Summing over the different contributions, using the 
additive properties of Stokes parameters yields a reduction of the observed polarization 
relative to the light emitted from the source: {U^ s + V^) 1 / 2 < (Z7| rc + V^,) 1 / 2 . Since 
the stellar sources of polarization we consider in this chapter are spatially extended, any 
observed (nonzero) degree of polarization provides a limit on the strength of gravity- 
induced birefringence that the star's gravitational field could induce. 

It is generally agreed that the polarized radiation from white dwarfs is produced at 
the stellar surface as a result of the presence of Megagauss dipolar magnetic fields [§HJIIIH|. 
To first approximation the flux of net circularly polarized light at wavelength A emitted 
toward the observer from the surface is directly proportional to the strength of the line 
of sight component of the magnetic field at the stellar surface. Therefore, this net flux 
can be written as 



where J(/x) is the intensity at wavelength A. The weightfunction w(/x) describes the 
strength of the vertical magnetic field component at position /x. Since our analysis en- 
compasses different white dwarfs with different types of dipole field geometries, w(/x) is 
calculated using an oblique dipolar rotator model which is introduced in the next section. 

To define a degree of circular polarization, we have to divide (j3.1|) by the total stellar 
flux emitted to the observer at wavelength A 




(3.1) 




(3.2) 
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The function also describes the limb darkening of the star. Since the surface of a 
white dwarf cannot be resolved the limb darkening cannot be measured directly. Because 
the broadband spectrum of Grw +70°8247 is well represented by blackbody radiation 
[100 and radiative equilibrium models |133^ 1101] it is reasonable to assume a simple law 
like the one describing the directly observed solar limb darkening. For this, we have 
chosen to use [103 

/x</,) 1 + Gu - l)g + (n 2 - l)h , (3.3) 



= 1) 

with 

<g + h< 1 

for Grw +70°8247. Since this kind of a linear interpolation is of a very general form 
we have taken the limb darkening model of Grw +70°8247 as a prototype for the white 
dwarfs within our sample. 

In practise, we have evaluated the degree of circular polarization (V\ >S r C /F\) numer- 
ically for increasing values of k* and, in the case of Grw +70°8247, of I*. In contrast 
to the earlier work of Solanki, Haugan and Mann (HH] (1999) we assume that the light 
emitted by the source is polarized proportional to B, i.e. to w(/i). 



3.2 Oblique dipolar rotator model 

The model we use was already described by Stift in 1975 |104j from whom we adopted 
most of our notations. 

Starting with the observer's system where the z-axis is defined by the line-of-sight, a 
surface point on the visible hemisphere of the star is given by z. This vector is transformed 
into the rotation system of the star by the matrix 



( 1 







cos i sin % 
- sin % cos % 



(3.4) 



which is followed by a rotation around the rotational z-axis 



/ cos0 sin0 
— sin <p cos <p 

V ooi 



(3.5) 



The position of the dipole system relative to the rotational system is uniquely described 
by the three Eulerian angles, e.g. by 

S = S&S X (3.6) 
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Figure 3.1: Geometrical representation of the dipole- relative to the rotation system of 
the star with offsets. The rotation axis is inclined to the observer by an angle %. In 
general, the offset dipole is tilted against the rotation axis as described in Fig. 3. 2. 



where S% and S x are generating rotations around the x-axis and S$ around the z-axis. 
Since the only measurable quantity in observations is the resulting tilt angle f3 between 
the dipole and rotation axis, the effect of the three Eulerian rotations can be summarized 
without loss of generality by a rotation with f3 around the x-axis. 

Given the offset coordinates of the dipole in the rotation system by x, we obtain the 
coordinates of the surface point z relative to the dipole system by 



r = SpiS^SiZ - x) , (3.7) 

so that the field strength in the dipole system is given by 

B dp = -V(m ■ r/r 3 ) . (3.8) 

From this we get the field strength in the observer's system through 

Bobs = SfS^SpB^p (3.9) 

where Sj , S^, Sp are the transpose matrices of Si, S^, Sp. The longitudinal field com- 
ponent in the observer's system is then of course the B & Si2 component of ([3.9)1 . 
The effective field B e for the visible hemisphere is then given by 

B e = f [ ~B Z I dA I [ fldA . (3.10) 



B 2 denotes the line-of-sight component of B & s , while the limb darkening is described, as 
above, by the empirical function I. 



3.3. GRW+7CP8247 
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3.3 Grw +70°8247 

Within the class of magnetic white dwarfs, Grw +70°8247 has played the role of a 
"Rosetta stone" in the analysis and interpretation of their spectra. Although the first 
measurement of circular and linear polarization in the spectrum of a white dwarf jHH] 
gave a strong indication for a magnetic field, the presence of broad absorption features 
in its blue spectrum remained unexplained until the 1980s. A solution of this puzzle was 
only given after detailed calculations of hydrogen energy levels at field strengths B > 100 
MG jHH EDI EU E2] so that, as a result, all spectral features could consistently be ex- 
plained in terms of stationary lines, i.e. lines whose wavelengths vary only slowly with 
the field strength, of atomic hydrogen in a centered dipolar magnetic field with a polar 
field strength of ~ 320 MG and a viewing angle (angle between the line of sight and the 
magnetic axis) % = 0° - 30° HSU Ell- 
in this way it was possible to determine the physical parameters of the stellar atmo- 
sphere and of the magnetic field from a best fit of synthetic to observed spectra. The 
observations can be best explained by an effective temperature of 14,000 K |133^ HOOj 
which implies a radius of 0.0076 R & (determined from photometry and parallax distance 
measurements), where R & is the solar radius, and a mass of about 1.0 M , where M 
is the solar mass |102j . Since the characteristics of the polarization curve have remained 
invariant for almost ~ 25 yr now 87 j we can infer a very long rotation period of much 
more than 25 years. This slow rotation is fitted into the obique rotator model by setting 
(3 = 0, implying a dipole aligned with the rotation axis. Thus at a given point in the 
observer system with polar coordinates (r,9,(f)), the Cartesian field components in the 
case 2 = 0° are simply given by 

B x = 3Bd sin 9 cos 9 cos 0/2r 3 , (3.11) 
By = 3Z?d sin 6 cos 9 sin <p/2r 3 , (3.12) 
B z = B d (3cos 2 9 - l)/2r 3 , (3.13) 

where B^ denotes the surface polar field strength for the case of a centered dipole. In 
our case, the z-axis is also the line of sight. Since we focus on measurements of circular 
polarization, only the longitudinal Zeeman-Effect is relevant for our purpose. 

Polarization measurements of Grw +70°8247 yielded a level of — 6±0.25% in the visible 
spectrum at 449 nm jHH EH EZj . Since the constraints on k"l and ll are stronger for larger 
observed degrees of polarization and shorter wavelengths we make use of Hubble Space 
Telescope spectropolarimetry in the ultraviolet, which revealed high levels of circular 
(12%) and linear (20%) polarization. This was measured between 130 and 140 nm and is 
related to the absorption feature at 134.7 nm jHZj. To be conservative, we assume a large 
absolute error of 1% on these measurements, and use a degree of circular polarization of 
11% at 134.7 nm. 

In 1999 Solanki, Haugan and Mann jHH| have used the NGT prediction to set strong 
upper limits on ll of the magnetic white dwarf Grw +70°8247 with a very simple model 
of the magnetic field distribution. Although it became clear in the meantime that NGT 
suffers from technical problems that render it mathematically inconsistent, we also carry 



52 



CHAPTER 3. MAGNETIC WHITE DWARFS 




Figure 3.2: Left figure: Projection of the conelike movement of the magnetic pole during 
the rotational period onto the celestial plane for (3^0. Right figure: Relative orientation 
of the line-of-sight to the dipole axis and the rotational axis. Adopted from Burleigh, 
Jordan and Schweizer [107J. 

out the analysis based on Eq. ()1.57|) for Grw +70°8247, in order to have a quantitative 
comparison between the old, unrealistically simple and the improved model incorporating 
an oblique rotator of the magnetic field. 

The results of our analysis can be seen in Fig. ()3.3j) and Fig. (p£ 4)1 Our search for the 
largest values of k* and Z* compatible with the observed degree of cicular polarization 
at 134.7 nm yields the constraints /c 2 < (0.074 km) 2 and Z 2 < (4.1km) 2 for a viewing 
angle i = 0°. This limit on Z 2 is to be compared with the former limit Z 2 < (4.9 km) 2 . 
These results were obtained for limb darkening coefficients (g, h) = (0, 1), which provide 
the most conservative constraints. Neglecting limb darkening yields /c 2 < (0.06 km) 2 and 
Z 2 < (3.7 km) 2 compared with the old constraint Z 2 < (4.6 km) 2 . For other values of the 
(g, h) pair the constraint on and Z* falls between these extremes. 

For increasing viewing angle % to the dipole axis the maximum degree of circular po- 
larization decreases, due to the dipole field geometry. From this it follows that constraints 
on k* become stronger for larger i, although the dependence turns out to be weak. E.g. 
for i = 30°, we obtain the tiny improvement /c 2 < (0.072 km) 2 with limb darkening and 
kl < (0.058 km) 2 without limb darkening. Hence the most conservative estimate is for 
i — 0, which we adopt. 
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Figure 3.3: Observed degree of circularly polarized light for increasing values of from 
NGT, normalized to full intensity 1 = 1. Left figure: Dipole inclination 2 = 0°. Right 
figure: Dipole inclination i = 30°. The different curves are for different limb darkening 
parameters. The horizontal line represents the observed value of 11%. 




Figure 3.4: Observed degree of circularly polarized light for increasing values of from 
metric-affine gravity, normalized to full intensity 1 = 1. Left figure: Dipole inclination 
% = 0°. Right figure: Dipole inclination i = 30°. The different curves are for different 
limb darkening parameters. The horizontal line represents the observed value of 11%. 
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3.4 RE J0317-853 
3.4.1 Introduction 

RE J0317-853 is certainly a highly unusual object which sets several records within the 
class of isolated magnetic white dwarfs. Discovered in 1995 as an extreme-ultraviolet 
source during the ROSAT Wide Field Camera (WFC) all-sky survey by Barstow et al. 
[108J (hereafter B95) the analysis of B95 revealed an unusual hot white dwarf with an 
effective temperature of ~ 50000 K and an exceptionally intense dipolar magnetic field 
of ~ 340 MG. B95 also reported nearly identical photometric and polarimetric variations 
of 725.4 s. Since the most plausible explanation is that these modulations are due to 
the rotation of the star, this would lead to the fastest rotation period ever measured 
so far among the isolated magnetic white dwarfs which have typical periods of several 
hours. Finally, the close proximity of the ordinary, DA white dwarf LB09802, separated 
by only 16" allowed B95 identify RE J0317-853 as a member of a double degenerate pair, 
which then enabled them to estimate the mass of RE J0317-853 at 1.35 M Q , close to 
the Chandrasekhar limit, which makes it the so far most massive known isolated white 
dwarf with a corresponding radius of only 0.0035 R & . They found that their data are best 
matched by a dipole model with a negative offset along the dipole axis of 20% viewed at 
an angle of 60°. Polarization measurements were not reported by B95. 

Ferrario et al. (hereafter F97) modelled phase-averaged spectropolarimetric data of 
RE J03 17-853 in 1997 and obtained a dipole model with a polar field strength of Bd = 450 
MG and a negative offset of 35% of the stellar radius along the dipole axis. Their best fits 
were derived for mean viewing angles to the dipole axis of 30° — 60°. Furthermore, they 
reported variations in the wavelength-averaged circular polarization data at a period of 
725 ± 10 s which confirmed the results from B95. More important for our work is that 
F97 also found circular polarization in the continuum up to a maximum of 8% in the 
wavelength range from 5600 A - 5800 A. To be conservative, we adopted for our purpose 
the maximum value of 5800 A and as for the B95 data the maximum limb darkening 
coefficients g — 0, h — 1. 

The so far most detailed investigation regarding the magnetic field geometry of RE 
J0317-853 was published by Burleigh, Jordan and Schweizer in 1999 |107j . Using phase- 
resolved, far-UV HST Faint Object Spectrograph spectra they obtained the best fit for 
a dipole with Bd = 363 MG and additional offsets perpendicular to the z-axis: x q = 
0.057, y Q = —0.04 and z q = —0.22. In this model the rotation axis is viewed at an 
angle of i = 50° to the observer while the relative angle between rotation and dipole axis 
was found to be (3 = 42°, leading to visible surface field strengths between 140 and 730 
MG. In addition they reported the limb darkening coefficients g = 0.3 and h = 0. In 1999 
Jordan and Burleigh (hereafter JB99 [_ 1 1 0J ) reported a peak level of circular polarization 
at 22% integrated over the wavelength range between 3400 A and 7000 A measured in 
spectropolarimetric data taken with the AAT (Anglo- Australian Telescope). This is up 
to now the highest level of circular polarization ever recorded for a magnetic white dwarf. 
Although it seems very likely from the data in JB99 that the real wavelength range in 
which the high level of circular polarization was measured is between 5600 A and 5800 
A this conjecture was so far not confirmed by the authors. 
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Longitudinal Magnetic Field Strength in MG 
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Figure 3.5: Successive rotational phases of RE J0317-853 in steps of 0.257T, beginning 
with <p = 0.257T (top left) to <fi = 2ir (bottom right). The dashed line markes the stellar 
equator whereas the solid line shows the projection of the magnetic equator on the stellar 
surface. The cross marks the position of the rotation axis and the dot the position of the 
magnetic pole. 

3.4.2 Birefringence analysis 

Due to its small radius and high degree of circular polarization, RE J03 17-853 is a very 
suitably object concerning limits on gravitational birefringence. Unfortunately the wave- 
length range in which the polarization was measured by F97 as well as by JB99 lies in 
the optical regime and additionaly is up to now not very well isolated in the case of the 
22% measurement of JB99 jllOj . So, in order to get here conservative and reliable upper 
limits on gravitational birefringence we use 7000A which marks the upper end of the 
measured spectral range. For comparison we also present our estimates in the case of 
5800 A as the upper wavelength end. 

For our analysis of the JB99 data we also need information about the magnetic field ge- 
ometry at that rotation phase, when the highest degree of polarization is emitted towards 
the observer. Since the highest degree is obtained for the highest average longitudinal 
magnetic field strength over the visible hemisphere we have plotted in Fig. ()3.5j) successive 
rotational phases of RE J0317-853 in steps of 0.257T based on the reported offset values 
x oS = 0.057, y Q = —0.04, z fi = —0.22, limb darkening coefficients g = 0.3, h = from 
()3.3j) and the angles i = 50° and j3 = 42° by Burleigh, Jordan and Schweizer |107j . The 
longitudinal field strength on the surface is colorcoded whereby the blue color denotes 
positive and the red color negative polarity. It can be seen, that in phase <ft — 0.257T 
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Figure 3.6: Integrated longitudinal (B z ) magnetic field strength of the visible hemisphere 
for successive rotational phases of RE J0317-853. 

the different polarities nearly cancel each other so that the disk averaged field strength 
is approximately zero. During the progression of the rotation the average field strength 
increases up to = 7r and then subsequently decreases. This is also quantitatively shown 
in Fig.flHH). 

So, by using the oblique rotator model, described in Sec. 3.2 we obtain an upper limit 
of kl < (0.038 km) 2 for a wavelength of 7000A and k\ < (0.036 km) 2 in the case of 5800 
A. For both wavelengths we have used the limb darkening coefficients g = 0.3 and h — 0. 
The F97 data yielded a result of kl < (0.05 km) 2 for i = 30° and kl < (0.051km) 2 for 
the case of i — 60° and g — 0, h — 1 with a polarization level of 8% at 5800 A. 

The table below summarises the various upper limits for k obtained for the data from 
B95, F97 and BJS99. 





i 




k (km) 


F97 (8%) 


30° 


(0°) 


0.05 


x E = °> VoS = ,ZaK = -°- 35 


60° 


(0°) 


0.051 


BJS99 (22%) 


50° 


42° 


0.038 


x oS = 0.057, y oS = -0.04 ,z oS = -0.22 






0.036 



Tab. 3.1: Upper limits on k for the F97 and BJS99 data. 
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Figure 3.7: Upper panel: Percentage of circularly polarized light, present in the Ha 
spectrum of PG 2329+267. Lower panel: The normalized spectrum. From: Moran, 
Marsh and Dhillon (1998) [TTTj . 



3.5 PG 2329+267 



The magnetic nature of the white dwarf PG 2329+267 was discovered in 1995 by Moran, 
Marsh and Dillon by observing the characteristic Zeeman splitting of Balmer lines |lllj . 
The subsequent analysis revealed a centered dipole inclined at i = 60° ±5° to the observer 
with a polar field strength of appoximately 2.3 MG. Rotation of the star was not reported, 
so that we take f3 = 0. The mass was determined at w 0.9 M which means that 
PG 2329+267 is more massive than typical isolated white dwarfs. The radius was not 
measured so far, so that we calculated it, using the mass-radius relation given by Weinberg 
Asuming a helium star we got R = 0.0156 R Q . This radius would decrease if we 
asume heavier elements than helium which, in turn would lead to a stronger gravitational 
birefringence, so that this is certainly a conservative limit on R. The results of the circular 
spectropolarimetric measurements in the Ha line are shown in Fig. (j3.7j) . One can see very 
clear the S-shaped profile in the upper panel, indicating a magnetic field, with a peak level 
in the a components of approximately 10% at 6530 A. As in the case of Grw +70°8247 
we have calculated the observable degree of circular polarization for increasing values of 
k and for different limb darkening parameters. For i = 55° we found k\ < (0.25 km) 2 if 
limb darkening is neglegted and /c 2 < (0.31km) 2 for maximum limb darkening. In the 
case of i = 65° we obtain k\ < (0.23 km) 2 and fc 2 < (0.28 km) 2 , respectively. 
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Figure 3.8: Successive rotational phases of 40 Eridani B in steps of 0.257T, beginning with 
<p = 0.257T (top left) to (p — 27r (bottom right). The poles of the grid match those of the 
magnetic dipole. 



3.6 40 Eridani B 

40 Eridani B is one of the most famous white dwarfs and represents together with Sirius B 
the first known examples of this class of degenerate stars. The mass and radius of 40 Eri B 
were determined by Koester and Weidemann in 1991 |113j and more recently by Shipman 
et al. HEI with the result of M = 0.501 ± 0.011M© and R = 0.0136 ± 0.00024i? o . Its 
magnetic field was measured in 1999 by Fabrika et al. |114j using time-resolved Zeeman 
spectroscopy. From their data they concluded the presence of a centered dipole with an 
average field strength of several kG. The rotational axis was reported to be inclined at 
i ~ 90° to the observer with a rotation period of about 4 hours. The magnetic axis 
inclination to the rotational axis is about the same, (3 ~ 90°. Fig. (j3.8j) shows successive 
rotational phases of 40 Eri B using the given parameters. The circular polarization 
measurements reported by Fabrika et al. |114j in the Ha line core (see Fig l3.9l on the 
next page) yielded a level of 10% at ~ 6564 A. As for REJ0317-853 we have taken the 
magnetic field geometry of the rotation phase <p = tt with the highest average surface 
field strength for our analysis. With the given parameters we obtained an upper limit 
k\ < (0.35 km) 2 neglecting limb darkening and fc 2 < (0.44 km) 2 with maximum limb 
darkening. 
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Figure 3.9: Top panel: Polarization obtained from individual spectra of the zero polar- 
ization "cross-over point" at 0.57T or 1.5ir. Bottom: Two spectra of different polarization 
in the region of the Ha line core. The magnetic shift is clearly visible. Middle: The 
result of the substraction of these spectra, i.e. Stokes V. From Fabrika et al. (2000) 

ma. 
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3.7 Comparison of the results 

After we have obtained several upper limits on k 2 for various objects (including the sun) 
it is interesting to ask, whether we can find a correlation between them. For this purpose 
we have listed in Tab. 3.2 the masses and radii in solar units of the investigated objects, 
together with the corresponding highest and lowest upper limits on k 2 obtained from 
our analysis and their mean values in the fourth column. The last column shows the 
corresponding ratio $ between the Schwarzschildradius r$ = 2GM/ c 2 of the star and its 
physical radius R, i.e. $ = 2GM/(Rc 2 ). 





M/Mq 


R/Rq 


k 2 (km) 2 


k 2 (km) 2 


$ = 2GM/Rc 2 (lO" 4 ) 


Sun 


1 


1 


0.69 - 1.95 


1.32 


0.042 


40 Eri B 


0.5 


0.013 


0.35 - 0.44 


0.395 


1.614 


PG 2329 


0.9 


0.0156 


0.23 - 0.28 


0.255 


2.42 


GRW 


1 


0.0076 


0.058 - 0.074 


0.066 


5.52 


RE J0317 


1.35 


0.0035 


0.036 - 0.051 


0.0435 


16.2 



Tab. 3. 2: Relevant data of considered objects. First and second column: Mass and Radius in 
solar units. Third column: Range of upper limits on k 2 obtained from our analysis. Mean 
values k 2 in the fourth column. The last column gives the ratio between Schwarzschild and 
physical radius xlO -4 . 

The first thing which attracts attention is that the limits on k 2 decrease for increasing $ 
of the stellar objects. This correlation is graphically shown in Fig. 13.101 on the next page 
where we plot from Tab. 3.2 the mean value of k 2 for each object against <£>, marked 
by solid diamonds. The best fit to this distribution can be found with gaussian linear 
regression method. This yields 

A; 2 ($) = 1.363982 • (exp(-$ a7 ) + 0.04 $-°- 18 ) . (3.14) 

The agreement between the values /c 2 ($) of the fittet curve and the measured k 2 are 
quite good as one can see in Tab. 3.3. Therefore, ()3.14j) can be interpreted in that 
way, that it gives approximately those upper limits on k 2 which one could expect for 
celestial bodies of mass M and radius R by using the polarization modelling technique 
used in this chapter. For example an object of about the mass and the size of the 
earth (M E = 5.97- 10 24 kg, R = 6.378 • 10 6 m) would give k% < (1.771501 km) 2 , while a 
supermassive black hole similar to the one suspected in the center of our galaxy with 
M BH = 2.5 • 1O 6 M and a Schwarzschildradius of r s = 2GM/c 2 = 2.969 • 10 9 m would 
give k 2 BH < (0.010396 km) 2 . Finally we note that this would also be the upper limit for 
any other object which has reached its Schwarzschildradius, since k 2 (§) depends only on 
the ratio rs/R*- 
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Figure 3.10: Mean value of k 2 plotted against the corresponding ratio $ in units of 10~ 4 . 
The distribution can be fitted by a smooth curve, given in the text. The different ranges 
of the k 2 values are of about the same size as the symbols and therefore omitted. 





k 2 ($) 


k 2 


Deviation from k 2 in % 


Sun 


1.320012 


1.320 


9 • 10" 4 


40 Eri B 


0.387055 


0.395 


2.01 


PG 2329 


0.259636 


0.255 


1.81 


GRW 


0.090102 


0.066 


36.5 


RE J0317 


0.034262 


0.0435 


21.23 



Tab. 3. 3: Comparison of /c 2 (<E>) values from the fitted curve with the measured mean values k 2 . 
The last column gives the relative deviation in percent of the fitted from the measured values. 
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3.8 Discussion and Conclusions 

We have used spectropolarimetric observations of the white dwarfs Grw +70°8247, RE 
J0317-853, PG 2329+267 and 40 Eridani B to impose new strong constraints on gravita- 
tional birefringence. To improve earlier investigations of Grw +70°8247 we employed a 
dipole model for the magnetic field structure which leads to, on average, an 18% sharper 
limit. 

However, for all white dwarfs we have considerd our constraints on k* and Z* strongly 
depend on the sophisticated form of the weight function w(fi). At the moment, u>(jit) 
only gives very conservative results from a simple relation between the magnetic field 
strength at a certain position on the visible stellar disc and the emitted degree of cicular 
polarization. So, the next step to enhance w(fi) could consist of calculating numerically 
the theoretical degree of circular polarization for each point on the stellar disc by using 
simulations of radiative transfer in stellar atmospheres. A further, clear improvement on 
the constraint on k+ should be possible within this approach. 

The possibiliy, given by ()3.14|) to calculate approximate values for upper limits on 
k 2 depending on the ratio between the Schwarzschild and physical radius is certainly a 
useful result. Those values could be used as a guideline for the birefringence analysis of 
other objects. However, the constant numerical values in (j3.14j) possibly have to changed 
slightly when measured k 2 values of other objects has to be fitted into this scheme. 



Chapter 4 

Cataclysmic Variables 



The method which we have used up to now for setting strong upper limits on gravity- 
induced birefringence is very limited. As mentioned before, this is mainly due to the 
lack of a reliable atmospheric model for magnetic white dwarf stars which could provide 
realistic estimates for the emitted degree of polarized radiation as a function of field 
strengths and opacities. 

In this chapter we circumvent this restriction by focussing our attention on magnetic 
white dwarfs which are members in close interacting binary systems. These are systems 
where two stars move around the common center-of-mass in a close orbit so that matter 
can fall along a ballistic stream from one star to the other. An important subclass of these 
cataclysmic variables are the AM-Her type binary systems where the white dwarf posses 
a strong magnetic field so that matter is prevented from forming an accretion disc and, 
instead, impacts directly on the stellar surface near the magnetic poles. During this pro- 
cess highly polarized cyclotron radiation is emitted whose rotational modulation provides 
a unique opportunity not only for setting upper limits on gravity-induced birefringence 
but perhaps also for the very first direct observation of this effect. 

We therefore begin this chapter with a brief introduction into the basic physics of close 
interacting binary systems. The characteristics of the polarized radiation are determined 
by the various parameters like the temperature structure and height of the shock front, 
formed when the accreting matter hits the stellar surface, as well as by the electron 
number density of the heated plasma and the angular extend of the emission region. A 
realistic estimate of the emitted degree of polarization therefore depends heavily on the 
quality of the shock front model. With this underlying physical structure, the theory 
of radiative transfer of polarized radiation in the case of cyclotron emission is reviewed. 
Although the theory is well known for decades now, we add a new feature by modifying 
the resulting Stokes vector when we allow for gravity-induced birefringence. 

Using this theoretical framework we present phase dependent light and polarization 
curves for the AM-Her system VV Puppis which are compared with polarimetric data 
taken by Cropper & Warner in 1986 |135j . Although the observed asymmetries are 
in agreement with our numerical curves, the flat topped character seen in observations 
cannot be achieved numerically. However, by allowing for sufficiently strong birefringence 
we obtain a good conformity with the observations. 
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4.1 Basic model of interacting binaries 

Approximately 50 — 60% of all stars are assumed to be members of binary systems. By 
involving nearly the whole spectrum of stars from red giants to black holes as possible 
binary components, these systems often provide examples of astrophysical processes which 
are among the most exciting currently known. From a more practical point of view the 
most precise mass determinations in astronomy are available from binary systems where 
the period of revolution is precisely measurable due to mutual occultations. 

The physics of binary systems depends of course on the masses of the involved stars 
and, more importantly, on their spatial separation. If the distance between both compo- 
nents is big enough so that they can interact at every evolutionary stage only gravitation- 
ally without any related mass transfer, then each member of the binary will pass through 
its individual development without being disturbed by its companion. An example for 
this is Sirius with the white dwarf Sirius B as the second star. However, the situation 
becomes more interesting if the stars move so close to each other that tidal effects are 
no longer negligible. In this chapter we are interested in such systems which consist of a 
white dwarf as the primary component and an, evolutionary younger and more massive, 
secondary component which can be represented by a main sequence or a red giant star. If 
both stars are not too close, their initial spherical shapes becomes teardrop-shaped with 
the tip pointing towards the companion. Basically a mass transfer stream between the 
components can now be initialized in two ways, i.e. either by an evolutionary expansion 
of the secondary component or by a continuous decrease of the distance between the 
stars. In both cases, the secondary component begins to fill its Roche lobe (the volume 
surrounding an object within matter is gravitationally bound to it) and matter is trans- 
fered to the primary component via the inner Lagrange point LI. The mechanism of how 
exactly the distance between the stars could be further decreased is not fully understood 
at the moment. One possibility is that they already originated in close proximity or that 
due to friction of the surrounding material from both components the rotational energy 
of the system was partially converted into internal energy. 

Once the mass transfer has started, the angular momentum of the system prevents the 
stream from falling directly onto the white dwarf surface. Instead the material forms, in 
most cases, an accretion disc where the gas slowly spirals inwards while, again by friction, 
the gravitational energy is converted into thermal energy. Finally, the matter falls onto 
the white dwarf from the inner edge of the disc. However, this is not a continuous process 
since the material is mostly accumulated in the accretion disc until the disc becomes 
unstable and, then, is suddenly deflated accompanied by a huge brightening of the disc 
due to conversion of gravitational to thermal energy. Such an event is also known as a 
dwarf nova. Such close binary systems which interact by means of a mass transfer stream 
are called cataclysmic variables. 

For our purpose, the most interesting among them are those where the white dwarf 
posseses a ultrastrong magnetic field of the of the order 10<_B<60 megagauss (MG), 
which prevents the formation of a normal accretion disc and leads, instead, to a well- 
defined accretion flow from the secondary component to near the magnetic pole(s) of 
the white dwarf. These objects are often called AM-Her systems, after the prototype 
AM-Herculi discovered in 1923. 
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Figure 4.1: Artist impression of an Am-Her type system. 



The name "Polar" for AM-Hers and related systems was introduced by Krzeminski & 
Serkowski in 1977 |117j because of the strong and variable circular (~ 10%— 30% |135| I115j ) 
and linear polarization which is typical for these objects. Certainly this attribute makes 
the polars very appealing for utilizing them for possible measurements of gravity-induced 
birefringence. Therefore, the next sections will give some insight into the basic physi- 
cal processes responsible for the emission of polarized radiation. This, in turn provides 
the basis for a critical analysis of the observed polarization properties with respect to 
gravity-induced birefringence. 
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Figure 4.2: Schematic representation of the accretion pattern in AM Hers (adopted from 
Wickramasinghe & Ferrario [HTj). 



4.2 Mass transfer and shock models 



The basic model which has been developed for AM-Her type systems is briefly sketched in 
Fig J4.21 It is mainly based on X-ray observations |118j . modelling of the continuum and 
emission-line radiation from magnetically confined accretion funnels , and polarized 
emission from shocks [120 . 

For an ordinary, nonmagnetic primary component the matter streams from the sec- 
ondary star and accumulates in a disc before impacting on the white dwarf surface. On 
the other hand the detailed accretion mechanisms are much more complicated and even 
controversial (e.g. |121j ) if the primary component posseses a strong magnetic field. 
However, it is generally accepted that after the material has left the companion star near 
the inner Lagrangian point, a stream is formed which falls freely towards the white dwarf 
on a ballistic trajectory. If the magnetic field of the accreting object is strong enough 
(> IMG) the increasing magnetic pressure will begin to dominate the material flow at 
some distance from the white dwarf in a region called the coupling region. Here, the 
radial velocity component is reduced and the flow aquires first toroidal and then poloidal 
velocity components as the stream couples onto the field lines. From that point the ma- 
terial falls again freely towards the accreting object in magnetically confined funnels that 
connect the coupling region with the white dwarf surface. Finally a static shock above 
one or both magnetic poles slows and heats the infalling material before it settles on the 
surface with a temperature of kT « lOkeV. 
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4.2.1 The shock region 

The shock region can be regarded as the main source of emission in the optical and the 
X-ray wave band. Calculations of the circular polarization light curve presented in this 
chapter are therefore based on a sophisticated model of the emission region, first invented 
by Wickramasinghe & Megitt in 1985 123\. 

The first attempts to model the shock fronts of AM-Her type systems, given by 
Meggitt & Wickramasinghe in 1982 [124J, were based on the assumptions of constant 
temperature accretion columns. Although successful in explaining some of the gross 
properties like the continuum energy distribution, this model, based on cyclotron opacity, 
failed to predict the correct degree of polarization (the predicted values were a factor of 
about 2-3 higher than observed). The model which is discussed here can be regarded as a 
further development of this early model since it includes a temperature structure (shock 
front) as well as free-free opacity as a pure absorption process. 

Assuming that the ionized material falls freely towards the magnetic poles, it moves 
a radial distance r with the velocity 

!/(,■) = (—) =5.2X10»(-) (— ](-) , (4.1) 

where M and R are the mass and the radius in (cm) of the white dwarf respectively. 
Before the material settles down on the surface, the density is increased in a strong shock 
by a factor 4 and the velocity decreases by the same factor across the shock front. In 
the case that this shock is formed near the white dwarf surface, the shock temperature 
is given by 



3fim H GM 8 / M W10 



cm 



^if^ 6x 10 "tej(— r ■ (42) 

where \i denotes the mean molecular weight. The postshock electron density is 

^- 3 ' 8Xl °il0WT)(3" /2 (S 1/2 ^ 3 ' < 4 ' 3 > 

where F s is the specific accretion rate. In this context it is important to note that the 
scope of this model only allows accretion within a small fraction of the stellar surface 
namely on a point like structure. This restriction will later be removed by a more general 
approach that also allows for more extended shock regions. 

The bulk of the translational energy is carried by the ions which achieve a Maxwellian 
velocity distribution by ion-ion collisions at a distance h ion ~ V ion t ion , where t ion denotes 
the ion- ion collision time scale |125j . while the electrons are mainly heated by Coulomb 
collisions with the ions. The electron temperature usually differs from the ion temperature 
and both vary with height in the shock |126| I127j . As the gas then settles on the stellar 
surface it is cooled by bremsstrahlung and cyclotron radiation. In the case that cyclotron 
radiation dominates, the electrons cool faster than they can be heated by collisions with 
ions, with the consequence that ions and electrons have different temperatures, so that 
the gas in the "two-fluid" regime. However, it is important to note that the accretion 
flow in AM-Her systems is not continuous. Occasionally these systems drop to low states 
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Figure 4.3: Orientation of the dipole axis d and the spin axis l s relative to the symmetry 
axis l p of the emission region and the line-of-sight. 

of reduced brightness which could be explained by interruptions or reductions in the 
material stream. The reason for this accretion modulation is not fully understood at 
the moment, but it may be related to solar-type magnetic activity in connection with 
starspots of the secondary component |122j . 

4.2.2 Extended emission regions 

The idea that the emission regions are pointlike sources located at the diametrically op- 
posite magnetic poles is insufficient and often contradicted by observations. For example 
the asymmetries seen in polarization light curves of most systems or the wavelength de- 
pendence of intensity and polarization in AM-Her itself are hardly reproduced by these 
models. Therefore, Wickramasinghe & Ferrario 132] I133j introduced emission regions 
which extend across and above the stellar surface. Consequently, this new concept al- 
lowed for field spread, density and temperature structure, and for displacement from the 
magnetic poles. 

The field geometry is that of a centered dipole. Using a cartesian coordinate system, 
the spin axis l s is inclined with respect to the line-of-sight by an angle % in the y — z- 
plane, while the symmetry axis of the emission region l p is identified with the positive 
z-axis making an angle A with respect to the spin axis. The emission region itself 
extends over a spherical cap which subtends an angle 2a at the center of the star. In 
this coordinate system, the dipole axis has a polar angle Qb with l p and an azimuthal 
angle $b with respect to a line perpendicular to l p in the same plane as l p and l s . We 
can therefore identify this line with the (positive) y-axis of the cartesian system in the 
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direction l p x (l p x l s ). The relative orientations of the various axes are shown in Fig J4.3l 
For our purpose, the best agreement with observations was achieved by assuming 
uniform electron temperature T e and fixed electron number density N e . The emission 
region has a constant radial thickness H so that the characteristics of the emission region 
can be summarized by the optical depth parameter 

A H = 2.01 x 10 8 (#/10 8 cm)(Ayi0 lc W 3 )(3 x 10 7 G/B) (4.4) 

in the radial direction ESQ . Of course A# varies over the emission region due to changes 
in B. The magnetic field strength B in the above expression is evaluated at the magnetic 
pole. 



4.3 Cyclotron radiation and birefringence 

Since the cyclotron emission is highly polarized, the AM-Her systems provide an exellent 
opportunity to test for gravity- induced birefringence. For this purpose we have used 
a FORTRAN program, written by D.T. Wickramasinghe which calculates the Stokes 
vector of cyclotron emission in the 'point source' approximation of Wickramasinghe and 
Meggitt |123| using a fixed electron temperature T e and a fixed optical depth parameter 
A. This section will present the underlying mathematical formalism with the additional 
new feature of a gravity-induced alteration of the Stokes vector. Extended emission 
regions, introduced in the previous section, are easily obtained from this concept as a 
homogeneous composite of pointlike sources. 

Consider an electron in helical motion in a uniform magnetic field directed along the 
z-axis in a cartesian coordinate system. The energy emitted per unit solid angle Q making 
an angle 9 with the magnetic field and per frequency interval du is given by 

e 2 oo 2 00 

e Q duj dQ = - — F n8n(y) (4.5) 

27FC n=l 

where €q is also the coefficient of spontaneous emission, with 

F n = (cot 6 - /3||Cosec #) 2 J^«) + /?i^ 2 «) (4.6) 

where J n is the Bessel function of order n and J' n {n^) = dJ n (n^) / d(n^) |128j . Further, 
(5\\ and /3± are the parallel and perpendicular components of the dimensionless velocity 
f3 = v/c, so that we define £ = sin 8/(1 — (3\\ cos 6). The argument of the 8 function is 

y = nuj 7 - u(l - /3\\ cos 6) (4.7) 

where u c = eB/mc is the electron cyclotron frequency and 7 = (1 — (3 2 )~ l l 2 the usual 
Lorentz factor. From this one can easily see, that the radiation spectrum consists of 
spectral lines occuring at frequencies 

nuj c 

" = 7(l-Acosfl) ' (48) 



which also gives a relation between 7 and harmonic number n. 
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In order to describe a uniform plasma with temperature T and electron number density 
N, one has to introduce a relativistic Maxwellian distribution in ()4.5j) and integrate over 

P\\ 

61 dw dVt = — — — ^ ^ J F/exp(— n^)^ /nuj c )d(3\\ duo dVt (4.9) 

where the dimensionless electron temperature \i = mc 2 /kT serves as the argument of the 
modified Bessel function of the second kind, K 2 (n). 

To obtain the emission coefficients €q and ey for Stokes Q and Stokes V, one defines 
a coordinate system where the z-axis serves as the line-of-sight and the uniform magnetic 
field lies in the x-z plane, making an angle 9 with the z-axis. Then, the complex Jones 
vector for the electric field of the radiation of a single electron is proportional to (|128j) 

/ (cot 9 — /3||Cosec 9) J n (n£) \ 

-i(3 ± j'M) . (4.10) 

v ; 

The emission coefficients €q and ey for the Stokes parameters Q and V are obtained by 
replacing Fj in (|4.9jl with 

F Q = (cote - /?||Cosec 9) 2 J 2 n M) - (3 2 ± j' 2 H) (4.11) 

and 

F v = -2(cot 9 - /3||Cosec 9) J n (n£)(3±j' n (n£) (4.12) 

respectively, while ejj = 0. Since < £ < 1 the Bessel function J n and its derivative can 
be replaced with the Wild-Hill approximations |l()5j which are more easy to implement 
in a computer program 

JnW) - |^xp(n W ) K + 5033/n] _ 1/6 
V27m(l + w) n 

j'nW) - ^frl > 3 + 1 - 193 /^] 1/6 [1 - V5n 2/3 ] (4.14) 
V 27rn(l + w) n 



with u> = ^1 — £ 2 - The cyclotron absorption coefficients are given by Kirchhoff's law, 
since the electrons are in thermal equilibrium, i.e. 

ei = n cyc B^, e Q = qcycB^, e v = v cyc B^ (4.15) 

with the Rayleigh- Jeans law 

B^ = uj 2 kT/AiT 3 c 2 . (4.16) 

In contrast to the first model, given by Meggitt & Wickramasinghe in 1982 |124j the new 
scheme also includes free-free opacity as a pure absorption process. The corresponding 
opacities are 

u 2 {2uo A + 2u 2 u 2 c - 3co 2 co 2 sin 2 9 + u 4 c sin 2 9) 
Ks = 2WV^f? Uc (4 - 17) 
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Uf, 



ufa* sin 2 6(u 2 - 3a; 5 
2civ 2 {u 2 - u 2 ) 2 



2uj 2 ujuj c cos 6 

2^2 

c I 



CU! 2 (iU 2 — LU 2 ) 



(4.18) 
(4.19) 
(4.20) 



with the plasma frequency cu p = (AnN e e 2 /m) 1 / 2 and the approximately collision frequency 
v c = 3.63iV e T~ 3 / 2 ln(2.95 x 10 n T/u;). It is assumed that the opacity can be treated as a 
sum of cyclotron and free-free components, hence 



"eye "T /%, Q 

The transfer equation now becomes 



<?cyc "T ^ff , 



(4.21) 
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(4.22) 



with the Faraday mixing coefficients 



/ = (^/c^ c )cos^/(a; 2 M 2 -l) 

h = (ul/cu c ) sin 2 9/2(u 3 /ul - uj/uj c ) 



(4.23) 
(4.24) 



For uniform conditions these coefficients are constant along a ray. In order to solve the 
transfer equation ()4.22|) set 



m 
R 
A 



0.5(/ 2 + h 2 ), n = 0.5(g 2 + v 2 ), p = qf - vh, r = qh + vf, 

[(m + n) 2 - p 2 ] 1/2 , at = (m + n + R)/p, a 3 = l/a x , 

(n-m + R) 1/2 , fx=(m-n + R) 1/2 , 

(/ - ?ai)/A, h = (/ - ga 3 )//i 

—(h + vai)/X,c 3 = —(h + va 3 )/fi . 



(4.25) 
(4.26) 
(4.27) 
(4.28) 
(4.29) 



Setting all intensities equal to zero at s = 0, the Stokes parameters of the radiation after 
passing through a distance s are (see |124j ) 



//-B^ = 1 — (p/2R)(a 1 cosh As — a 3 cos^s) exp(— ks) 

Q/Bu = — (p/2R)(bi sinh As — 63 sin/xs) exp(— ks), 

U/Bu = (p/2i?)(cosh As — cos/is) exp(— ks), 

V/Bu = — {p/2R)(c\ sinh As — c 3 sin /is) exp(— ks) 



(4.30) 
(4.31) 
(4.32) 
(4.33) 



The question of how gravitational birefringence could be implemented in this scheme de- 
pends heavily on the shock height above the stellar surface where the polarized radiation 
is emitted. Since the phase shift A$ in (jl.filj) is proportional to 1/R where R denotes the 
radial distance from the stellar center, an emission region which is extended of approx- 
imately one stellar radius above the surface would require an integration process along 
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the shock height where a gravitationally modified transfer equation has to be solved for 
each plane parallel slab. An upper limit on the shock height is given in the 'single-fluid' 
approximation when bremsstrahlung is the dominating cooling mechanism. In this case 
the height h hr is simply proportional to the product between the free-fall velocity V s (the 
index 'ff', of course, must not be confused with the index for free- free absorption) and 
the specific cooling time t s for electrons due to bremsstrahlung. The explicit expression 
for h hl , given by Wickramasinghe & Ferrario jHZ] is 

7 / 10 16 cm~ 3 \ ( M \ ( 10 9 cm\ , . 

Assuming a white dwarf mass of M with a corresponding radius of 0.014 R & this yields 
a height of ~ 10 8 cm with iV e = 10 15 . Since we focus on a 'two-fluid' model when the 
electrons are cooled by cyclotron radiation, this means that we have for the cooling time 
^cy C "C t a and, therefore, h cyc <C h hI . For this reason it is justified to view the polarization 
as emitted from near the surface because h cyc <C 10 2 km is certainly a sufficient upper 
limit. We can therefore assume that birefringence acts directly on (|4.32j) and ()4.33|) 
according to 

/ cosA$ \ / U~ 



grav 



V„„ I \ cosA$ 




(4.35) 



Our objective is to include birefringence given by these equations in the calculations 
of synthetic polarization light curves from AM-Her type systems. In this way we can 
set sharp upper limits on this effect with a new method, independent from techniques 
presented in chapter two. 



4.4 Polarimetry of VV Puppis 

VV Puppis serves as a key system for our purpose of either detecting weak gravitational 
birefringence or setting strong upper limits on it by means of gravitationally modified 
polarization curves. We therefore present here some of the most important results that 
have been achieved for this system within the last decades and which provide the basis 
for the subsequent analysis. 

VV Puppis was discovered in 1931 by van Gent |129| as a faint (V = 14.5 — 18) 
variable with a period of 100 min. Sinusoidal variations in radial velocity of the H and 
He II line, correlated with the 100 min photometric variation led Herbig in 1960 to the 
conclusion that the observations could best be explained by means of a binary system 
where the emission lines originate on the brighter (primary) component |130j . Herbig 
also concluded that changes in the shape of the light curve could have its origin in 
variations of the emission area on the visible stellar surface which can be realized by self 
eclipses of the relevant area due to rotation. This picture was refined in 1977 by Tapia's 
discovery of strong linear and circular polarization with a maximum of ~ 16% |116j in 
each polarization state of the optical light curve, confirming the suggestion by Bond and 
Wagner (IAUC 3049) that this object is similar to the "original" AM-Her system. Two 
years later, Visvanathan and Wickramasinghe identified a series of absorption features 
in the spectrum with the 6th, 7th and 8th harmonics of cyclotron absorption in a nearly 
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Figure 4.4: Intensity and polarization light curves of VV Puppis. From: Copper & 
Warner (1986) |I55] . 

uniform magnetic field of ~ 3 • 10 7 G. This discovery renewed the general interest in 
astrophysical cyclotron emission models, since it was the first detection of resolvable 
cyclotron harmonics in the optical spectrum of a stellar object. However, these results 
had an uncertainty of one harmonic in the identification since the observed features did 
not have well defined cores. As a consequence the corresponding uncertainty in B was 
~ 20%. Concerning the visibility of cyclotron harmonics, numerical computations, based 
on the calculations presented in section 4.3 revealed that very special conditions are 
required to produce consistency with the observed spectrum |132j . For the shock front 
model related to VV Puppis one can therefore conclude a viewing angle i ~ 70° — 90° 
with respect to the magnetic field, low optical depths A ~ 10 5 and high temperatures 
of T ~ 10 keV in a magnetic field of ~ 3 • 10 7 G. The fact that VV Puppis is the only 
stellar object to show resolvable cyclotron harmonics in its optical spectrum is difficult 
to explain, although it might be possible that magnetic field broadening often renders 
the features undetectable |132j . 

Fig l4.4l shows phase dependend, wavelength averaged polarization curves of VV Pup- 
pis from Cropper & Warner (1986) |135j . As one can see, the light curves are strongly 
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Figure 4.5: Visualization of the longitudinal magnetic field components of VV-Puppis 
during successive rotational phases. The black spot which is visible from cf) « 0.75 to 
<p Pd 1.25 markes the position of an emission region with size a = 10° and displacements 
= 90°, Q B = 30° from the dipolar axis. 



modulated at the orbital period and the accreting pole is visible only for ~ 40% of the 
period. A common feature is the slow rise after eclipse egress and a more rapid decline be- 
fore eclipse ingress. Although linear polarization is present during the entire bright phase, 
with a second, lower peak at phase 0.7, the strong linear polarization peak of ~ 9% is 
expected to be caused by a magnetic field configuration that can be seen only briefly at 
a viewing angle of ~ 90° with respect to the line-of-sight. In this case linear polarization 
due to cyclotron radiation is strongest and would therefore explain the observed feature. 
The circular polarization reaches a maximum value of ~ 10% during the bright phase 
after a more gradual rise. The steeper descent of the circular polarization is followed by a 
sign reversal and a short dip, lasting 0.05 of a cycle. One of the most important results is 
that VV Puppis shows steady negative circular polarization throughout the entire eclipse 
phase, which could be explained by accretion onto a second region of the white dwarf 
[136J, that is visible at all rotation phases |137j . However, the negative polarization dip 
at eclipse ingress is probably related to the shape of the primary accretion region. This 
conjection is supported by synthetic light curves of only one accretion region as we will 
later show. 

On the basis of the light curves in Fig |4.4[ Cropper & Warner reported for the system 
geometry an inclination of i — 77° ± 7° and a magnetic colatitude of A = 155° ± 6°. 
In 1987 Wickramasinghe & Ferrario used i = 75° and 5 = 150° for the calculation of 
phase dependent light curves of a displaced accreting pole. Fig |4.5l shows a visualization 



4.4. POLARIMETRY OF VV PUPPIS 



75 



of the longitudinal magnetic field component during a rotation cycle of VV Puppis. The 
emission region with size a = 10° has an azimuthal and polar angle of = 90° and 0_b = 
30°, respectively, marked by the black spot. One can see that in this configuration the 
emission region merely grazes the limb during the bright phase, which is extraordinarily 
useful for our purpose, since gravitational birefringence is most pronounced for sources 
located at the stellar limb. Further, the asymmetry in the linear pulse and the sign 
reversal in circular polarization, prior to eclipse ingress can now be explained by means 
of the system geometry. Calculating the angles 9 between the field directions in the spot 
and the line-of-sight one finds, that the angles are mostly greater than 90° at eclipse 
ingress and mostly less than 90° at eclipse egress. The observed reversal in the sign of 
circular polarization occurs just shortly before eclipse ingress when 9 > 90° while the 
stronger linear pulse as well as the maximum intensity occurs at 9 ~ 90° as expected 
from the beaming properties of cyclotron radiation. 

4.4.1 Gravitationally modified lightcurves 

The programm which we have developed for calculations of gravitationally modified po- 
larization curves is based on a FORTRAN routine written by D.T. Wickramasinghe which 
calculates cyclotron emission in the " point source" approximation of Wickramasinghe & 
Meggitt |123j for fixed optical depth parameter A and fixed electron temperature. Since 
the point source approximation is insufficient to account for most of the observed polar- 
ization curve characteristics as explained in section 4.2.2, our code describes an extended 
emission region as a composition of ~ 3 x 10 2 point sources, depending on the size of 
the emission region. This emission region is placed on a rotating sphere according to the 
coordinate system described in sect. 4.2.2, so that the polarization values (circular and 
linear) taken at each rotational phase yields the final polarization curve. 

We first present light curves without gravitational modification in order to have a 
comparison to the observed curves of Cropper & Warner |135j . Fig l4.fi! shows the phase 
dependent polarization and light curves for two different models. While both models 
have the same inclination i = 75° and magnetic colatitude A = 150° for the position 
of the spot, the left column shows the curves for polar cap size a = 5° and a dipole 
axis orientation with 9b = 10° and = 90°. The curves in the right column have 
a = 10°, 9b = 30° and $b = 90°. We take note that these curves are in good agreement 
with the curves presented by Wickramasinghe & Ferrario in 1988 |133j (hereafter WF88) 
for extended emission regions using similar parameters. At this point it is important to 
note, that these curves (including those in WF88 |134j ) require an additional constant 
background of unpolarised light in order to reproduce observations at least approximately. 
This constant, unpolarized continuum probably comes from a component of the binary 
system which is always visible, like the accretion stream, white dwarf photosphere or 
something similar - up to now the genuine origin is unknown as well as how to correctly 
estimate it. Currently the only restriction is placed by the observed polarization levels 
and, so, the continuum has to be introduced by hand to get the observed levels right. 

Since our code only allows for emission from one pole, the circular polarization curve 
does not show a constant negative polarization throughout the hole period which is seen 
in the observations. Basic characteristics like the asymmetry in the linear pulse and the 
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negative dip in circular polarization prior to eclipse ingress stand comparison to what is 
observed by Cropper & Warner. Nevertheless we note that the observations shows a more 
gradual rise in intensity and polarization whereas the rapid declines are in agreement with 
our models. 

These "pure" curves (without any birefringence influence) could possibly better fit- 
ted to observations by a different geometrical shape of the emission region. While the 
present curves are based on a cylindrical symmetric source, a more elongated, arc-like 
structure could perhaps account for the observed gradual rise. Observational hints for 
arc-like emission regions has been presented by Beuermann et al. in 1987 jT38i and 
also by Cropper |139j . We will come to this point again after the discussion of the 
gravitationally modified polarization curves. In accordance with WF88 we emphasize 
that the characteristics of the polarization curves depend sensitively on for a given a 
and 9b- For instance, for = 0° we obtain no polarization reversals while for $_b = 180° 
we observed two symmetrical reversals at eclipse ingress and eclipse egress. 
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Figure 4.6: Phase dependent circular polarization curve of VV Puppis for an emission 
region centered at latitude A = 150°. The orbital inclination is % = 75° and uo/u c = 6.2. 
The curves in the left column correspond to a = 5° and 8b = 10°, $_b = 90°, while the 
right column corresponds to a = 10° and 9 B = 30°, = 90°. 
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One of the main problems that these curves and also the curves presented by WF88 
have is the sinusoidal form of the circular polarization and intensity, in contrast to what 
is observed. For this reason it is certainly interesting to see how gravity-induced bire- 
fringence modifies the circular polarization curve. Fig l4.7l on the next page shows the 
phase dependent circular polarization using the same system parameters as in Fig |4.6l 
and M* = 0.4 M Q , R+ = 0.014 R Q but this time additionally with increasing values of 
the metric-afline coupling constant k. As expected, for k 2 = the curves are identical 
with the unmodified polarization regime while for a sufficiently big k 2 , the light is almost 
completely depolarized. In between these two "boundary" states, one can adjust k 2 to 
a value (0.13 km) 2 < k 2 < (0.14 km) 2 so that the circular polarization curve is nearly 
flatten as seen in observations. At the same time, the negative dip which arises prior 
to eclipse ingress is only marginally influenced by birefringence, although it but also 
decreases somewhat with increasing k 2 . 

In order to see if this range for k 2 is in agreement with upper limits, derived with the 
methods of the chapter 3, we have made two consistency checks. By using the oblique 
dipolar rotator technique we got an upper limit of ky V Pup < (0.35 km) for VV Puppis 
without limb darkening and ky V Pup < (0.42 km) with maximum limb darkening and an 
observed polarization level of 10%. These limits are similar to those of 40 Eridani B from 
the previous chapter, i.e. k 2 mEriB < (0.35 km) and k 2 mEriB < (0.44 km), respectively 
with M/^EriB = O.5M and Rmehb = 0.0136 R & comparable to Myv Pup = O.4M and 
Rvv Pup — 0.014 R Q . The other white dwarfs of chapter 3 have higher masses and/or 
smaller radii so that their upper values on k 2 are correspondingly lower. The same is 
true for the k 2 value of the sun. However, one has to note that this dipolar rotator 
technique assumes that the emitted degree of polarization is direct proportional to the 
longitudinal magnetic field strength, whereas we have seen that this relation becomes 
much more complicate in the case of cyclotron radiation. It is therefore more reliable to 
use the "plain" model proposed by Solanki, Haugan and Mann in 1999 [HE] which assumes 
that the polarized radiation is emitted homogeneously over the visible stellar surface. 
Since this is a rather conservative restriction the resulting upper limit is congruously 
k 2 < (0.4 km) 2 . Nevertheless it is important that the range of A; 2 -values for VV Puppis is 
consistent with the limits from the dipolar rotator model as well as with the plain model. 

Since the amount of Stokes U in the linear polarization signal is always < 5% the 
linear polarization curve is hardly affected by gravity-induced birefringence which only 
acts on Stokes V and Stokes U as one can see in Fig J4.81 Nevertheless, the first smaller 
peak becomes sharper and more pronounced for a = 5° while for a = 10° the peak is 
washed out. 
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Figure 4.7: Gravitationally modified, phase dependent cicular polarization curves. The 
physical parameters are the same as in Fig Mfil Upper curves: a = 5°, 6 b = 10°, = 
90°. Lower curves: a = 10°, 9 B = 30°, $ B = 90°. 
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Figure 4.8: Gravitationally modified, phase dependent linear polarization curves. The 
physical parameters are the same as in Fig l4.(il For comparison the same /c-values as for 
circular polarization have been plotted. Upper curves: a = 5°, 6b = 10°, $_b = 90°. 
Lower curves: a = 10°, 6 B = 30°, $ B = 90°. 
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4.5 Comparison of the results 

Finally, we want to check how the limits on k 2 for VV-Puppis, obtained by lightcurve 
fitting and the polarization modelling technique from chapter 3 fit into the empirical 
scheme, given by ([3.14| . For this purpose, Fig. 14.91 shows again the curve from Fig. 13.101 
additionally with the ky VP values, obtained in this chapter. The limit which is obtained 




Figure 4.9: A; 2 ($) curve of chapter 3 together with the VV Puppis limits. The upper star 
(0.6|0.385) marks the limit obtained by the polarization modelling technique. The lower 
star (0.6|0.15) marks the limit from lightcurve fitting. A better fit to the VV Puppis 
value is provided by the dashed curve. 

by the polarization modelling technique deviates slightly from the curve. Nevertheless 
the fitting can be improved with a minimally different function 

k 2 (§) = 1, 026262 • (exp(-$ L1 ) + 0.1 $" a18 ) , (4.36) 

which yields the dashed line in the above figure. The bigger deviation of the lower point 
(0.6 10. 15) from the curve is not surprising, since its intention is not to serve as an upper 
but as a lower limit on k 2 . Certainly more datapoints measured at different celestial 
bodies are needed to see which curve can best approximate the upper limits on k 2 . 
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4.6 Conclusions 

The main difference between the results of this chapter and those of chapter 3 is obvious. 
While for "ordinary" magnetic white dwarfs we merely set upper limits on k 2 we now 
additionally have to assume a lower bound on k 2 in order to link the presented theory of 
cyclotron radiation in AM-Her type systems with observations. However one has to be 
very careful when making an interpretation of these results since the underlying model 
of interacting binary systems is afflicted by many uncertainties: 

1. The fact that a constant unpolarized continuum background has to be introduced 
by hand to reproduce observations is unsatisfactory. A certain aspect of the physics 
of binary systems is therefore probably not understood and although it is rather 
unlikely that this aspect could lead to a modification of the present lightcurves 
similar to the gravitationally modified curves, it cannot be completely ruled out. 

2. More importantly, Wickramasinghe & Ferrario [140 have shown that an arc-shaped 
emission region is also able to replace the sinusoidal shape seen in the circular 
polarization curve by a flat top. 

At a first glance, the second point seems to discard any claim for a lower bound on k 2 . 
Indeed, we are here confronted with two alternative models which try to explain the same 
observed polarization properties with different approaches. Currently there is no reliable 
and convincing possibility to judge in favour of one or the other approach. Although the 
agreement of the conventional approach with observations seems to be very appealing, the 
new question which now arises is, if the assumption of an arc-shaped emission region for 
VV-Puppis is justified any longer, taking into account the new results from gravitational 
birefringence. One of the main arguments supporting the arc-shaped emission regions is 
the flat top character of the calculated circular polarization curve as well as the correctly 
predicted degree of polarization - both also provided by gravity-induced depolarization. 
Concerning future projects it is therefore highly desirable to look for specific predictions 
of the conventional model which do not involve polarized radiation and result, 
not affected by birefringence. 

However, as an important conclusion of this chapter we can say that, as long as 
gravitational birefringence cannot be excluded, the interpretation of polarized radiation 
from compact astrophysical objects is always afflicted with many uncertainties since the 
real nongravitational source properties may be hidden behind a birefringent curtain. 



Chapter 5 
Future Projects 



Since astronomical tests of the Einstein equivalence principle with respect to gravitational 
birefringence are a relatively new area of research, we are left with several promising but 
so far unexplored ideas which have the potential not only for setting sharp upper limits on 
this effect but also for possible direct detections. We present in this chapter a selection of 
two elaborated but not yet finished examples where gravity-induced birefringence could 
either serve as a possible important contribution for a satisfying interpretation of current 
data or could show up as a new effect in atomic spectra under certain conditions. 

The first example is concerned with the problem of circular polarization from active 
galactic nuclei (AGN). Under the reasonable assumption that synchrotron radiation is 
the main source for polarized radio emission, the degree of circular polarization measured 
in those AGN is on average an order of magnitude higher than expected from synchrotron 
theory. The suggested classical solution in terms of circular repolarization is not free from 
contradictions which leads to the conclusion that a gravity- induced crosstalk between 
linearly and circularly polarized light could provide an important contribution for solving 
this puzzle. The second example is concerned with gravitational red shift measurements 
in nonmetric theories of gravity. While in metric theories of gravity the ticking rates of 
(atomic) clocks are independent on their internal structure, i.e. the red shift is universal, 
this universality is violated for clocks which are affected by nonmetric gravitational fields. 
Defining three atomic clocks by different transitions of hydrogen one gets three different 
red shift predictions which could be tested in astronomical spectra of compact objects or 
under well defined laboratory conditions. 

In the last part of this chapter we briefly discuss two further ideas for promising 
future projects. We suggest to compare the polarization properties of multiple images 
generated by gravitational lenses to look for traces of gravity-induced alterations due to 
the propagation of light through different gravitational potentials. The last issue picks up 
the current controversy about an intrinsic anisotropy of space time where the underlying 
data could be explained in terms of cosmological birefringence. 
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5.1 Circular Polarization of AGN 



Active galactic nuclei (AGN) are characterized by a huge variety of high energy phe- 
nomena and nonthermal radiation processes which cannot be ascribed to normal stellar 
activity. Instead it is commonly accepted, that the main central engine is driven by a 
super-massive black hole, whose strong gravitational field in combination with magnetic 
fields and relativistic jets feeds the processes which supply the kind of radiation which 
is central to our current investigation. Of special interest are in this regard the com- 
pact, extragalactic radio sources which belong to the class of active galaxies. Typical 
representatives are Quasars and BL Lac objects, named after the prototype BL-Lacertae 
with strong, variable polarization and weak emission lines (Blazar). A further important 
example is provided by our own galaxy with the central radio source in Sagittarius A*. 
The optical brightness of all these objects is often clearly exceeded by emission in the 
radio regime with a nonthermal spectrum and non-zero degrees of polarization. 

The most plausible, dominant mechanism for such radio emission in AGN is usu- 
ally considered to be synchrotron radiation from an ensemble of relativistic electrons. 
Concerning the often differently used notations it is important to mention, that usually 
the name cyclotron radiation is reserved for emission from low-energy electrons whereas 
synchrotron radiation traditionally describes emission from highly relativistic electrons 
because it was first observed in 1948 in electron synchrotrons. While it is usual to use 
the notion 'cyclotron radiation' with respect to magnetic white dwarfs without regarding 
different energies, the traditional distinction between cyclotron and synchrotron radiation 
is used in most of the relevant literature concerning AGN and therefore also used here. 

Our special interest in this context applies to the origin of circular polarization ob- 
served in only a few AGNs. Analysis of the collation of circular polarization data pre- 
sented by Weiler and De Pater in 1983 |141j revealed that out of 50 radio galaxies and 
43 quasars, circular polarization was detected for only 6 radio galaxies and 15 quasars. 
Although at a first glance it seems very likely that a synchrotron process has to be con- 
sidered as the source, the small degrees of circular polarization observed in many AGN 
has a frequency dependence and time variability that is not consistent with the simple 
predictions based on the intrinsic polarization of synchrotron emission (see 142 also for 
a huge list of references). Likewise the observed degrees of circular polarization are often 
indeed small but nevertheless too high for a simple synchrotron source. It was therefore 
suggested that the circular polarization is due to a propagation effect where the elliptic 
eigenmodes of a relativistic plasma effects cyclic conversion of linear polarization into 
circular polarization |143j . However, this approach is again not free of difficulties so we 
suggest that a conversion mechanism due to gravitational birefringence might solve these 
problems and contribute a significant fraction to the observed circular polarization. 

In the following we give a brief overview on some of the relevant data, mainly fo- 
cussed on observations of the galactic center in Sagittarius A* (Sgr A*) and discuss the 
conventional, propagation induced polarization mechanisms. We outline the idea how 
gravity-induced birefringence could account for the observed polarization. 
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5.1.1 Polarization properties of synchrotron emission 

Following the current understanding, the radio emission of AGN is mainly produced by 
highly relativistic electrons which move on helical paths along magnetic field lines. A 
single electron moving with the synchro-cyclotron frequency u c emitts its radiation in a 
narrow cone of half-angle r y~ l along the instantaneous trajectory, where 7 denotes the 
usual Lorentz factor. For a distant observer this gives a continuous spectrum with a peak 
at the frequency u m oc 7 3 co> c . Considering an isotropic and homogeneous ensemble of rela- 
tivistic electrons within a magnetic field having a power-law energy distribution such that 
the particle density between E and E+dE can be written as N(E) dE oc E~ r dE, one finds 
that the total intensity observed from the source is I{y) oc a(T)(B sin 9p r+1 ^ 2 u~^ r ~ 1 ^ 2 , 
where a(T) is a slowly varying function of T and 9 denotes the pitch angle with respect 
to the magnetic field. The spectral energy distribution, i.e. the dependence of intensity 
I v or flux F v on the frequency can often be approximated by a power law I„ oc v~ a at 
least over limited frequency ranges. The exponent a, called the spectral index is related 
to T by T = 2a + 1. 

In the optically thin case, the polarization plane is perpendicular to the projected 
magnetic field with a linear polarization degree of mi = (r + l)/(r + 7/3). In the 
optically thick regime, the total intensity spectrum reaches a z/ 5//2 dependency, while 
the degree of linear polarization now becomes ttil = 3/(6r + 13), with the electric vector 
being maximum parallel to the projected field. The degree of linear polarization can reach 
theoretically up to ~ 70% but this high level is only rarely observed. Currently, there 
are two accepted mechanisms which can cause depolarization: Unresolved inhomogenities 
within the source and the presence of Faraday rotation internal to the source. The last 
point will be discussed in more detail later. 

Besides linear polarization, the theory of synchrotron radiation predicts that there 
should also be a small amount of circular polarization |148j . Considering again a single 
electron gyrating in a magnetic field B which is directed along k, the circular polarization 
is right-handed when the direction of motion passes close to the line-of-sight (LoS) on 
the opposite side of B , while it is left-handed if it passes close to the LoS on the same 
side as Bo. Therefore, in an ensemble of electrons with an isotropic distribution of pitch 
angles one finds approximately as many electrons contributing to right- and left-handed 
polarization within the emission cone so that, to first order, the circular polarization 
cancels out, leaving the emission linearly polarized. However, Legg & Westfold |148j 
proposed that the total emission from an ensemble of gyrating electrons could have a net 
circular polarization if the number of electrons contributing to the observed radiation with 
right- and left-handed polarization is not equal. For an isotropic distribution of electrons, 
Legg & Westfold showed that this condition is satisfied. The standard equation for the 
degree of circular polarization for an optically thin, homogeneous synchrotron source as 
given by Melrose in 1971 |151j reads 



with 9: angle between the field and the line of sight; v. emission frequency in Hz and 
vb is the electron gyrofrequency. f{o) is a weak function of the spectral index a; for 





sc. 
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optically thick emission in the limit of strong Faraday rotation, f(a) varies monotonically 
between 0.6 and 2.0 for a between and 2 |151j . Using typical source parameters like 
B lmG, sin^ « 1 and v = 1GHz one can expect a level of m c ~ 0.1%. Usually, m c is 
smaller than 0.1% in AGN with only few cases for which CP approaches 0.5% [15Qj. In 
an optically thick source the degree of CP remains similar, only modified by a term of 
order unity dependent on the distribution of relativistic electrons |151j . Equation ()5.1|) 
implies that CP is mainly produced by particles with low Lorentz factor 7, or by emission 
nearly parallel to the magnetic field (small 9). Furthermore, it is important to note that 
CP produced by a synchrotron process is expected to show a v~ x l 2 dependence as can be 
seen from (|5.1|) . 

5.1.2 Observations 

This section will present a brief summary of the most important results concerning cir- 
cular polarization measurements in AGN over the last few years. We do not lay claim to 
completeness and we apologies for omissions. 

Sagittarius A* 

Studies of stellar proper motion in the vicinity of the nonthermal radio source Sagittarius 
A* |144j revealed a highly compact object with a mass of ~ 2.5 x 1O 6 M on a scale less 
than 0.01 pc. Hence, the most conservative interpretation for this is given today by a 
super-massive black hole with a synchrotron emission region, fed through accretion. 

The first detection of CP radiation from Sgr A* was made by Bower, Falke, & Backer 
in 1998 145]. Using the Very Large Array (VLA) they measured an average fractional 
polarization of m c = -0.36% ± 0.05% at 4.8 GHz and of m c = -0.26% ± 0.06% at 8.4 
GHz, respectively. The average spectral index (jn c oc v a ) was a = —0.6 ± 0.3. Since 
the special off-axis design makes the VLA a poor instrument for CP measurements, it 
is important that this result was later confirmed by Sault and Macquart |146j who used 
two archival and one new ATCA (Australia Telescope Compact Array) observation. 

Multifrequency observations of CP and LP in Sagittarius A* were reported recently 
in 2002 by Bower et al. |147j . They obtained 13 epochs of VLA observations at 1.4, 4.8, 
8.4 and 15 GHz in summer 1999 and 11 epochs from ATCA at 4.8 and 8.5 GHz in the 
same year. Their measured mean fractional polarizations are given in the following table, 
together with data taken from the VLA archives. 



Data from 


1.4 GHz 4.8 GHz 8.4 GHz 15 GHz 

(%) (%) (%) (%) 


VLA 1999 
ATCA 1999 
VLA archive 


-0.21 ± 0.10 -0.31 ± 0.13 -0.34 ± 0.18 -0.62 ± 0.26 
-0.37 ± 0.08 -0.27 ±0.10 
-0.31 ± 0.13 -0.36 ± 0.10 



Tab. 5.1: Mean fractional circular polarization from Bower et al. 2002 |147j . 



For the low frequencies they repoted a spectral index a ~ —0.5 while for higher frequen- 
cies the spectrum was best matched by a ~ 1.5. Changes in CP were accompanied by 
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small changes in the total intensity, suggesting that the processes driving these phenom- 
ena are fundamentally linked. Concerning LP, Bower et al. reported no detection with 
an upper limit of 0.2% at 8.4 GHz, less than 0.2% at 22 GHz and less than ~ 1% at 
112 GHz. Aitken et al. [T4*3] have claimed 10% of LP at v > 150GHz, on the basis of 
low-resolution James Clerk Maxwell Telescope (JCMT) observations, but this detection 
is so far not confirmed. 

Quasars and BL Lacs 

A recent ATCA-survey by Rayner et al. |150j for CP in radio-loud Quasars, BL Lacs and 
Radio Galaxies has revealed fractional CP at 5 GHz between 0.05% and 0.5% in 11 out 
of 13 sources at a spatial resolution of 2 arcsec. Also, VLBI measurements from Homan 
& Wardle in 1999 [E2I yielded localized CP of 0.3% up to 1% in the jet-cores of 3C273, 
PKS 0528+134 and 3C279 (1%), while a few cases may be as high as the local linear 
polarization. Brunthaler et al. |153j measured CP in the compact radio jet of the nearby 
spiral galaxy M81*. They reported a value of m c = 0.54% ± 0.06% ± 0.07% at 8.4 GHz 
and m c = 0.27% ±0.06% ±0.07% at 4.8 GHz with an error separation into statistical and 
systematic terms. Similar to Sgr A* they detected no LP im M81* at a level of 0.1%. 

5.1.3 Problem 

Although LP observed in AGN is generally accepted as being produced by a synchrotron 
process, the correct interpretation of the high levels of CP in AGN is still unclear. The 
theoretically predicted degree of LP due to synchrotron emission is about 70% - a level 
which is rarely observed so that several depolarization mechanisms like Faraday rotation 
or strongly tangled magnetic fields have been discussed in order to explain that the 
observed LP is often at least an order of magnitude lower. By taking these various 
mechanisms into account, the low LP is not really a surprise. But, what is a surprise are 
the high levels of CP given these stringent limits on LP. Applying the same depolarization 
mechanisms also to CP, the observed levels should also be an order of magnitude lower 
than the predicted 0.1% from equation ()5.1|) . Additionally, one must take into account 
that the m c oc v~ l l 2 dependence as expected from synchrotron emission is generally 
not observed |142j . For this reasons it is commonly recognized, that simple synchrotron 
models cannot account for the full polarization characteristics without depolarization 
or repolarization, discussed later here, in the source or the accretion region. We will 
therefore briefly recall the circular repolarization mechanism, first invented by Pacholczyk 
in 1970 |154j where CP is explained as being due to a propagation effect. We point out 
some difficulties of this approach and suggest a perhaps more suitable propagation effect: 
Gravitational birefringence. 

Circular repolarization 

For a given polarization mode with Stokes parameters /, Q, U and V it is possible to 
define a second polarization mode with the same / but opposite Q, U and V as an 
equivalent solution of the Maxwell equations. While these so-called orthogonal modes 
can travel independently with the same propagation velocity through empty space and 
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homogeneous isotropic media, in astrophysical magnetized plasmas, however different 
polarization modes have different propagation velocities. In such a plasma for a given 
propagation direction with respect to the magnetic field, one can always find two such 
eigenmodes which travel through this medium with different propagation velocities, but 
without changing their polarization vector. Such a medium with two refractive indices is 
called birefringent. In general, birefringence will be elliptical, i.e. the eigenmodes have 
linear and circular contributions. The transfer equation for a Stokes vector, subjected 
to magnetooptical effects due to different refractive indices in an anisotropic medium 
without absorption is of the form (Melrose & McPhedran 1991 [106 , p. 188) 
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Here, s denotes the distance along the ray path, Ak the difference in wavenumber between 
the eigenmodes and ip the polarization angle. T is the axial ratio of the polarization ellipse 
of one of the eigenmodes. Regarding the physical relevance, we basically distinguish 
between circularly and linearly polarized eigenmodes. 

In the first case, circular polarization implies a polarization ellipse with an axial 
ratio equal to unity, i.e. \T±\ = 1, '+' and '-' denoting the orthogonal modes. Hence, 
Pq — Pu — 0, so that ()5.2j) reduces to 



(5.4) 



Here, the transfer equation leaves / and V unaffected, while U and Q changes only such 
that the polarization angle ip = 1/2 arctanC//Q rotates according 

Ts =- 2 Pv , (5.5) 

which is commonly known as Faraday rotation. Otherwise, if the natural modes are 
linearly polarized this corresponds to T = or T = oo so we get from ()5.2j) 



PuV, ±U = - PQ V 



(5.6) 



This means that if the natural modes are linear or also elliptical then radiation that is 
initially linearly polarized develops a circularly polarized component as the polarization 
changes in a periodic manner along the ray path |155j . Following Kennett & Melrose 
we can therefore write 



V(u) = £/ (^)sin(A 3 RRM), 



(5.7) 



RRM 



3 x 10 4 
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The relativistic rotation measure RRM describes quantitatively the strength of circular 
repolarization, namely the phaseshift A$ between the two natural modes, depending on 
the special details of the distribution of relativistic particles and the local magnetic field 
B. n r denotes the particle density, L the path length within the source region and £l 
the minimum Lorentz factor of the plasma. 

5.1.4 Gravitational birefringence and repolarization 

Although the conditions required for the observed circular polarization of some syn- 
chrotron sources due to circular repolarization are indeed not impossible, they are rather 
very restrictive. We present two situations in which this mechanism might operate, point 
to some difficulties and discuss the possible alternative mechanism of gravitational bire- 
fringence which could circumvent several problems. 

In the first situation where circular repolarization might operate, the relativistic par- 
ticles dominate within the plasma, so that the natural modes are linearly polarized. 
Requiring either RRMA ~ 1 in a small part of the source or RRMA ~ 1CT 3 over a large 
part leads to significant circular polarization, comparable to the degree emitted from 
synchrotron radiation. However, the frequency dependence is more like m c oc u~ 3 than 
m c oc <jj~ x l 2 predicted for a uniform, optically thin synchrotron source |143j . The second 
situation assumes that the plasma is a mixture of cold and relativistic particles causing 
the natural modes to have a small circular component. The relative phase shift A$ 
between the two natural modes is then determined by the cold plasma, and provided 
that A$ <C 1, the resulting circular polarization is of the order of the eccentricity of the 
modes. This situation implies m c oc u^ 1 . 

Even though these scenarios seem to be quite realistic, the existing data do not appear 
consistent with m c oc u~ a with either a = 0.5 or a = 1 [143] . Another problem concerns 
the ratio of linear to circular polarization. For a homogeneous optically thin synchrotron 
source and a highly relativistic plasma one expects m//m c > 100 |158j while Ryle and 
Brodie [157J reported sources with even mi/m c > 1. This problem is also known as the 
"circular polarization excess". 

To summarize up we can say, that despite the different possibilities which have been 
discussed so far, the mechanism for the production of circular polarization in AGNs is 
still not known with absolute certainty. The question therefore arises if gravity- induced 
birefringence could serve as an important source for circular polarization which could 
solve some of the problems, mentioned above. Such a gravitational conversion of linearly 
to circularly polarized light might be possible for rays emitted from the same pointlike 
source and is independent of the special conditions within the plasma, e.g. the amount 
of cold electrons and the optical thickness of the source. However, one obstacle is that 
the birefringence models we have used up to now are appropriate for stars with spher- 
ical symmetry where one can, in principle, assign to each pointsource on the surface a 
heliocentric angle 9 = arccos fi. Consequently, the quality of our limits on birefringence 
strongly depends on the amount of available information about the source. It is clear 
that the situation is now quite different in the case of an active galactic nuclei. Of course 
a first primitive ansatz could consist in putting the nuclei inside an imaginary sphere 
and locate the relativistic jet as the source of linearly polarized light at the limb of this 
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sphere. Using the metric-affine phase shift formula ()1.61|) with the parameters of Sagit- 
tarius A*, e.g. M S grA = 2.5 x 10 6 M , A = 6 x 10~ 2 m and Rs gr A = O.Olpc we get a phase 
shift of A$ = 0.056 x 10~ n k 2 . It is currently not possible to decide, whether this value 
of A$ is sufficient to convert a reasonable amount of linearly into circularly polarized 
light because of the complete lack of a reliable estimate for k 2 of a supermassive black 
hole. Furthermore it is currently unclear how the observed wavelength is correlated to 
the emission radius. A further important aspect which could test for a possible influence 
of birefringence is the frequency dependence of the polarized radiation, described by the 
spectral index. As already mentioned above, the synchrotron model as well as the circu- 
lar repolarization model fail to predict the correct frequency dependence. The hope is, 
that a combined model which incorporates repolarization as well as birefringence might 
be able to provide here a better fit to the observations. 



5.2 Gravitational redshift measurements 

Since electromagnetism was the only fundamental interaction besides gravity that was 
known in 1915 and for which this generalization was relevant, Einstein himself suggested 
the first test of this new principle by showing that an electromagnetic wave propagating 
between points with different gravitational potentials as a consequence must suffer a red 
shift. Therefore, measurements of the gravitational redshift belong to the first classical 
tests of the EEP. 

In every metric theory of gravity, the frequencies of two identical atomic clocks at rest 
at different locations x x and x 2 in a static gravitational potential differ according to 



v\ - v 2 



^ 0o(Xl) 1 , (5.9) 



V2 \l #00 (X 2 ) 

so that for weak gravitational fields, e.g. goo — 1 + 2$ (c = 1), where $ is the Newtonian 
gravitational potential to first order, this yields 

z = $(x x ) - $(x 2 ) . (5.10) 

Since the derivation of ()5.10|) is only based on the validity of the weak equivalence principle 
and the conservation of energy, this prediction is independent of the special form of the 
field equations of a certain theory and therefore valid for every metric theory of gravity. 
Congruously, this result was put in the form of a conjecture by Will 



Universal gravitational redshift conjecture (UGR): 

Any complete, self-consistent, and relativistic theory of gravity that embodies the 
universality of gravitational redshift is neccesarily a metric theory. 



A simple plausibility argument supporting this statement can be given by means of 
the metric postulates (see Sec. 1.2.3, p. 7). Since in a local free-falling reference frame 



5.2. GRAVITATIONAL REDSHIFT MEASUREMENTS 



91 



the physical laws are those of special relativity, the frequencies of free-falling atoms are 
functions only of the universal atomic constants and therefore independent of the external 
gravitational field. Consequently, a comparison of the frequencies of the free-falling atoms 
involves only a comparison of their space-time trajectories, but since these trajectories are 
universal, the measured red shift is universal. Nonmetric theories of gravity violate, per 
definition, one or more of the metric postulates and therefore do not predict a universal 
redshift. On the contrary, Will |159j has shown that in nonmetric theories of gravity the 
ticking rates of atomic clocks in a gravitational field and, so, the gravitational redshift 
are affected in a manner that depends on their internal structure. By applying the 
gravitationally modified Dirac equation in the THefi formalism to the hydrogen atom, 
he found that the predicted redshifts of clocks, defined by transitions between principal 
levels, between fine-structure levels and between two hyperfine levels are given by 



with one clock being far away from the gravitational potential ($(xi) = 0, xi — > oo). In 
order to explore the physical consequences of these predictions, we focus attention in a 
first step on weak gravitational fields ($ < 10~ 5 ) so that the results can be used for tests 
within the solar system ($ ~ 10~ 6 ) and also in good approximation for white dwarfs. 
In the weak field limit, the functions T, H, e and /i are expanded in power series in $ 
according to 



By requiring that a = 1 we achieve that T = 1 — 2$ to first order, so that the theory is 
in agreement with Newtonian gravitation at lowest order. Now, the redshift predictions 
(l5~TTl) - (l5~T3l take the form 
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(5.24) 

(5.25) 
(5.26) 
(5.27) 



The exact expression for the T's and T's are different for different nonmetric theories. 
In metric theories of gravity, all T's and T's vanish identically. In a further project we 
would therefore like to derive the values of these parameters for metric-afhne gravity 
and see how the validity of this theories could be constrained by experiments. Since the 
possible measurable differences in the predictions between universal gravitational redshift 
and those of nonmetric theories becomes more pronounced the stronger the involved 
gravitational fields are, it is obvious to use redshift measurements of solar spectral lines 
or of white dwarf lines for such tests. On the Sun Manganese lines show rich hyperfine 
structure which could be used for setting strong limits on nonmetric theories of gravity. 
Unfortunately in the case of white dwarfs these measurements are seriously complicated 
by extreme pressure broadening which currently renders the signals of finestructure and 
hyperfine structures hardly detectable. 

A possible alternative could be provided by measurements of various redshifts under 
controlled laboratory conditions. Measuring the energy levels of hydrogen in gravita- 
tional potentials different from those on the earth's surface, for example onboard the 
International Space Station (ISS), give the chance of revealing differences in the levels 
of hydrogen to those of standard textbooks. Since the relevant parameters which could 
influence the measurements can be determined with high accuracy, the probability of 
detecting possible deviations is not insignificant. 
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5.3 Further ideas 



The consequences of the idea that gravitational birefringence might have an important 
impact on the interpretation of astronomical polarimetric data are so far a rather un- 
explored topic. We have shown that gravity- induced depolarization has the potential to 
bridge the gap between former AM-Her models and observations and that it also might 
help to solve the puzzle of circular polarization from active galactic nuclei. However, 
since we are involved in a relatively new area of research, many other systems or issues in 
astrophysics may also be better understandable in terms of gravity-induced birefringence 
or could, at least, provide strong upper limits on this effect. We therefore provide here 
two brief examples of such systems which it might be worth to turn one's attention to 
and which we have not yet explored as deeply as the foregoing topics in this chapter. 

Gravitational lensing: In a gravitational lense system, light emitted by a distant 
astronomical object is deflected by the gravitational field of a massive object (star, galaxy, 
cluster of galaxies) which lies along the line of sight, allowing for multiple source images. 
Depending on the relative orientations of source, lens and observer multiple source images 
correspond to different trajectories of the light around the lense. Consequently, light rays 
of different images could have propagated through different gravitational potentials, so 
that the polarizations of these images becomes different if gravitational birefingence is 
not negligible, taking into account the time delay. A systematic survey of the polarization 
properties in gravitational lens systems therefore not only has the potential for strong 
upper limits on this effect but also for a direct detection of birefringence in principle. 

Polarization rotation over cosmological distances: We have already shown in this 
chapter that gravity-induced birefringence might help to understand the polarization 
properties of active galactic nuclei. Since many of these objects are at high redshifts, 
the polarization properties of light could therefore also be altered by propagating over 
cosmological distances in a nonmetric background gravitational field. 

Recently, several authors claimed to have found a systematic rotation of the plane of 
polarization of light emitted by distant radio galaxies [16"H| llfil] [162 . By measuring the 
angle 6(X) between a fixed reference direction in space and the plane of polarization of a 
radio wave of wavelength A, they found that their data could be fitted by 



The linear dependence of the angle 9 on A 2 is a characteristic feature of Faraday rotation. 
The fitting parameter a depends upon the magnetic field and the electron density along 
the line of sight. Therefore, the angle x describes the orientation of the polarization 
plane after Faraday rotation is taken out, i.e. the orientation of the plane before the 
rotation. Now, Nodland and Ralston 162, claimed that the observed x angles could be 
reproduced by assuming that the polarization plane of a wave emitted by a galaxy is 
initially oriented at a fixed angle relative to the galaxy's major axis, and then undergoes 
a rotation, specified by an angle /3, that depends on the direction of the galaxy in the 
sky. Defining an angle 7 between the line of sight to the galaxy and a fixed direction is 



0(A) = a\ 2 + x 
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space they found 

/? = -A~VcoS7 , (5.29) 

where r is the distance to the galaxy and A denotes a rotation measure. This dependency 
of (3 on a fixed direction in space, provided by 7, is an indication of anisotropy. The 
importance of this effect lies in the fact that it can probe the fundamental structure of 
space time at very early times and over length scales commensurate with the size of the 
universe. Because of the clear importance that the confirmation of an intrinsic anisotropy 
of space would have, much effort have been put into explaining these results. At least 
three mechanisms have been proposed which can lead to this cosmological birefringence 

• Modification of conventional electrodynamics by introducing pseudoscalar axion 
fields [IH3]. 

• Modified dispersive Maxwell equations due to semi-classical space time with polymer- 
like structure at microscales in loop quantum gravity [164 . 

• Rotation of the polarization plane due to a nonmetric gravitational background 
field. 

The polarization rotation therefore serves as a test for alterations of various basic phys- 
ical theories which should be thoroughly studied to analyse the distinguishing features 
and measurable differences between them with the hope of finding true modifications of 
fundamental physics. A high priority in this list of possible alterations should be given 
to the coupling of torsion to electromagnetism which appears natural in the framework 
of metric-afline gauge theories of gravity. The main advantage of this approach is that 
its predictions are already testable and even perhaps verifiable in other, more accessible 
branches of astrophysics, as was presented in this thesis. 

It is important to note that the importance of further tests of these different modifi- 
cations is unaffected by the current controversy over the correctness of the interpretation 
of the data |165| 11661 H67j . Even if Nodland and Ralston |162j could be proven wrong, 
these studies can be useful to uncover the source of any polarization rotation effect in 
the future. 



Chapter 6 
Conclusions 



The main purpose of this thesis was the development of new astronomical tests of the 
Einstein equivalence principle (EEP) in terms of gravity-induced birefringence. So, the 
very first question we had to answer in the beginning was how justified it is to proceed 
with testing a principle which belongs to the physical predictions with the most accurate 
empirical underpinning. However, this principle was invented at a time when physicists 
discussed about the existence of atoms and, so, it is obvious that the Einstein and also the 
Weak equivalence principle are blind to the microscopic structure of matter. Regarding 
all kinds of charges and spins in modern quantum field theory it would be more than 
remarkable if this statement holds up for all times. Therefore the desired unification of 
quantum mechanics with general relativity provides the strongest clue for violations of 
the Einstein equivalence principle as was shown in chapter one. 

An appropiate framework for developing and analysing experimental tests of the EEP 
is given by lagrangian based nonmetric theories of gravity. Within this class we have 
focused our attention on two prototypical representatives, the nonsymmetric gravitation 
theory (NGT) and the metric-affme gauge theory of gravity (MAG), both build upon a 
non-Riemannian geometry of space time. While it was known for several years that NGT 
predicts violations of EEP in terms of gravitational birefringence where the strength of 
birefringence is detemined by a material dependent coupling constant £ 2 we have shown 
that this is also the case for MAG. The latter couples torsion to the electromagnetic 
field in the \9~ formalism, giving rise to a new coupling constant k 2 . In this context it is 
important to note that the question of how the electromagnetic field exactly couples to 
gravity, in particular to torsion, is an age-old, still unresolved problem. Certainly, there 
may be other couplings than those we have proposed, which might seem no less natural 
so that a lot of work is left for the future. However, our objective of constraining possible 
EEP violations in these theories is therefore equivalent with setting strong limits on I 2 
and k 2 . 

The search for traces of gravity- induced birefringence in polarized light from an as- 
tronomical object is promising only if we have a sufficient knowledge about its source 
properties so that we can compare the observed polarization signal with its theoretically 
predicted 'original' properties without any birefringence influence. For this reason we 
have used in our first project polarization measurements in solar spectral lines where 
the underlying theory about the generation of polarization profiles within the solar pho- 
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tosphere is well developed. Using the Stokes asymmetry technique and the new profile 
difference technique, our new limits on £q are 3-7 orders of magnitude smaller than pre- 
vious results. Since birefringence is most pronounced for short wavelengths these limits 
could be improved with the same methods by using UV polarization measurements where 
we can expect a maximum gain relative to the current analysis of a factor 4. Currently, 
however, it is not possible with our methods to decide whether the observed asymme- 
tries of solar Stokes profiles are partially due to gravitational birefringence since possible 
influences must be less than or equal to the asymmetries induced by nongravitational 
mechanisms. 

At this point a very important issue must be mentioned: Conventional models which 
try to explain the creation of Stokes profiles in the solar atmosphere are not axiomatic 
models from first principles. Instead the main guideline is the agreement with observa- 
tions to which the model has to be fitted. Consequently, the predictions of these models 
which are based on interpretations of polarized signals must not be valid if gravitational 
birefringence has a non negligible influence. 

The limits on k 2 were the first so far made and therefore could not be compared with 
other, older values. Although the idea of a possible utilization of future space missions 
is very appealing because of its high potential for further improvements of our current 
limits, the technological realization of such an experiment is currently out of reach. 

In our second project we have used polarization measurements from isolated magnetic 
white dwarfs to constrain birefringence. By using a dipole model for the magnetic field 
geometry we were able to sharpen previous results, given by Solanki, Haugan and Mann 
for GRW +70°8247 by ~ 18%. In addition to GRW we found three other white dwarfs 
which complied with our restriction of available wavelength resolved polarimetric data 
and also of known mass, radius and magnetic field geometry. The resulting limits indicate, 
that k 2 for a particular star only depends on the ratio between its Schwarzschild radius 
and the physical, stellar radius. Whether k 2 also depends on the chemical composition of 
a celestial body like £ 2 cannot be decided conclusively. The source of the torsion field is 
not specified in most papers concerning metric-affine gravity and, therefore, they do not 
make firm predictions of differences between the torsion fields generated by, for example, 
different kinds of stars. It is certainly possible that such differences exist. For this reason 
it would be interesting to compare limits on k 2 obtained from a large sample of stars with 
different chemical abundances. 

The polarization characteristics of AM-Her systems provide a fundamentally different 
analysis for our purpose. By calculating the gravitationally modified polarization curves 
for VV Puppis we could achieve a good agreement with observations only for a nonzero 
k 2 . However, since other models which assume different source properties are able to 
achieve the same good conformity with observations, gravitational birefringence can only 
serve currently as one possibility among other more conventional approaches. 

This thesis represents the first systematic, extensive search for gravitational birefrin- 
gence in astrophysical spectropolarimetric data. Therefore several promising projects are 
left for further investigations. In this regard, one the most interesting open question is 
concerned with the origin of circular polarization from active galactic nuclei. Gravity- 
induced birefringence seems to be the ideal candidate since its influence is not afflicted 
by entangled magnetic fields or different compositions with respect to the amount of 
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relativistic or cool electrons in hot plasmas, which are the main sources of problems for 
the conventional repolarization approach. Nevertheless, red shift measurements, gravita- 
tional lenses and the question of intrinsic cosmic anisotropy have also the potential for 
making major progress in this field of research. 

The quest for possible violations of the Einstein equivalence principle still remains an 
open problem. Although we have several theoretically convincing arguments as well as 
experimentally gained hints, the results are so far not conclusive. However, the results of 
this thesis have shown that this possibly could change, not too far away in the future. 
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Appendix A 

Gravity-induced birefringence in the 
Xg-formalism 



In this chapter the propagation of polarized light is investigated under the assumption 
of a relatively weak gravitational background field, that vary on length and time scales 
which are large compared to the light's wavelength. Since the electromagnetic lagrangian 
density of special relativity has the form (jl.l2|) with x a ^ lS = ^(rf^rf 5 — rj^r]^ 1 ) it is 
possible to find a quasi-Lorentzian coordinate system in the weak-field limit where x af3lS 
has the form 

x *f3js = \{jfnjfs _ ^Y 7 ) + Sx a ^ 5 (A.l) 

with 5x aPlS < 1. 

The amplitude and phase representation of a plane, electromagnetic wave is given by 

E = A £ e i *, B = A B e** . (A.2) 

The main feature of this representation is that the temporal and spatial derivatives of the 
vector amplitudes Ag and A# are small compared to derivatives of the rapidly varying 
phasefunction $. Since the gravitational background field varies only slowly in space and 
time too, its derivatives are also small compared to derivatives of By considering the 
propagation of a high-frequency electromagnetic wave this means that all derivatives other 
than those of the phasefunction can be ignored. Therefore, locally the gravitational field is 
treated as a homogeneous medium, which is known as the high-frequency approximation. 
So, the electromagnetic field equations which follow from (|1.12|) are 

X a07S F lS ,0 = O . (A3) 

Defining electric and magnetic fields in the usual way as F 0i = Ei and = e^Bi 
together with the decomposition (jA.ljl one gets the modified Maxwell equations describing 
electromagnetic fields in a background gravitational field. Written in a somehow more 
transparent way this reads as 

V • E + terms proportional to 8\ and E or B (A. 4) 

and 

<9E 

V x B — - — h termes proportional to 5x and E or B . (A. 5) 
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The homogeneous Maxwell equations are unaltered 

V x E + — = (A.6) 
at 

and 

V x B = . (A.7) 

The objective target of the following analysis is the derivation of the eikonal equation 
which describes the propagation of a locally plane electromagnetic wave and, therefore, 
provides information about the local coordinate velocities of wave propagation. For this 
purpose (|A.2J) is inserted into the field equations (|A.4jl - (|A.7|) neglegting all derivatives 
other than those of the phase function $. Denoting k^ as the gradient of this function 

k„ = d$ = (d$/dt, V$) = {-uj, k) (A.8) 

the homogeneous Maxwell equations yield 

k x A E , , , 

A B = (A.9) 

UJ 

and 

k-A B = (A.10) 

respectively. While the latter of these implies that the magnetic field vector of a locally 
plane wave is transverse to the direction in which the wave propagates, it follows from 
(TPl) that 

k ■ A E = Terme proportional zu 5x und A E oder A B (A. 11) 

i.e. electric field vector is transverse to the direction of propagation if and only if no 
gravitational field is present! Together with equation (|A.9j) this implies for A E expressed 
via the two independent components of A B 

A E = — — k x A B + terms proportional to 5x and A B (A. 12) 

k 2 



so that this equation together with (jA.5|) and (jA.9|) implies the eikonal equation 

( 1 — — ] A B = terms proportional to 5\ and . (A.13) 
V k J 

Since the magnetic amplitude A B can be regarded as the superposition of two inde- 
pendent polarized components with the coordinate phase velocity ui/k, finding the two 
independent polarization states is equivalent to solve a two-dimensional eigenvalue prob- 
lem. For this purpose one makes a decomposition of 5x a ^ & in a set of SO (3) tensor 
objects 

= -5 X 0i0j (A. 14) 

rfj = _ e jlm 5x 0ilm ( Al5 ) 

( ij = - e ilm e jpq 5x lmpq (A.16) 
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where e ljfe is the Levi-Civita antisymmetric symbol. This decomposition of 5x al3 ^ s is 
now rotated from the original, quasi-Lorentzian coordinate system (t, x, y, z) into a new 
system where the background gravitational field is represented to a set of (t, x', y', z') 
coordinates and the light propagates in the z'-direction. In this representation A# has 
only x' and y' components so that equation (jA.13|) reduces to a system of two equations 

AA l B -BA 2 B (A. 17) 

-BA l B +CA 2 B . (A.18) 

The matrix which defines the structure of the right-hand sides of (jA.17|) and (jA.18|) is 
real valued and symmetric when 5\ a ^ S is real. The coefficients A, B and C depend on 
the location in space-time and on the direction in which the wave propagates. This can 
be expressed in terms of the tensor components ^ , Q J and 7* in the (t, x', y', z') 
coordinate system. In particular Haugan & Kauffmann [20] has shown 

A = f 2 ' - 27 2 ' 1 ' - C 1 ' 1 ' (A.19) 
B = ^' 2 ' + ( 7 2 ' 2 '-y' 1 ') + C 1 ' 2 ' (A.20) 

C = e' 1 ' + 27 1 ' 2 ' - C 2 ' 2 ' • (A.21) 
The eigenvalues of (JA.17J) and (jA.18jl are given by 

\± = ^±y(A-Cy + 4B 2 . (A.22) 

From this it follows that the corresponding eigenvectors define the polarization states 
which propagate with well-defined phase velocities 

c± = l-\\± + 0(5 X 2 ) (A.23) 

Denoting the fractional difference between c + and c_ by 5c/ c one gets a local, dimen- 
sionless observable 

^ = y(A-Cy + 4B* . (A.24) 

Therefore, during the propagation of wave with circular frequency to this effects yields 
an accumulated phase shift 

A$ = lu J ^ dt + 0(5 X 2 ) (A.25) 

e.g. the difference in the local velocities leads to a change in the relative phase between the 
two independent components which in turn implies an alteration of the initial polarization 
state Since this analysis revealed that the speed of an electromagnetic wave depends on 
its orientation within the background gravitational field this implies that gravity-induced 
birefringence is a direct consequence of a violation of EEP. 
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This scheme is now applied to the special case of metric-affine gravity. As explained in 
the introduction the nonmetricity independent, tensorial part of the torsion as spherically 
symmetric solution of the metric-affine field equations was given by Tresguerres in 1995 

[52 EB] 



T a = k 



0i) + ( m Ar gftjffi 

\ r 2 3^ xr 3 



A9 C 



(A.26) 



with the dilatation charge N , the torsion mass w! and ko = 1 for a = 0, 1 and ko = —1/2 
for a = 2, 3. Setting 

/ m Ar nb^N^ 



A 



3X 



(A.27) 



this gives 



rp3 



- 9° A 9 1 — A9° A 9 1 

r 

-9° A9 1 + A9° A9 1 
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1 

~2 
1 

~2 
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"2 
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A) 9° A9 1 

r J 

- + A) 9° A 9 1 

r ) 



-(9 U A 9 

f \ 

1 



2 A 2 ) + A (6° A 9 2 + 9 1 A 9 2 
1 /l 



+ A ) 9° A 9 



r 



- A )9 1 A9 2 



(e° A 9 3 - 9 l A 9 3 ) + A {9° A 9 3 + 9 1 A 9 3 



+ A) 9 A 9 + 



A) 9 A9 



(A.28) 
(A.29) 
(A.30) 
(A.31) 
(A.32) 
(A.33) 



which can be expressed in terms of the most general structure of the torsion in the case 
of spherical symmetry [37] 



T 1 

rp3 



a(r) 9° A9 1 
(3{r) 9° A 9 1 

7 (1) 9° A 9 2 + 7 ( 3) 9 1 A 9 2 
7(i) 9°A9 3 + 1{3) 9 1 A 9 3 



(A.34) 
(A.35) 
(A.36) 
(A.37) 



with 7(2) = 7(4) = 0. According to (jl.59|) and equation (3.11) in jHJJ together with the 
usual scalar-valued two-form F, representing the electromagnetic field one can write 



T AF 
-T x A F 
-T 2 A F 



a 9° A 9 1 A F = a F 23 9° A 9 1 A 9 2 A 9 3 

(39° A9 1 AF = (3F 23 9° A9 1 A9 2 A9 3 

7 ( i) 9° A 9 2 A F + 7( 2 ) 9° A 9 3 A F + 

70) 9 1 A 9 2 A F + 7(4) 9 1 A 9 3 A F 

7(i) Fi 3 9° A 9 2 A 9 1 A 9 3 + 7(2) F 12 9° A 9 3 A 9 1 A 9 2 + 

7(3) F 03 9 1 A 9 2 A 9° A 9 3 + 7(4) F 02 9 1 A 9 3 A 9° A 9 2 

[7(2) F 12 - 7(1 ) F13 + 7(3) Fos ~ 7(4) ^02] 9° A 9 l A 9 2 A 9 3 

[7(2) B 3 + 7 ( i) B 2 + 7(3) B 3 - 7(4) E 2 ] 9° A 9 1 A 9 2 A 9 3 
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— T 3 A F = -7(2)0° A9 2 AF + 7(1) ^° A # 3 AF 

-7(4) 9 1 Ad 2 A F + 7(3) A9 3 A F 
= -7(2) Fi 3 6° A6 2 A9 1 A 6 3 + 7(1) F 12 6° A 6 3 A 6 1 A d 2 

-7(4) ^03 1 A6 2 A9° A 6 3 + 7(3) F 02 6 1 A 9 3 A 6° A 6 2 
= [7(2) Fi 3 + 7(i) ^12 - 7(4) ^03 - 7(3) ^02] 0° A 9 1 A 9 2 A 6 3 
= [-7(2) B 2 + 7(1) 5 3 - 7(4) E 3 - 7(3) £ 2 ] 6° A 1 A Q 2 A 9 3 . 

Since the lagrangian density we refer to reads 

5C EM = k 2 *(T a A F)*(T a A F) (A.38) 

this leads to 

*(T a A F)*(T a A F) = [*T° A F] 2 — [*T { A F} 2 

= {a Br) 2 - (f3 B 1 ) 2 - ( 7(2) B 3 + 7(1) B 2 + 7(3) E 3 - 7(4) £ 2 ) 2 - 
(7(i) B 3 - 7 (2) -B 2 - 7(4) £3 - 7(3) -E2) 2 

= (a 2 - /3 2 )-B 2 - ( 7 (2) Bl + 27(4)7(2) £2^3 + 27 (2) 7( 3 ) B 3 E 3 - 
2 7(2)7(4) B 3 E 2 + 7( 2 X ) -B 2 + 27(4)7(3) B 2 E 3 - 27(4)7(4) 5 2 E 2 + 
7(3) #3 - 2 7( 3 )7(4) £ 2 £ 3 + 7 (4) £ 2 ) - (7 (1) B 3 - 27(i)7( 2 ) S 2 S 3 - 
27(1)7(4) B 3 E 3 - 27(1)7(3) 5 3 E 2 + 7 ( 2 2) B\ + 2 7 ( 2 ) 7 ( 4 ) 5 2 E 3 + 
27(2)7(3) 5 2 E 2 + 7( 2 4) El + 27(3)7( 4) £ 2 £ 3 + 7( 2 3) E\) 

= (a 2 - f)B\ - ( 7( 2 i) + 7(2))£ 3 2 + 2(7(1)7(4) - 
7(2)7(3))5 3j B 3 + 2(7(1)7(3) + 7(2)7(4))5 3 £ 2 - 
(7( 2 i) + 7( 2 2 ))^2 2 - 2 (7(i)7(3) + 7(2)7(4) )B 2 E 3 + 
2(7(1)7(4) - 7(2)7(3) )B 2 E 2 - ( 7 2 3) + 7 ( 2 4 ))^ 3 2 ~ (7( 2 3 ) + 7( 2 4 ))^ 2 2 

Rearranging the terms yields 

*[T a A F]*[T a A F] = (a 2 -m 2 -(7( 2 i)+7( 2 2))[52 2 + 53 2 ]-(7( 2 3)+ 7( 2 4)) 

[E 2 2 + E 2 3 ] + 2(7(1)7(4) - 7(2)7(3) )[B 2 E 2 + B 3 E 3 ] + 
2(7(1)7(3) + 7(2)7(4)) [B 3 E 2 - B 2 E 3 ] 

which can be written in terms of matrix elements 

*{T a A F)*(T a A F) = (nB 2 l +( 22 B 2 2 +( 33 B 2 3 -& 2 E 2 2 -& 3 E 2 3 +2 l22 E 2 B 2 + 

2733 E 3 B 3 + 27 32 E 3 B 2 + 2'j 23 E 2 B 3 . 



The symmetries of 5x a/3,y5 implies that ^ and £ u are symmetric. 
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The explicit matrices are 



( (« 2 - P 2 ) 



/ 

(7 2 3)+ 7 ( 2 4)) 
\ (7 2 3)+7 ( 2 4)) 






,2 i „,2 










-(7(i)+7( 2 2) ) 

-(7( 2 i)+7( 2 2) ) 






„2 i „,2 



(A.39) 



(A.40) 



and 



\ 

' " | (7(1)7(4) - 7(2)7(3)) (7(1)7(3) +7(2)7(4)) ( A - 41 ) 
-(7(1)7(3) +7(2)7(4)) (7(1)7(4) -7(2)7(3)) / 

The expressions A — C and B from (|A.22|) can now be expressed in terms of the spherical 
components of the tensors C? 3 and 7 y . In terms of components in the (t,x',y',z r ) 
coordinate system Haugan & Kauffmann [201 have shown that 



A-C 



((^ + ei) + «(7^ - 7^) + (cr + a,)) 



(2) 



.(2) 



.(2) 



and 



-^W#-a0+2(7? ) +75)+^ 2) -C5)) 



(A.42) 



(A.43) 



Only I = 2 components appear. Expressions for A — C and B in terms of components in 
the original (t, x, y, z) coordinate system follow from the transformation law for spherical 
tensor components [156J, e.g. 



(A.44) 



where 0, 9 and if) are the Euler angles specifying the rotation from (t, x, y, z) to (t, x', y', z') 
and the rotation matrix D given in terms of spherical harmonics. It is now useful to 
introduce a local quasi-Lorentzian (t, x, y, z) coordinate system at each point of the light 
ray's path oriented so that the x axis lies in the ray's plane because the spherical tensors 
introduced above are simple in these local coordinate systems. Specifically, (Iff and 
jiff are nonzero only for m — 0. 

At each of these points along the ray path, the local (t, x, y, z) coordinate system 
is rotated about the y axis through an angle 6 to obtain a local (t, x', y', z') system so 
that now the ray runs in the z' direction. The local value of 5c/ 'c in (jA.24|) is expressed 
in terms of A — C and B which are, in turn, according to ()A.42j) and (|A.43j) expressed 

(2) (2) (2) 

in terms of £Y_2> C±2 an d 7±2- Since the Euler angles of the rotation from (t,x,y,z) to 
(t, x', y', z') are 9 and <ft = ip = 0, the transformation law ()A.44j) implies 



±2 



sin 2 ^ 



(A.45) 
Again, 

according to Haugan & Kauffmann the form of (jA.42|) and (jA.43|) and the transformation 



with the same relationship between an d Cq 2 ^ 



and between 7^2 and . 
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law (|A.44|) implies that B is proportional to 7 2 ' ) while A — C is proportional to +(q 2 ^)- 
In order to get explicit values for 5c/ c one can express Q , £o 2 \ Co "*) Co 2 "* an d To , To 
by the matrix elements and 7 l - ? 



t(0) 
SO 


_ til 
— s 


, ^22 , £33 _ 

+ s + s — 


9(% 2 , ..2 \ 
47(3) + 7(4) J 


t(2) 
SO 


= e 1 


-^ 22 +e 33 : 


> = "(7(3) +7(4)) 


a(0) 

so 


= c 11 


+ c 22 + c 33 = 


(a 2 -/3 2 )-2( 7( 2 1) +7( 2 2)) 


a(2) 
SO 


= c 11 


_1 (C 22 + C 33 


) = (a 2 - (3 2 ) + (7 ( 2 1} + 7( 2 2) 


(0) 

7o 


= 7 11 


+ 7 22 + 7 33 = 


2(7(1)7(3) - 7(2)7(4)) 


(2) 

7o 


= 7 11 


-^(7 22 + 7 33 


) = -(7(1)7(4) -7(2)7(3)) 



Since, according to (|A.36|A.37j) we have 7(2) = 7(4) = and, therefore, B = 0. So, in this 
case we have 

^-C = -^(d 2) + Cf)sin 2 ^ . (A.46) 
Using the expressions (|A.28J) - (|A.33J) one gets 

(2) 1 / 1 2A a9 \ . , t . 

(o 2) = -- + - A {^ + - + A ) (A.48) 

which yields 

f« +rf ) = _^ + ^ = _M = _^ . (A ,49) 

(y (y ry ryO 

In this equations we have dropped the last two terms in (|A.27|) since the dilatation charge 
N vanishes if the nonmetricity field does and the small observed value of the cosmological 
constant A means that effects of its term can be neglected on galactic and smaller scales. 
Together with (jA.26|) and k 2 from (jA.38|) this yields 

A-C = = sin 2 # . (A.50) 

V6r 3 

Therefore the fractional difference between the velocities of the two polarization states is 
given by 

5c /-k 2 m 9 „ . , _ . 

— = S in 2 . A.51 

c r d 

For the total phase shift A<3> which accumulates between the source and the observer one 
now has to calculate 

5c , 1- , o f sin 2 6 



c J r 3 



dt = -V6uk 2 m / — — dt . (A.52) 
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Using the definitions explained in Fig.Bl with the ray parametrization x(t) = b + k t 
and sin 2 #(t) = b 2 /R 2 {t) (R(t) = r), we can write 



A^ = -VQ^k 2 mRl(l-f, 2 ) J —t- 



dt 



to=-RoA< 



(i^(l-^))6/2 



(A.53) 



where we have used R(t) = \/b 2 + t 2 and b = R(t)\J~(l — fi 2 ) with /i = cos9(t). R denotes 
the radius of the star. 

The integral could be solved with the substitution 



R = at 2 + bt 2 + c, a = l,b = 0,c = Rl(l-fi 2 ) . 



(A.54) 



So that we have 



h = 



At 



oo 



dt 

Jtf = 3(4^(1 _ ^ ))(jR 2 (1 _ ^ + t 2)3/2 



12i2§(l- ^) J 3 



^3 4i^(l - A* 2 )(i2g(l - A* 2 ) + t 2 ) 1 / 2 



Evaluating the above term for J 5 yields 



t =RoH 



R 2 (l + fi) 



+ 



5 " 3i? 4 (l-^) 1 3i? 4 (l-^)(l + /x) 
so that the accumulated phase shift becomes 



(A.55) 



■\/~6ujk 2 muR^il — /x 2 )/ 5 



' 2 k 2 muj 
3 Rq 



H- 



i + u 



Finally we get 



>2 2n k 2 m /( / u + 2)( J u- 1)' 



(A.56) 



' 3 XRl \ n + 1 J 

Whether k 2 also depends on the chemical composition of a celestial body like £ 2 cannot 
be decided conclusively. The source of the torsion field is not specified in most papers 
concerning metric-affine gravity and, therefore, they do not make firm predictions of 
differences between the torsion fields generated by, for example, different kinds of stars. 
It is certainly possible that such differences exist. 
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Tests using the Solar Probe spacecraft 




In this chapter we derive the formula which gives the accumulated phase shift between a 
point, located at a distance R from the sun and an observer in an infinite distance. The 
principles of this calculation are mainly based on the calculation that Gabriel et al. 
used to derive the phase shift formula ([1.57)1 . 

We consider a ray which originate on the Sun's surface and use the ray parametrization 
x(t) = b + kot. As usual the unit vector ko, for convenience not explicitly shown in the 
figure specifies the direction of the (unperturbed) ray. By demanding that k • b = we 
define b as the impact vector which connects the center of the sun with the closest point 
on the ray. Those points on the ray which are inside the Sun for the case b < R are 
irrelevant in the calculation. The integration 

1 °° 

A$(ji) = ~u J fisin 2 6{t)dt (B.l) 

begins at the Sun's surface, t = (R% — 6 2 ) 1 / 2 = R cos9 = Rqh with b = R(t) sin#(£). 
The integration extends to an infinite observer in t± = oo. For convenience we use for 
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t- 4 

the speed of light c = 1 in our calculations. Using the NGT relation Vt = -^^y this yields 



A* M = i. / V*? = % J (R ^-^ r * • (B.2) 

to=Roii to=RofJ. 

This integral could be solved with the substitution 

R = at 2 + bt 2 + c, a = l,b = 0,c = R 2 (l-n 2 ) . (B.3) 

So that we have 

dt 2at + b f 1 3a \ Qa 2 ( 2 2at + b\ 

—5 = 1 — 9 H = + -ttt — f= arctan (B.4) 

if A \ 2 R 2 A-RJ A 2 Vv/A v^A y 

with 

A = 4ac - 6 2 (B.5) 

This solution is valid as long as A > 0, i.e. 4i?o(l — f 2 ) > which is always fulfilled. 
Performing the integration in (jB.4j) yields 



dt 2t ( 1 



, if 4flg(l-/*a) V2(^(l-^ 2 )+t 2 ) 2 

to=fMKO 



+ (B.6) 



3 \ 6 

+ 



4i^(l-^)(i2g(l-^ 2 )+t 2 ); 16^(1-^ 



2\2 



1 t 

arctan ■ 



37T Rq[J, 1 1 3 



16^(1 - /i 2 ) 5 / 2 2i? 2 (l - /1 2 ) V,2i2g 4i2g(l - ^ 2 

6 u 
arctan ■ 



16^(1 -/i 2 ) 5 / 2 (l-/i 2 )V2 

37T yU 



16^(1 -/i 2 ) 5 /2 ARlil-n 2 ) 
3/i 3 



8^(1 -/i 2 ) 2 8^(l-/i 2 ) 5 /2 



arcsm /i 



so that we get the usual function for the phase shift between the solar surface and an 
observer in infinite distance 



A^W_ = ^ 2 (l-/i 2 ) / I (B.7) 



t =nR 



7T^q / 37T /i 3/i 



A.Rjj V 16 ( 1 _ /^ 2 ) 3/2 4 8(1 -/x 2 ) 
3 . \ 



8 (1 -^2)3/2 



arcsm fj, 
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The phase shift which is accumulated between the solar surface and R is given by 

A*(a*W-^ = ^(l-^ 2 ) / f (B.8) 



t =/iR ^' 

ir£% (R ( 1 



A \2 \2{R 2 {1 - ^) + Riy 
3 



4J%(1 - fi*){I%{l - ft) + &) . 



3 R a 
o\n arctan 



l-fi 2 ) 2 R 3 i? (l - /x 2 ) 1/2 4i?g 

3/i 3 



8(1 - /i 2 )i?jj 8jR 3(! _ ^2)3/2 



arcsm \i 



One can see that, if the signal travels in the opposite direction from R towards to = fiRo 
we have to switch the integration limits which gives a negative sign for the phase shift. 

At last, we get the phase shift between R and the infinite distant observer by sub- 
stracting (IB.8jl from (IB.7jl . 



I I °° 7 ^7 

7ti% 0/ 0n \ r at r at 



A$(/i)^oo = ^^(l-/i 2 ) J =3 - y =3^ (B.9) 



tt£% ( 3vr i? 



A V 16(1 - /i 2 ) 3 / 2 i?jj 2 \,2(i?g(l-/i 2 ) +i? 2 ) ; 



4i? 2 (l - /i 2 )(^(l - /i 2 ) + i? 2 ) ; 8R 3 (l - ^ 2 ) 3/2 



arctan 



Ro(l-^ 2 , 



Since for our purposes the best spacecraft position is given when it passes the solar limb, 
we can focus on the case \i = which yields 

A* I = <(^L-*( ! + ? ) (BIO) 

' £?„ A \l6Rl 2 ^2(i? 2 + i? 2 ) 2 + ARl{Rl + R 2 ) J 1 ' Uj 

3 fl\ 
■ arctan ■ 



8R 3 Ro, 
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